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ABSTRACT 

The  Electromagnetic  Vorticity  Meter  (EMVM)  directly  measures  small-scale  vorti- 
city  in  the  ocean  using  the  principles  of  electromagnetic  induction.  The  scientific 
motivation  for  developing  the  EMVM,  the  design  concept,  and  the  theory  of  operation 
are  discussed.  Complete  descriptions  are  given  of  the  instrument’s  design,  construction, 
and  performance  through  several  evolutions  of  the  sensor,  and  extensive  observations  are 
presented  of  sensor  performance  during  laboratory  and  field  tests.  The  tests  demonstrate 
a  vorticity  noise  level  of  10~^  s”^ . 
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1.  INTRODUCTION 

An  instrument  and  technique  have  been  developed  for  the  direct  measurement  of 
small-scale  ocean  vorticity.  This  report  reviews  the  theory  on  which  the  instrument, 
called  the  Electromagnetic  Vorticity  Meter,  or  EMVM,  is  based,  describes  its  evolution 
and  testing  from  1982  to  the  present,  and  presents  results  obtained  during  laboratory  and 
field  tests. 

Section  1  defines  vorticity  and  discusses  its  importance  in  oceanography.  Section  2 
describes  the  concept  and  operating  principles  of  the  EMVM,  including  detailed  equa¬ 
tions  explaining  its  operation.  Section  3  contains  detailed  descriptions  of  the  first  four 
versions  of  the  instrument,  and  Section  4  describes  the  fifth  version.  Finally,  the  results 
of  vorticity  observations  conducted  with  the  EMVM  in  Puget  Sound  are  given  in  Sec¬ 
tion  5. 

1.1  What  is  Vorticity? 

Vorticity  is  a  property  of  the  physical  state  of  a  fluid.  Qualitatively,  it  is  the  spin  or 
swirl  of  a  fluid  parcel.  Mathematically,  it  is  the  curl  of  velocity,  V  x  v.  Measuring  the 
full  vector  vorticity  is  a  formidable  task.  The  curl  is  composed  of  nine  derivatives, 
although  only  six  are  independent  for  divergence-free  flow. 

Vorticity  is  fundamentally  important  in  fluid  mechanics  and  particularly  important 
in  physical  oceanography.  Because  few  processes  generate  or  dissipate  it,  once  there,  it 
can  be  long  lived.  In  fact,  some  theories  treat  vorticity  as  a  conserved  quantity  that  does 
not  change  as  the  fluid  moves.  In  addition,  because  of  its  organized  motion  (i.e.,  swirl),  it 
stands  out  from  the  general  background  motions.  Often  the  description  of  a  particular 
flow  feature  is  most  simply  stated  in  terms  of  the  vorticity  distribution  and  dynamics. 

In  the  open,  undisturbed  ocean,  only  boundary  friction  from  wind  at  the  surface  or 
drag  at  the  bottom  and  mixing  (often  caused  by  breaking  internal  waves)  can  impart  vor¬ 
ticity.  This  characteristic  is  evident  from  the  momentum  equations  for  irrotational  (i.e., 
high-frequency)  internal  waves: 

«f  =  -pyp  +  F’  (1.1) 

and 

V{  =  -py/p  +  G,  (1.2) 

where  u  and  v  are  the  horizontal  velocity  components,  p  and  p  are  the  pressure  and  den¬ 
sity  functions,  F  and  G  are  the  nonlinear  and  forcing  terms,  and  the  subscripts  denote 
differentiations.  Cross  differentiation  of  the  momentum  equations  yields 

{Vx-Uy)t  =  Gx-Fy.  (1.3) 

The  quantity  Vj^  -  Uy  is  the  vertical  component  of  the  curl  of  velocity,  i.e.,  the  vertical 
component  of  vorticity.  In  the  absence  of  nonlinearities  and  forcing  (i.e.,  F  =  G  =  0),  the 
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vertical  vorticity  should  be  small  or  zero.  High-frequency  internal  waves  have  little  vert¬ 
ical  vorticity,  and  lower-frequency  ones  have  vertical  components  that  are  less  than  the 
planetary  value  (denoted  as  for  the  Coriolis  parameter,  =  10“^  s“'  at  mid-latitudes). 

Thus,  mapping  the  vorticity  distribution  can  yield  information  about  what  dynamic 
processes  a  water  parcel  has  experienced.  In  studying  ocean  mixing,  for  example,  one  of 
the  difficulties  is  distinguishing  normal,  ambient  structures  that  are  not  undergoing 
significant  stirring  from  structures  that  appear  similar  but  are  undergoing  active  stirring 
and  dissipation.  One  difference  is  that  mixing  events  contain  large  contributions  from 
small-scale  vorticity. 

1.2  Need  for  Vorticity  Observations 

Vorticity  plays  an  important  role  in  the  cascades  of  energy  between  the  large  scales 
of  generation  and  the  smaller  scales  over  which  energy  dissipation  occurs.  Vortical 
motion  is  possibly  the  explanation  behind  much  of  the  variability  observed  in  the  ocean 
that  cannot  be  assigned  to  internal  waves.  On  mesoscales,  say  tens  of  kilometers,  vorti¬ 
city  is  the  most  important  dynamic  variable  in  physical  oceanography.  Features  such  as 
ocean  eddies  often  exhibit  vorticity  comparable  to  the  planetary  value,  but  at  scales  of  a 
meter  and  less  the  principal  signals  are  due  to  mixing  and  nonlinear  processes. 

Historically,  there  has  been  a  need  to  measure  small-scale  vorticity  for  general  stu¬ 
dies  of  turbulence  as  well  as  specifically  for  oceanography.  Saffman  (1981)  claims  that 
"the  invention  of  an  accurate  non-intrusive,  fine-resolution  vorticity  meter  would  be  of 
inestimable  value...,"  while  Muller  et  al.  (1986)  state  that  in  order  to  test  recent  theories 
of  oceanic  finestructure,  "we  need  to  develop  techniques  to  measure  potential  vorticity  on 
small  scales." 

Too  little  is  known  about  the  vorticity  field.  This  situation  results  from  technical 
difficulties  in  observing  vorticity,  not  from  a  conscious  intent  to  ignore  it.  Muller  (1984) 
says,  "The  vortical  mode  also  exists  at  scales  traditionally  associated  with  internal  waves. 
This  basic  fact  has  been  somewhat  repressed  in  the  oceanographic  community."  He 
argues  that  small-scale  vortical  structures  exist  and  may  be  responsible  for  the 
nonintemal-wave  "finestructure"  invoked  to  explain  moored  current  and  temperature 
measurements.  (The  part  of  the  observed  current  and  temperature  variance  that  does  not 
conform  to  behavior  required  of  linear  internal  waves  has  been  denoted  as  finestructure.) 
Muller  states  that  we  need  to  learn  how  to  measure  potential  vorticity,  which  is  a  mea¬ 
sure  of  the  ratio  of  absolute  vorticity  to  the  effective  depth  of  the  vortex.  (Absolute  vorti¬ 
city  is  the  sum  of  planetary  and  relative  vorticity.) 

Another  justification  for  conducting  this  development  is  that  measurements  of  velo¬ 
city  and  vorticity  are  vitally  important  to  understanding  the  role  of  motions  at  scales  of 
1-10  m.  Gargett  et  al.  (1981)  report  a  change  in  the  slope  of  vertical  shear  spectra  at 
about  10-m  vertical  wavelength.  The  horizontal  wavelengths  associated  with  this 
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spectral  "knee"  are  not  known.  The  region  spans  a  gap  between  internal  waves  with  low 
vertical  wavenumbers  and  motions  with  high  vertical  wavenumbers  governed  by  tur¬ 
bulence  and  viscous  dissipation.  New  ocean  observations,  perhaps  of  vorticity,  are 
needed  to  understand  more  about  the  structure  and  physics  of  this  region.  In  addition, 
Muller  makes  the  point  that  the  vortical  mode  is  needed  to  generate  a  complete  descrip¬ 
tion  of  the  internal  wave  regime. 

Lin  and  Pao  (1979)  have  demonstrated  in  laboratory  experiments  that  stratification 
inhibits  vertical  motion.  Thus  a  stratified  fluid  quickly  exhibits  quasi  two-dimensional 
motion  in  the  form  of  horizontal  vortices  which  have  vertical  vorticity.  Their  collapse 
occurs  on  the  order  of  six  buoyancy  periods.  One  hypothesis  is  that,  regardless  of  the  ini¬ 
tial  spectrum,  there  will  be  an  inverse  cascade  of  vorticity  variance.  Thus,  if  there  is  no 
additional  generation  of  vorticity,  the  overall  variance  levels  will  decrease,  but  the 
remaining  energy  will  move  to  lower  wavenumbers. 

Finally,  because  the  natural  background  of  small-scale  vorticity  is  expected  to  be 
small,  observations  of  enhanced  levels  may  indicate  an  unnatural,  artificial  source.  The 
torque  needed  to  spin  up  a  vortex  about  the  vertical  axis  is  not  easily  provided  within  the 
upper  ocean.  Wind  stress  or  water  moving  past  a  surface  are  potential  sources,  and  vorti¬ 
city  is  a  definitely  a  signature  of  wave  instability  and  breakdown. 

In  spite  of  the  advantages  of  directly  measuring  vorticity,  no  sensor  with  that  capa¬ 
bility  at  appropriate  length  scales  and  useful  accuracy  had  been  developed.  This 
deficiency  was  addressed  through  the  development  and  use  of  the  Electromagnetic  Vorti¬ 
city  Meter. 

In  order  for  the  EMVM  to  be  scientifically  useful,  it  must  have  a  noise  level 
significantly  below  the  expected  ambient  level  of  vorticity  in  the  open  ocean.  It  is  possi¬ 
ble  to  estimate  the  energy  levels  of  the  current  finestructure  and  infer  the  needed  sensor 
performance.  From  observations  during  the  Internal  Wave  Experiment  (IWEX),  where 
the  presumption  of  finestructure  was  first  invoked  (Briscoe,  1975),  the  nonlinear  contam¬ 
ination  is  estimated  to  be  about  1  cm/s  rms  at  high  frequency  (-1  cph).  A  simple  model 
for  the  relative  vorticity  of  a  circular  vortex  of  peak  speed  1.4  cm/s  and  radius  L  to  that 
speed  is 

C= — ~ — xlO'^s  .  (1.4) 

The  factor  of  two  comes  from  the  contributions  of  both  velocity  components.  In  IWEX, 
the  coherence  at  high  frequency  dropped  rapidly  at  less  than  10-m  separations,  so  L  must 
be  around  10  m  or  less.  Then  ^  would  be  about  0.003  s”^ ,  a  level  slightly  higher  than 
demonstrated  sensor  performance. 
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2.  CONCEPT  AND  OPERATING  PRINCIPLES 
2.1  Concept  of  EMVM 

The  EMVM  operates  on  the  principles  of  motional  induction  which  govern  the  elec¬ 
tric  fields  induced  as  water  moves  through  a  magnetic  field.  These  principles  are  the 
basis  of  the  conventional  electromagnetic  current  meter  which  generates  a  magnetic  field 
and  measures  the  potential  differences  between  spaced  electrodes.  In  this  case,  only  the 
first  derivatives  of  the  potential  field  are  observed.  In  the  EMVM,  we  measure  the  first 
(V(j))  and  second  (V^(|))  derivatives  of  the  electric  potential  field.  The  first  differences 
yield  the  usual  linear  velocity  components;  the  second  provide  the  determination  of  a 
component  of  water  vorticity. 

The  first  mention  of  measuring  V^(t)  to  estimate  a  component  of  small-scale  vorticity 
appears  to  be  by  Grossman  et  al.  (1957).  Baker  (1971)  described  the  construction  and 
use  of  a  "vortex  probe"  based  on  the  two-dimensional  Laplacian  of  the  potential  field. 
Tsinober  et  al.  (1987)  described  a  very  small  (=1  mm)  probe  used  to  study  turbulence  in 
salt  water.  They  reported  good  success  with  the  method  and  offered  the  following  list  of 
advantages: 

•  It  offers  the  possibility  of  distinguishing  between  velocity  components  because 
of  the  vector  nature  of  the  basic  relation  and  the  sensitivity  to  the  direction  of 
the  flow. 

•  It  is  linear  and  does  not  require  any  calibration  procedure  (in  fact,  no  calibration 
procedure  can  be  applied  to  this  method). 

•  It  is  insensitive  to  the  physical  properties  of  the  fluid  medium. 

•  It  yields  an  instantaneous  measure  of  quantities  fluctuating  in  time,  i.e.,  its 
response  time  is  extremely  small. 

•  It  operates  in  flows  of  a  complex  nature,  and  it  allows  the  determination  of  the 
separation  and  reattachment  points  and  the  moment  and  position  of  the  tur¬ 
bulence  onset  or  relaminization. 

•  It  can  be  used  in  nonisotropic  turbulence  for  determining  the  spectral  interval 
where  the  turbulence  becomes  locally  isotropic  by  using  high-pass  filters  and 
checking  if  the  relations  between  turbulence  components  are  valid. 

2.2.  Governing  Equations 

The  governing  equation  for  the  motionally  induced  potential  differences  in  seawater 
is 

V<l)  =  vxB-J/a,  (2.1) 

where  (j)  is  the  electrostatic  potential,  v  is  the  water  velocity,  B  is  the  magnetic  field,  J  is 
the  electric  current  density,  and  a  is  the  electrical  conductivity  (assumed  constant).  The 
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EMVM  works  by  measuring  the  electric  field  induced  by  the  vorticity  of  seawater  in  a 
magnetic  field.  The  magnetic  field  induced  by  J  is  assumed  to  be  negligible  and  the 
external  magnetic  field  to  be  uniform  and  stationary.  Since  V-J  =  0,  then,  by  taking  the 
divergence  of  Eq.  (2. 1),  we  obtain 

V^(t)  =  V  •  (vxB)  =  B  •  (D- V  •  VxB .  (2.2) 

where  ©  is  the  vorticity  of  the  water  (i.e.,  Vxv).  The  last  term  in  Eq.  2.2  can  be  shown 
to  be  small  compared  with  the  next-to-last  term,  and  can  be  ignored.  The  magnetic  field 
produced  by  the  EMVM  is  unidirectional  or  symmetrical  about  an  axis,  so  that  only  one 
component  of  the  vorticity  is  obtained.  Thus 

V2(t)=:B3©3,  (2.3) 

where  the  subscript  3  denotes  the  component  parallel  to  the  magnetic  field,  which  in  the 
instrument’s  usual  configuration  will  be  vertical.  If  we  can  measure  V^(|)  and  know  B, 
then  we  can  determine  the  component  of  vorticity  parallel  with  B. 

2.3  Measurement  Approach 

We  have  no  sensor  for  measuring  V^(j)  directly,  but  we  can  estimate  it  by  a  finite- 
difference  scheme.  An  array  of  electrodes  is  used  to  sense  voltages  that  result  from  vorti¬ 
city  in  the  presence  of  an  applied  magnetic  field.  In  its  simplest  configuration,  the  instru¬ 
ment  measures  the  x  component  of  V^<t)  (i.e.,  as  shown  in  Figure  2.1.  The  vol¬ 

tages  across  three  collinear  electrodes  are  measured.  Each  outer  electrode  is  referenced 
to  the  center  (ground)  electrode,  resulting  in  the  signed  voltages  shown.  The  output  of 


Vx  =  [e(.^)-e(->^)]-GAIN 

Figure  2.1.  Measurement  scheme  for  determining  A^(t)/Ax^,  the  x  component  of  V^(j). 
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the  electronics  circuit  is  the  difference  in  the  voltages  between  the  electrode  arms.  Thus 
the  device  is  insensitive  to  a  uniform  electric  field  such  as  that  caused  by  a  uniform  flow 
of  water  or  the  lateral  translation  of  the  sensor  through  the  surrounding  seawater.  Rather, 
a  voltage  will  arise  only  if  the  seawater  has  a  relative  rotary  motion.  It  makes  no 
difference  whether  the  sensor  or  the  water  rotates. 

The  measured  voltages  can  be  expressed  in  terms  of  the  potential  (j).  This  potential 
difference  is  given  as  (refer  to  Figure  2.1) 


0 

1 

s 

II 

(2.4) 

and 

1 

T 

II 

-e- 

T 

1 

0 

(2.5) 

The  finite  difference  approximation  of  d^<^/dx^  is  given  by 

A^d)  m  +  ¥-l)-2W) 

Ax^  f 

(2.6) 

This  can  be  expressed  in  terms  of  the  measured  voltages  as 

A^d)  _  eil)-e(-l) 


Thus  the  approximation  for  d^(|)/9x^  in  terms  of  the  output  Vx  is 


dH  _  Vx 

dx^  /^xGain 


(2.7) 

(2.8) 


This  analysis  can  be  extended  to  the  generalized  three-dimensional  case  in  which  three 
independent  measurements  are  made  and  combined  to  yield  a  three-dimensional  approxi¬ 
mation  of  Figure  2.2  illustrates  this  process. 


2.3.  J  2-D  Case:  Flow  and  magnetic  field  uniform  in  z  direction 

Consider  the  case  of  a  body  of  seawater  in  solid-body  rotation  about  its  vertical  axis 
in  a  uniform  vertical  magnetic  field.  A  cylindrical  coordinate  system  is  used  (r,0,z). 
Denote  the  rotation  rate  as  '/2©2  (the  V2  is  there  because  the  rotation  rate  is  twice  the  vorti- 
city  for  solid-body  rotation).  This  situation  is  equivalent  to  rotating  an  electromagnetic 
sensor  in  stationary  seawater  in  the  presence  of  the  vertical  magnetic  field,  a  method  used 
to  evaluate  our  EMVM.  For  rotation  about  the  r  =  0  axis,  the  horizontal  velocity  in 
solid-body  rotation  is  'Ato^r  le,  where  6  is  the  angle  in  the  cylindrical  coordinate  system 
(positive  CCW  from  x  axis).  The  electric  potential  function  fromEq.  2.1  is 

r 

^{r,Q)  =\{V2Bz(Si^r-Jx/<5)dr.  (2.9) 

0 


The  absence  of  vertical  variations  in  the  magnetic  field  or  flow  structure  means  that 
and  dJ^/dz  are  zero.  Continuity  of  electric  current  requires  then  that  dJ^/dr  is  also 
zero.  Axial  symmetry  requires  that  7^  =  0  at  r  =  0.  Hence,  J  =  0  everywhere.  Then, 
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z 


V^(t>(x,y,z)  =  (t)(Jl,0,0)  +  (!)(-«,, 0,0)  +  (f)(0,Jl,0)  +  ())(0,-^,,0)  +  (t)(0,0,J!-)  +  (f.{0,0,-i!,) 
-  64>(0,0,0)  / 

=  |ie(i!.,0,0)  -  e(-5,0,0)  +  e(0,Jl,0)  -  e(0,-£,0)  +  e(0,0,i!.) 

-  e(0, 0,-11))]  / 


Figure  2.2.  Three-dimensional  finite-difference  approximation  of  V^(|). 


(|)(r,e)  = 


(2.10) 


A  set  of  three  electrodes,  equally  spaced  on  a  line  in  the  horizontal  plane  and 
located  at  r-l,r,  and  r  -i-  Z  would  sense  an  average  difference  voltage  of 


^  lB,(i),[(r  +  lf  +  (r-lf-2r^]  (2.11) 

2  o 


=  . 


It  can  be  shown  that  any  three  equally  spaced  electrodes  on  a  line  in  the  horizontal 
plane  would  yield  the  same  voltage  regardless  of  orientation  or  reference  to  the  origin  of 
the  arbitrary  coordinate  system.  An  easy  way  to  confirm  this  is  to  consider  a  triplet  of 
orthogonal  electrodes  positioned  so  that  the  center  electrode  is  at  r  =  0.  Each  axis  contri¬ 
butes  the  same  difference  voltages,  and  this  condition  does  not  change  as  the  electrode 
array  is  moved  away  from  the  origin. 
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A  second  horizontal  axis  would,  therefore,  observe  the  same  voltage.  Thus,  a  two- 
axis  unit  would  measure  an  average  voltage  of  A  set  of  three  electrodes 

equally  spaced  along  the  vertical  axis  will  sense  zero  potential  difference  and  contribute 
nothing  to  the  sum  obtained  with  two  sets  of  electrodes  oriented  horizontally,  say  along 
the  X  and  y  axes.  Hence,  a  sensor  with  two  axes  in  the  horizontal  plane  would  observe  an 
average  voltage  of  An  average  voltage  is  used  so  that  the  observed  potential 

does  not  increase  linearly  with  the  number  of  sensor  axes. 

2.3.2  3-D  Case:  Uniform  rotary  flow  in  dipole  magnetic  field 

Consider  the  case  of  an  infinite  volume  of  seawater  in  solid-body  rotation  about  the 
vertical  (z)  axis  in  a  vertical  dipole  magnetic  field.  This  is  a  realistic  example  which 
simulates  what  an  EMVM  (with  its  own  dipole  magnetic  field)  would  observe  in  a  patch 
of  vorticity  with  a  spatial  scale  that  is  large  compared  with  that  of  the  sensor.  Also,  it 
corresponds  to  the  rotation  of  an  EMVM  that  produces  its  own  dipole  magnetic  field  in  a 
comparatively  large  tank,  such  as  was  done  in  the  laboratory  to  check  the  instrument’s 
performance.  More  complex  magnetic  fields  can  be  constructed  by  superposition  from 
this  elemental  solution. 

The  symmetry  of  the  problem  places  constraints  on  electric  current  paths.  For 
example,  the  electric  currents  must  be  antisymmetric  about  the  z  =  0  plane  (i.e.,  mirror 
images).  Consider  the  line  in  the  z  direction  that  passes  through  the  origin.  For  a  volume 
in  solid-body  rotation  about  this  line,  the  vertical  electric  currents  must  change  sign  at 
the  origin.  That  is,  the  currents  are  an  odd  function  with  respect  to  z.  This  fact  can  be 
seen  by  considering  the  induction  terms  along  the  z  axis.  The  induction  term,  v  x  B  •  1^,  is 
zero  on  the  z  axis.  Hence, 

0 

<t)(0)-(l)(-oo)  =  -  \j^dz  (2.12) 

—  oo 

and 

oo 

<t)(0)-(l)(oo)  =  (2.13) 

0 

Because  the  magnetic  field  is  of  finite  size,  it  is  necessary  that  the  potential  be  a  constant 
at  infinity;  that  is,  (])(<»)  must  be  the  same  value  as  <j)(-  <»).  It  is  convenient  to  take  this 
value  to  be  zero.  The  only  form  for  that  will  satisfy  the  above  expression  is  for  it  to  be 
an  odd  function  of  z.  In  this  way,  the  potential  at  the  origin  is  the  same  regardless  of 
whether  the  integration  toward  the  origin  is  along  the  positive  or  the  negative  z  axis. 
Another  consequence  is  that  the  electric  currents  converging  on  the  origin  along  the 
z  axis  are  compensated  by  electric  currents  flowing  outward  on  the  x-y  plane.  Thus  the 
horizontal  current  along  the  x  axis  must  also  be  an  odd  function  of  x  and  of  the  opposite 
sign  as  that  for  z.  That  is,  if  is  positive  for  z  >  0,  then  7,  must  be  negative  for  x>  0. 
The  induction  term,  vxB  •  1;^,  integrates  to  zero  along  the  line  from  infinity  to  the  origin. 
Hence  it  makes  no  contribution  to  the  potential  at  the  origin. 
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The  consequence  is  that  integrating  from  infinity  along  the  z  axis  yields  a  potential 
at  the  origin  that  is  of  opposite  sign  to  that  obtained  by  integration  along  the  x  axis.  The 
potential  at  the  origin  must  be  independent  of  the  path  of  integration.  The  only  way  for 
the  potential  be  to  independent  of  the  path  of  integration  is  for  to  be  zero  on  the  x  axis, 
for  Jy  to  be  zero  on  the  y  axis,  and  for  to  be  zero  on  the  z  axis.  This  argument  does  not 
exclude  electric  currents  flowing  in  other  parts  of  the  space,  but  a  formal  solution  given 
later  demonstrates  that  J  =  0  eveiywhere.  Hence,  there  can  be  no  electric  currents 
flowing  as  a  consequence  of  solid-body  rotation  in  a  dipole  magnetic  field.  The  result  is 
that  for  vorticity  the  electric  current  field  is  zero. 

A  vorticity  of  0.1  s"*  corresponds  to  a  voltage  of  approximately  25  nV  over  5-cm 
arms  in  a  10"^  T  (1  G)  magnetic  field.  The  greatest  difficulties  in  detecting  such  a  small 
signal  are  due  to  the  drift  of  the  electrode  offset  voltages,  electrode  noise,  and  the  1// 
behavior  of  electronic  noise.  Figure  2.3  shows  the  spectrum  of  the  noise  associated  with 
the  electrodes  and  electronics.  The  l//behavior  of  the  electronic  noise  makes  the  detec¬ 
tion  of  low-frequency  vorticity  difficult. 


Figure  2.3.  Typical  noise  level  associated  with  input  amplifier  using  Analog  Devices  OP- 27 
integrated  circuit  operational  amplifiers.  This  configuration  is  similar  to  that  used  for  the  first, 
second,  and  third  versions  of  the  EMVM. 
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A  more  formal  derivation  of  the  finite-difference  equations  can  be  based  on  Gauss’s 
theorem.  Assume  that  the  vorticity  and  magnetic  field  are  axially  symmetric  about  the 
vertical  axis,  for  example,  solid-body  rotation  in  a  vertical  dipole  magnetic  field. 
According  to  Eq.  2.2,  (j)  =  B  •  o .  Integrating  this  equation  over  the  volume  V  and 

applying  Gauss’s  theorem  yields 

jv^-ndA=iB-GidV ,  0  14) 

A  V  t  ; 


where  n  is  the  unit  vector  normal  to  the  surface  A.  Consider  a  body  that  is  a  right  circular 
disk  of  radius  R  in  the  horizontal  or  x-y  plane  and  thickness  dz  centered  about  z  =  0.  The 
left-hand  integral  can  be  expressed  in  terms  of  contributions  from  the  surface  at  r  =  /?  and 
from  the  top  and  bottom  areas: 

JV(|)-nJA  =  j  9(|)/9r<iA -F  J  9<|)/9zdA-F  J  d^/dzdA.  o  1^^ 

AAA  A 

^  ^  lop  ^  bot 

The  contributions  from  the  top  and  bottom  surfaces  cancel  because  of  the  assumed  sym¬ 
metry  of  the  flow  and  of  the  magnetic  field  about  the  vertical  axis.  The  potential  and 
radial  derivative  at  R  can  be  expressed  in  terms  of  MacLaurin  series  which  can  be  mani¬ 
pulated  to  yield 

2 

3(|)(r)/9r  =  ~[  (t)(T)  - (t)(0)  ]  -  d^/dr  I  o  +  I  o  +  higher  order  terms .  (2. 16) 

r  6 

From  the  earlier  analysis,  we  know  that  (t)(r)  =  ®  r^.  The  term  0(|)/9r  I  q  is  zero,  and 

the  third-  and  higher-order  derivatives  are  zero  (and  will  be  small  in  general  for  vorticity 
patterns  that  are  larger  than  the  sensor).  Thus  we  obtain 

J[(|)(R)-(t)(0)]cM  =  f-jB-coc/\/  •••  . 

^  V 

Dropping  the  higher-order  terms  and  expressing  the  integrals  as  averages  gives 

<  (\>iR)-m  >A,  =  <  B  •  0  >v , 

Ar 


(2.17) 

(2.18) 


where  <  >^4^  represents  an  average  over  the  surface  Ar  and  <>y  represents  an  average 
over  the  volume  V. 

The  ratio  of  the  volume  (i.e.,  V  =  nR^dz)  to  the  area  (i.e.,  Ar  =  InRdz)  of  the  disk¬ 
shaped  body  is  'A/?;  thus  Eq.  2.18  becomes 

d2 

<  (t)(R)-<l)(0)  ^  <  B  •  0  >v. .  (2.19) 

For  a  thin  disk  {dz  <<R),  the  "volume"  average  is  more  nearly  an  average  over  the 
surface  between  the  origin  and  the  radius  R,  and  the  "surface"  area  is  more  nearly  an 
average  around  the  circumference  at  the  radius  R.  If  electrodes  are  added  off  the  x-y 
plane,  the  volume  changes  as  does  the  bounding  surface. 
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The  simplest  approximation  of  a  spherical  array  is  six  electrodes  at  distances  of±R 
on  the  X,  y,  and  z  axes  and  one  reference  electrode  at  the  origin.  The  corresponding  area 
and  volume  are  for  the  sphere  of  radius  R  about  the  origin.  Recall  that  the  vertical  poten¬ 
tial  gradient  is  zero,  and  electrodes  separated  in  that  direction  make  no  contribution  to 
the  average.  Therefore,  the  average  potential  difference  is  2/3  that  of  a  single  axis.  The 
volume  average  of  B  •  eo  times  the  ratio  of  volume  to  area  is  1/3  B  ■  (O.  Hence,  for  the 
2-D  case  where  the  volume  integral  is  trivial. 


p2 

<  ^(R)  -  (|)(0)  >^^  =  -^  <  B  •  (0  >  V,  (2.20) 

in  confirmation  of  Eq.  2.19. 

The  utility  of  the  vertical  electrodes  appears  when  flow  is  not  two-dimensional  and 
perpendicular  to  the  axis  of  B. 

We  interpret  our  measurements  in  terms  of  the  integral  expression  because  it  expli¬ 
citly  accommodates  spatially  variable  B  and  ©.  In  summary,  the  average  potential  from 
a  set  of  equally  distributed  electrodes  (minus  the  potential  at  the  origin  of  the  zirray)  is 
interpreted  as  equal  to  the  average  scalar  product  of  magnetic  field  and  vorticity.  We  use 
the  robust  relationship  given  in  Eq.  2.19  for  converting  average  electric  potential 
differences  into  a  component  of  the  average  vorticity  parallel  to  the  known  applied  mag¬ 
netic  field. 


2.4  Modeling 

Section  2.3  derived  expressions  between  potential  measurements  and  vorticity  with 
a  minimum  of  mathematical  complexity.  In  this  section,  analytical  solutions  are  derived 
for  uniform  flow  and  uniform  solid-body  rotation.  In  this  analysis,  we  will  assume  that 
the  magnetic  field  is  generated  by  a  thin  coil.  Using  the  assumptions  made  earlier,  the 
governing  equations  are 

V<t)  =  vxB-J/a  and  V^<^  =  B-'Vxv.  (2.21) 


The  goal  is  to  develop  analytical  expressions  for  ([)  and  J  based  on  specified  forms  of  B 
and  V. 


The  magnetic  field  vector  B  produced  by  an  electric  current  (/)  in  a  thin  coil  of 
radius  a  can  be  expressed  as  two  infinite  series  of  Legendre  polynomials  in  a  spherical 
coordinate  system  (r,  0,  y).  The  first  is  for  the  domain  within  a  sphere  defined  by  the  coil 
(denoted  hereafter  as  the  "inner"  region). 


1  to  < 


«+l 


Br  =  (-1)  '  (n  +  l)n 

nodd 


n 

JL 

n  +  1 

a 

n-\ 


P„  (co50)  (2.22a) 


Ba  = 


2a  AA  n 

n  odd 


n 

n  +  1 

a 

n-\ 


Pl  (COS0)  (2.22b) 
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and  for  the  domain  beyond  the  coil  sphere  (the  "outer"  region). 


^  n-r  1. 

Ilf  1  to  <»  I 

=  I  (-1) '  («  +  i)n 

nodd 


n  +  1 


P„  (cos6)  (2.22c) 


n/  1  to~  — 

=  -^  E  (-1)  '  n 

nodd 


n  +  1 


pI  (cosG)  , 


(2.22d) 


where  p  is  the  magnetic  permeability  and  P„  and  Pj  are  Legendre  and  associated  Legen¬ 
dre  polynomials,  respectively.  The  product  term  here  will  denote  for  n=l,  for 

n  =  3,  etc.  Only  the  values  of  n  used  in  the  summation  (here  n  being  odd  only)  will  be 
used  in  the  product  term.  The  magnetic  field  computed  by  Eq.  2.22  is  displayed  in  Fig¬ 
ure  2.4. 

The  boundary  conditions  are  (1)  that  <])  is  bounded  at  the  origin,  (2)  that  (j)  is  zero  at 
infinity,  and  (3)  that  the  radial  components  of  the  V<t)  and  J  vectors  are  continuous  across 
the  boundary  between  the  inner  and  outer  regions.  We  will  now  examine  the  solution  for 
two  simple  flow  regimes. 

2.4.1  Uniform  flow 

Assume  a  flow  of  uniform  horizontal  velocity  (U),  with  no  vertical  component  or 
variation.  Then 


y  =  U  f  =  U  sin6  cosy  r  +  U  cos0  cosy  0  -  [/  siny  y . 


(2.23) 


For  uniform  flow,  the  velocity  field  is  irrotational  (V  x  v  =0),  so  the  problem  reduces  to 
solving  for  the  Laplacian  of  the  potential,  V^<j)  =  0,  with  the  conditions 

(r  =  a ) (|),  =  (l)o  and  (r  =  a)  vxB-V(t)  ,r=  vxB-V(t)  ^r,  (2.24) 

where  the  subscripts  i  and  o  denote  inner  and  outer  regions,  respectively.  First,  we  will 
examine  the  radial  component  of  v  x  B  in  order  to  deal  with  the  second  boundary  condi¬ 
tion.  In  spherical  coordinates. 


y  xB  =  -VyBQ?  +  VyBfQ  +  (VrBQ  -VQBr)y. 


For  the  inner  region, 

(vxB)P  =  ^‘|”  (-1) 

for  the  outer  region, 


n  +  I 


(2.25) 


sinyP^  ;  (2.26a) 


(vxB)P  =  --^'x“(-l)  2 

n  odd 


sinyP^ .  (2.26b) 


The  second  boundary  condition  (at  r=  a)  can  then  be  written  as 
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Figure  2.4.  Magnetic  field  produced  by  an  electric  current  in  a  thin  coil.  Distances  are  normal¬ 
ized  by  a,  the  radius  of  the  coil.  The  current  is  running  counterclockwise  through  the  coil,  as 
viewed  from  above.  The  field  has  been  transformed  into  a  x-z  coordinate  system.  The  contour 
values  and  scales  are  arbitrary,  and  gray  shading  denotes  negative  values,  (a)  Contours  of  the 
vertical  component  of  B.  (b)  Contours  of  the  horizontal  component  of  B.  (c)  Vectors  of  B. 
Note  that  the  vectors  are  nearly  uniform  in  both  direction  and  magnitude  near  the  center  of  the 
coil. 
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^12^  f^l 

dr  dr  2a  n  n  +  1  i  n  \  j 


The  general  solution  for  the  inner  region  is 


<!>/  =  Z  E 

n  =0/n  =  0 


W  ye  ,  ci  yO 
mn  ^  mn  ^ mn  ^  mn  5 


where  C  and  S  are  as  yet  undetermined  coefficients  and 

Ymn  =  Pn  ( COS0 )  cosm y  and  Y%„  =  Pn(  cos0 )  sinm y 


(2.28) 


(2.29) 


are  spherical  harmonics.  In  a  similar  manner,  the  general  solution  for  the  outer  region  is 


=  E  E 


/~'0  ye  ,  CO  yO 


n=0m=0 


So,  the  first  boundary  condition  can  be  expressed  as 


n  =  0  m  =  0 


(2.30) 


2  2  =  2  X  ,  (2.31) 

n=0/n=0  /z=0m=0  ^ 

while  the  second  boundary  condition  gives 

XV’  —  /~'i  ye  ,  ci  yO  _  _  ^  (w  +  1  )  ye  ,  co  yO  1 

n=0m  =  0  ^  a  \.  J 


+  -^'e'”(-i)  '  n  -77 

2a  n  n  +  1 


(2.32) 


After  dividing  through  by  a,  the  last  term  can  be  expressed  as 

iM'lf-,  yO 

^  2j  ^  In  ^  In  ^ 

^  n  odd 

where,  for  n  odd. 


(2.33) 


Ai„  =  (-1) 


n  and  =  sinyP^  .  (2.34) 

n  +  1 


Since  the  only  external  constraint  on  our  problem  comes  from  the  v  x  B  portion  of  the 
boundary  condition,  and  it  contains  only  terms,  we  can  reduce  our  boundary  condi¬ 


tions  to 


CO  oo 

Xci  vO 

In  -  2/  ■Jln-'^ln 


(2.35a) 
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and 


i  «  =  -  i  (n  +  1 )  Sf, z"/l  1,  .  (2.3Jb) 

n=0  n  =  0  n  odd 

Equating  similar  degrees  of  Yi„  and  setting  all  even  degrees  of  5i„  to  zero,  we  obtain 


and  n  = -(n  +  1  )5?„  + -^  A,„ 

We  now  have  two  equations  and  two  unknowns,  which  yield 

_  UiU_ 


S\n  =  = 


2  (  2n  +  1 ) 


The  expressions  for  the  potential  0  in  the  inner  and  outer  regions  are  then 


and 


n  odd 


UlLl 

=  E  (-1) 

^  n  odd 


n 

r_ 

n  +  l 

a 

sinyP„ 


n  +  l 

2  J_ 

n 


n 


n 

n  +  l 

r 

n  +  l 


(2.36) 


(2.37) 


(2.38a) 


sinyPji .  (2.38b) 


The  0  field  is  shown  in  Figure  2.5.  For  uniform  flow,  the  potential  field  has  no  curvature 
in  the  center,  and  thus  no  vorticity  would  be  measured.  This  is  expected  for  uniform 
flow.  However,  to  infer  velocity  from  V  0,  the  current  density  must  be  known.  Equation 
2.21  is  rewritten  as 


V0  =  V  X  B 


J/g 

vxB 


(2.39) 


where  Hd  is  a  correction,  or  "head"  factor.  The  above  division  process  expresses  the  fol¬ 
lowing  form:  Hd  =  J/o  •  /  /(v  x  B)  •  /,  where  i  is  a  unit  vector  in  the  ith  direction.  This 
factor  will  be  computed  next. 


The  three  components  of  the  electric  current  density  are 

30 


Jr  =  g 


je  =  g 


-Me- 3^ 


VyBf  - 


1  ^ 

r  30 


and 


(2.40a) 

(2.40b) 
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Figure  2.5.  The  potential  field  <j)  generated  by  uniform  flow  in  the  presence  of  the  magnetic 
fields  shown  in  Figure  2.4.  The  field  has  been  transformed  into  an  x-z  coordinate  system.  Again, 
gray  shading  denotes  negative  values.  The  view  is  with  the  flow  coming  out  of  the  page.  The 
fine-scale  perturbations  are  an  artifact  of  the  Legendre  polynomial  solution.  The  contour  values 
are  arbitrary. 


(2.40c) 

(2.41a) 

(2.41b) 

(2.41c) 
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^Tc-i)  ^  n  ^ 

„  odd  "  +  1 


(n  +  1)P„  + 


TJ|7J  “"T" 

(n  +  l)Pj[_i  —ncosQPlt 
nsinG 


«  +  l  n—\ 

aUlil  n  r 

JiY=cosY-:^  X  (-1)  n  ~T7  “ 


(n  +  l)coseP„+  , 

«sin9 


,  -  cUiil  ^  .  X  9  „ 

JoY=cosY— ^  X  (-1)  n 

n  odd 


n 

a_ 

n  -F  1 

r 

/  1  \  rv  «  sin^6  +  1  -.1 

(n  +  l)coseP„- - — — Pi 

nsinG 


(2.4  Id) 


(2.41e) 


(2.41f) 


The  J  fields  are  plotted  in  Figure  2.6,  while  Hd,  the  ratio  of  J/g  and  v  x  B,  is  plotted  in 
Figure  2.7.  Near  the  center  of  the  coil,  the  values  of  Hd  are  nearly  uniform  at  0.5.  Thus 
the  velocity  estimated  from  V  (j)  measurements  alone  is  only  half  what  exists. 


2.4.2  Solid-body  rotation 

Let  us  now  assume  a  flow  of  solid-body  rotation,  centered  about  the  z  axis.  For  a 
spherical  coordinate  system,  the  vorticity  co  of  the  system  is 

coz  =  ©cosGr  -  ©sinGG .  (2.42) 


The  problem  now  is  to  solve  the  Poisson  equation,  V^<t)  =  B  •  co,  with  the  constraints  that 
the  potential  (|)  is  bounded  at  the  origin  and  is  zero  at  infinity  and  that  the  potential  as  well 
as  the  radial  component  of  J  is  continuous  across  the  interface  between  the  inner  and 
outer  regions.  The  solution  scheme  will  proceed  as  follows.  First,  we  will  rewrite 
V^({)  =  B  •  (D  as 

( (|)p  +  (|)// )  =  B  •  CO ,  (2.43) 


where  =  B  •  ca  and  =  0,  and  the  subscripts  P  and  H  denote  the  particular  and 
the  homogeneous  solutions.  We  will  first  solve  for  ^p,  beginning  with  the  right-hand  side 
of  our  inhomogeneous  equation,  where,  after  regrouping  terms  and  simplifying,  we 
obtain 
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Figure  2.6.  Vector  plots  of  the  electric  current  density  J  produced  by  a  uniform  flow.  The 
fields  have  been  transformed  into  an  jc-y-z  coordinate  system,  (a)  View  in  the  x-y  plane,  looking 
down  from  above.  The  uniform  flow  is  from  left  to  right  (positive  x  direction).  Within  the  coil,  J 
is  in  the  same  direction  as  vxB.  Since  V(t)  =  vxB  -  J/a,  the  value  of  V (|)  is  less  than 
expected  from  vxB  alone,  (b)  View  in  the  y-z  plane,  looking  from  the  positive  x  direction. 
There  is  no  vertical  component  of  J  in  this  slice. 
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Figure  2.7.  Contours  of  the  head  factor  Hd.  The  light  gray  shading  denotes  negative  values 
(whose  contours  are  not  shown),  while  the  dark  gray  indicates  the  region  where  this  component 
of  the  head  factor  is  between  0.45  and  0.55.  (a)  Contours  in  the  x-y  plane,  (b)  Contours  in  the 
x-z  plane.  Note  that  for  most  of  the  region  within  the  coil,  the  head  factor  is  nearly  0.5. 


TR9503 


19 


•  UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY- 


=  („  +  i)n  ^ 

n  even 


-  Pn- 
a 


(2.44) 


The  product  term  is  defined  here  to  be  equal  to  1  when  n  =  0.  Since  this  equation  forms 
the  right-hand  side  of  a  Poisson  equation  and  contains  terms  of  the  form  P„  and  r/a,  we 
will  try  a  form  of  (|)p  that  has  these  same  characteristics.  Note  that  the  right-hand  side  has 
no  dependence  on  the  angle  y.  So,  let 

frl"*" 

I  ^  I  '  i  l-k  X  .  /\  \  A 


n=0 


P„(COS0)  , 


(2.45) 


where  the  C„  are  undetermined  coefficients  and  m  is  a  constant  power  to  be  determined. 
The  Laplacian  of  (|)  in  spherical  coordinates  is 


V  ^ip  = 


1  a 

sin6 

sin0  39 

59 

^  1  .  _ 1 _  ^  n  Af,\ 

59  r^sin^G  57^ 


If  we  use  the  assumed  form  of  the  solution,  the  Laplacian  of  (t),p  becomes,  after  consider¬ 
able  algebra, 


V^<t>/P  =  T  2  Q  [2(«  m)  +  {n  +  m){n  m  -  \  )-n{n  +  1  )] 


Pn-  (2.47) 


Now,  equating  the  left-hand  side  of  our  Poisson  equation  with  the  right-hand  side,  for 
n  =  0,  we  obtain 

r  f 

—  [2m  +  m(m-l)]  —  .  (2.48) 

(X  2aCI 


In  order  that  the  powers  of  r  cancel  on  each  side  of  the  equation,  m  is  required  to  equal  2. 
Thus, 

^[4  +  21  =  ^  or  C„  =  ^.  (2.49) 


Because  the  right-hand  side  of  our  equation  is  exactly  zero  for  odd  values  of  n,  then 
C„  odd  =  9-  Iio  general,  C„  can  be  expressed  in  the  form 

c„  =  aV(")=  (2-50) 

2a  2 


where 


1  n  -I-  1  ^  n  -  1 

■^<"*=I(2n+3)n  „ 


(2.51) 
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So,  plugging  our  value  for  C„  back  into  our  original  expression  for  <l),p,  we  obtain 

lirOtOoo  ~  mJ.1  in  _  ll 


=  -2±i-n[A^|  [i]"*' 


P„(cose).  (2.52) 


Now,  let  us  repeat  the  previous  work,  but  for  r/a  >  1  outside  the  thin  coil.  The  left-hand 
side  is  now 


co-b=--^^|"(-i)2  («)n  [  —  1  [-1 

2a  n  r 

n  even 


Pn .  (2.53) 


In  a  manner  similar  to  before,  let 


^oP  =  Yd  P^n 

n=0 


P„(cose). 


The  Laplacian  then  becomes 

V^(j)„p  =  ^- 2 \(,n+m){n+m  +  \)-2{n+m)-n{n  +  \)  — 

y.  Jr 


n  +  m  +  2 


(2.54) 


Pn  •  (2.55) 


Matching  terms  and  realizing  that  D„  =  0  for  all  odd  n  and  for  n  =0,  then  for  n  =  2  we 
obtain 

r%  f  ^  2  + ftt-t-2  _  ft3 


[(2  +  /n)(3  +  m)-2(2  +  m)-2(2+  1)]  — 


_a_ 

2  +  m  +  2 

(Op/ 

R. 

r 

2a 

r 

(2.56) 


In  order  that  the  powers  of  r  cancel  on  each  side  of  the  equation,  m  is  required  to  equal 
-1.  Thus 


V[2-2-6]  -  =-^  ^ 

r  2a  r 


(2.57) 


Z)2=- 


_  _  _af_  cou7  _  _  (oau/ 
~  6  2a  ~  12 


In  general,  D„  can  be  expressed  in  the  form 


(-1)^  g{n), 


where 


.X  1  n  „  n-l 

8in)  =  -—— - —  n  - 

2  (2n  -  1 )  n 


(2.58) 


(2.59) 


(2.60) 
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So,  finally, 


_  coau/ 


^  n  even 


n-1  a 


P„(cos9).  (2.61) 


Since  the  inner  and  outer  particular  solutions  of  the  potential  ( ([)//>  and  (j)^^  )  do  not  match 
at  the  interface  r  =  a,  we  will  add  a  solution  of  Laplace’s  equation  such  that  the  com¬ 
bined  potential  (homogeneous  and  particular  solutions)  is  continuous.  As  an  additional 
constraint,  we  will  force  the  induced  current  density  to  be  continuous  across  the  interface 
as  well. 


For  the  inner  region,  the  solutions  are 


^iP  = 


_  coau/  (-1)  2  ”  +  1 


2n  +  3 


Pfi  (cos0)  (2.62a) 


=  I,  Art  —  Pn  (cose)  , 


(2.62b) 


where  will  be  determined  to  satisfy  the  boundary  condition. 
For  the  outer  region,  the  solutions  are 


coau/ 


r  2; 

^  n  even  ^ 


a 


Pn  (cosG)  (2.63a) 


<t>oH  =  X 


P„  (cos0) . 


(2.63b) 


where  5„  will  be  determined  to  satisfy  the  boundary  condition.  Now,  we  must  ensure 
that  the  potentials  match  at  the  interface,  or 

i^iP  +  <t>)i/)r  =  a  =  i^oP  +  ^oH)r  =  a  •  (2.64) 


For  n  even,  we  match  degrees  of  Legendre  polynomials,  so 


n  +  \ 
2n  +  3 


(2.65) 


COa[i/  /  1  \  2  - Ik -  |-r 

4  ^  2n-l  n 
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(2n  +  l) 


A  -B  +  . 

"  4  ^  (2n-l)(2n  +  3) 


n 


n  - 1 


n 


(2.66) 


Next,  we  ensure  that  the  radial  derivatives  match  at  the  interface.  From  the  governing 
equation,  we  have 

J  =  a(vxB-V<t)).  (2.67) 

So,  to  ensure  that  the  radial  component  of  J  is  continuous,  we  have  the  constraint  that 


(v  X  B  )i  -  (((jj-p  + 


r  = 


J  r  =  a 


(v  X  B)^  (^oP  +  ^(?//) 


(2.68) 


r  =  a 


Our  model  flow  is  solid-body  rotation  with  vorticity  O);  thus  it  has  an  angular  rate  of  co/2 
and  a  velocity  of 


v  =  Y  (rsine)Y. 


So 


\XB  =  -VyBQT  +  VyBrQ  . 

The  radial  component  of  v  x  B  for  the  two  regions  can  be  written 

(^XB)  P  =  (-1)^  {4nUSn^  +  n-2) 

^  4  ^  (2n-l)(2n  +  3)  ^ 


(2.69) 


(2.70) 


n  - 1 


n 


Pn  (2.71a) 


n  - 1 


n 


Pn  •  (2.71b) 


,  -  COU/  ,  ..0  (4n^ +4n^-3n-l)  ^ 

(vxBU  =  -f  JJ-l) 

Now,  matching  degrees  of  Legendre  polynomials  across  the  interface, 

„  4  -  ( „  4.  n  B  t  I  'i'^  ( 2n  1 1  K  t  2b  - 1 )  „ 

-bA.  -  (n  +  DB.  -  ^  (-1)  (2„.ix2„43)  H 

Now  we  have  two  equations  and  two  unknowns.  Using  Eq.  2.66  and  Eq.  2.72,  we  obtain 


n  - 1 


.  (2.72) 


n 


(2n-l) 


n 


n  - 1 


and 


im\U  2  T-T 

4  <  (2«  +  3)n 


n 


n-\ 


(2.73a) 


(2.73b) 
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All  the  odd  degrees  of  A„  and  are  exactly  zero.  The  final  forms  of  the  inner  and  outer 
homogeneous  solution  are 


and 


^  (mill 

=  4  S 

n  even 


n  p. 

n-1 

r_ 

(2n-l) 

n 

a 

cos(e)  (2.74a) 


COatt/  2  JgjdL  |-r 

4  (2n.3)n 


n-l 

n 


r 


n  +  \ 


P„cos(e).  (2.74b) 


Finally,  combining  the  homogeneous  and  the  particular  solutions  of  the  potential  results 
in 


([),  = 


coau/  ^  +  1  |-r 

^  ^  ^  ^  2n  +  3 


n  even 


n  - 1 
n 


^  n  +  2 


r_ 

a 


Pfi  -Pn+2 


(2.74c) 


and 


<1>^  =  E  (-1) 

^  n  even 


n 

n-1 

n-\ 

7  n 

A 

Pji  -Pn-2 

2n  - 1  n 

n 

r  \ 

1  J 

.  (2.74d) 


The  potential  field  (j)  for  the  velocity  field  in  solid-body  rotation  is  shown  in  Figure  2.8. 


\ 

(^1 

(fe 

\ 

-3  -2-10  1  2  3 

X  Distance 


Figure  2.8.  The  potential  field  <[)  generated  by  flow  in  solid-body  rotation  in  the  presence  of  the 
magnetic  field  shown  in  Figure  2.4.  The  field  has  been  transformed  into  an  x-z  coordinate  sys¬ 
tem.  Again,  contour  values  are  arbitrary.  The  view  is  with  the  flow  coming  into  and  out  of  the 
page.  Unlike  for  uniform  flow,  (j)  here  is  everywhere  positive. 
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To  compute  the  induced  electric  currents,  we  plug  our  expressions  for  v,  B,  and  (}) 
into  our  governing  equation  vector  expression  of  a  (v  x  B  -  V  (j)),  and  obtain 

n  +  l 


^  n  even 


n  - 1 

r_ 

n 

a 

(2.75) 


sin0  P} 


n  +  l 


(n  +  l)(n  +  2) 
2n  +  3 


Pn-Pn+2 


+  sin0 


00)1x7 

4 


0  to  oo 

E  (-1) 


n  even 


n  +  l  „ 

n  - 1 

r_ 

2n  +  3  n 

n 

a 

(2n  +  3)  P„  + 1  + 


P'  -P' 

^  n  ^  i 


71+2 


0. 


It  can  be  shown  that  both  the  terms  in  brackets  are  identically  zero.  Thus,  the  current 
density  J  is  zero  in  the  inner  region  when  the  flow  is  in  solid-body  rotation. 

Now,  for  the  outer  region,  the  vector  addition  of  a  (v  x  B  -  V  (j))  yields 


J  = 


00)0  |X/ 

4 


2  to  oo 

Z 

n  even 


(-1) "  n 


n-l 

n 

r 

(2.76) 


.  njn-lL 


Pn-Pn-2 


.  „  O(0]X/  2^““  ,  ,^2  _« 


sin0 


n  - 1 

£- 

n 

r 

Again,  both  the  terms  in  brackets  are  identically  zero.  Thus,  the  current  density  J  is 
zero  in  the  outer  region. 


2.5  Response  Function  for  Sensor 

As  shown  earlier  in  Section  2.3.2,  the  EMVM  determines  the  area  average  of  vorti- 
city  in  the  circle  encompassing  a  2-D  electrode  array.  An  isolated  vortex,  say  of  magni¬ 
tude  fiv  radius  ry,  completely  within  the  sensing  area  will  produce  an  apparent  vorti- 
city  of  Qy(ry)/(rj),  where  is  the  radius  of  the  electrode  configuration.  Thus,  the  ampli¬ 
tude  response  function  for  the  sensor  is  unity  for  vortices  with  a  size  much  larger  than 
and  decreases  as  (ry/r^)  for  small  vortices.  The  power  response  of  the  sensor  exhibits  a 
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fallofif  of  the  formyl  (or  k~^,  if  the  frequency  scale  is  converted  to  wavenumber  using 
Taylor’s  hypothesis,  /  =  kU). 


For  vortices  with  a  size  the  same  order  of  magnitude  as  that  of  the  sensor,  the 
response  is  determined  by  the  average  amount  of  overlap  of  the  vortex  and  the  sensing 
area.  The  computation  of  overlap  area  depends  on  the  relative  sizes  of  the  two  circles 
(rv<rj  or  r^>rs)  and  on  the  overlap  geometry  (barely  intersecting  or  broadly  overlap¬ 
ping).  There  are  four  cases  for  overlapping  circles,  as  depicted  in  Figure  2.9a.  Suppose 
the  sensor  is  located  at  the  origin,  while  the  vortex  is  centered  at  x  =  xo,y  =  0.  The 
points  of  intersection  of  the  two  circles  are  jc,-  ±  y, ,  where 


r]-rl  ■¥  xl 


2Xn 


(2.77) 


and 


yi  = 


r  —  r  +  X 
'  S  f  y  *  2 


2Xr, 


(2.78) 


Let  Ts  be  the  area  within  the  sensor  circle  on  the  vortex  side  of  the  chord  connecting  the 
points  of  intersection  (see  Figure  2.9b)  and  be  the  area  within  the  vortex  circle  on  the 
sensor  side  of  the  same  chord.  As  shown  in  Figure  2.9b,  the  area.  A,  enclosed  by  a  circle, 
a  chord,  and  the  centerline  perpendicular  to  the  chord  is 


A  =  xVr^-jc^  +  r^sin  * 


2L 

r 


(2.79) 


Therefore,  when  Xi  >  0, 

Ts  =  -  Xi  Vr?  -  x}  -  r,  sin“* 

When  Xi  <  0,  the  last  case  shown  in  Figure  2.9a, 


Ts  =  ^r]  +  {-Xi)  -  (-Xi)^  +  sin  ’ 


-Xi 


(2.80) 


(2.81) 


which  simplifies  to  the  identical  expression  for  Xi  >  0  in  Eq.  2.80.  Likewise,  for  Ty,  when 

Xq  >  Xi, 

'  XQ-Xi 


Ty  =  ^rl  -  (XQ-Xi)  -  {XQ-Xif  -  rl  sin  * 


(2.82) 


For  the  one  case  where  xq  <  x,-. 


Ty  =  +  (Xi-xo)  xjrl  -  {x-xq)'^  +  rl  sin  ’ 

which  reduces  to  the  Xo>Xi  case.  In  all  cases,  the  total  overlap  area  is  Ts+Ty. 


Xi-Xo 


(2.83) 
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Xo>Xi 


I  ^  >  0 


^  Xo<  ^ 


^  ^  >^  <  0 


XiO  Xq 


Figure  2.9(a).  Four  geometric  cases  of  a  vortex  overlapping  the  sensor  head.  The  cases  are  dis¬ 
tinguished  by  the  relative  size  of  the  circles  and  the  amount  of  overlap. 


Figure  2.9(b).  The  area,  A,  enclosed  by  a  circle,  a  chord,  and  the  centerline  perpendicular  to 
the  chord  is  x  -x^  +  T^sin“^(A:/r).  The  overlap  area  of  the  sensing  area  and  a  vortex  of  the 
same  size  is  Ty  +  T,,  where  T,  =  rt  r//2  -  [Xi^rj  -xj  +  xhirr^x/r,)],  where  x,-  is  the  horizon¬ 
tal  distance  from  the  center  of  the  sensor  head  to  the  point  of  intersection.  Likewise, 

Ty=‘K  ry^/2  -  [(xo  -Xi)  Vtv  -(to  -x,)^  +  (xq  -x,)^  sin“^[(xo  -x,)/rv],  where  xq  is  the  center 
of  the  vortex  relative  to  the  center  of  the  sensor. 
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A  MATLAB  program  was  written  which  computes  the  average  effect  of  a  vortex 
centered  at  all  possible  locations  on  an  equi-dimensional  grid  with  separations  of  rg/NV, 
where  NV  is  the  number  of  vortices  in  the  interval  (0,  r^+r^).  An  example  of  the  posi¬ 
tions  occupied  in  the  first  quadrant  when  NV  is  18  and  is  shown  in  Figure  2.10. 

When  Ky  «  rg,  for  example,  the  response  will  be  about  r^/r^  minus  the  effect  induced  by 
some  of  the  vortices  being  on  the  rim  of  the  sensor  so  that  only  a  fraction  of  their  area 
overlaps  the  sensor’s  area.  Because  the  area  expressions  assume  the  vortex  is  centered  on 
the  X  axis,  a  vortex  at  each  quantized  location  in  jc  and  y  is  expressed  in  a  rotated  coordi¬ 
nate  system  such  that  its  center  is  located  on  the  x  axis.  The  overlapping  area  is  com¬ 
puted  for  each  vortex  using  Tg  +  Ty  as  expressed  above.  Owing  to  symmetry,  it  is 
necessary  to  average  only  over  a  45°  section  of  the  sensor.  This  was  iterated  for  a  range 
of  vortex  sizes,  and  the  result  is  shown  in  Figure  2.11. 

It  is  possible  to  check  the  algorithm  against  an  analytic  solution  for  the  special  case 
where  the  vortex  radius  equals  that  of  the  sensor.  As  two  circles  of  the  same  radius,  R, 
pass  over  one  another,  the  total  overlap  area.  A,  is  twice  the  overlap  area  of  one  circle  on 
the  other. 
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-1.5  -1  -0.5  0  0.5  1  1.5 

X 

Figure  2.10.  An  example  of  one  initial  condition  for  the  MATLAB  program  that  calculates  the 
response  of  the  sensor  to  different  spatial  distributions  of  vorticity.  In  this  example,  the  positions 
and  sizes  of  the  analyzed  vortices  are  shown  for  the  first  quadrant  only,  where  NV  =  IS  and 
fix  =  5y  =  Vg/NV. 
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Figure  2.11.  Theoretical  spatially  averaged  amplitude  response  of  the  EMVM  sensor  to  an  iso¬ 
lated  vortex,  where  is  the  ratio  of  the  radius  of  the  vortex  to  the  radius  of  the  sensing  area. 
For  vortices  much  larger  than  the  sensor,  the  response  approaches  unity.  For  vortices  much 
smaller  than  the  sensor,  the  response  approaches  (ry/r^)^. 


A=2Ty  =  kR^  -  2r  -r^  - 


2R^  sin"’ 


(2.84) 


The  average  area  overlap.  A,  can  be  determined  by  integrating  the  above  expression  over 
the  grid  of  possible  vortex  positions  resulting  in  overlap  (in  this  case,  a  circle  of  radius 
2R)  and  dividing  by  the  area  of  that  grid.  That  is, 

_  2R  2R  2R2R 

A  =  - — -  I  J  A(x,y)dxdy  =  — -  J  \A{x,y)dxdy  .  (2.85) 

AnR  -2R-2R  t^R  0  0 


Transforming  to  polar  coordinates, 

JL 

2  2R 

A  =  — —  j  I A  (r,  Q)rdrdQ , 

tiR  0  0 


(2.86) 


where  ©  is  the  azimuthal  angle.  However,  since  A  (r,  0)=A  (r). 
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After  integration. 


2R 


IK  n 


(2.87) 


A  =  7iR^  - 


\r{R^-r^)^-^  +  ■^R^r'^R^-r^  +  "“sin  '(“) 
4  8  8  A 


2R 


(2.88) 


-4 


ri  ^ 

sin  ^ 

_r. 

+  -^R^-r^ 

2  4 

R 

4 

2R 


J  0 


=  Tzr'^  -  ~TiR^  -  ~7c/?^  =  ~nR^  . 

4  2  4 

Normalized  to  the  sensor  head  radius,  nR^  , 

-At  =  0.25 .  (2.89) 

iiR^ 

A  value  of  0.25  agrees  with  Figure  2.11  for  =  1;  that  is,  the  expected  response 
for  a  unit  vortex  with  its  center  randomly  placed  in  and  around  the  sensor  circle  is 
0.25  s“* .  Of  course,  when  the  vortex  and  sensing  area  completely  coincide,  the  response 
is  1 ;  however,  the  response  goes  to  zero  as  the  displacement  between  the  centers  of  the 
two  circles  goes  to  2R. 

2.6  Relationships  Between  Electric  Field  and  Velocity  Components  for  2-D  Sensor 

A  2-D  EMVM  sensor  head  consists  of  five  electrodes  connected  by  tubing  to  five 
locations  above  a  permanent  magnet.  The  geometry  is  shown  in  Figure  2.12.  The  volt¬ 
ages  induced  by  flow  between  the  center  electrode  and  each  outer  electrode  are  mea¬ 
sured.  The  four  independent  voltages  allow  the  determination  of  an  equal  number  of 


Figure  2.12.  Electrode  geometry.  Electrode  sites  are  numbered 
1-4  starting  from  the  upper  left  facing  the  sensor  with  the  mag¬ 
netic  field  component  pointing  out  of  the  page. 
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flow  characteristics.  One  can  choose  which  quantities  to  compute,  such  as  velocity  and 
two  gradient  terms.  The  potential  gradient  in  the  water  induced  by  a  flow  U  through  a 
magnetic  field  B  is 

V(t>  =  UxB-J/a,  (2.90) 

where  ^  is  the  electric  potential,  J  is  the  electric  current  density,  and  a  is  the  electrical 
conductivity  of  the  seawater. 

Consider  a  right-hand  coordinate  system  with  x  and  U  positive  to  the  left,  y  and  V 
positive  into  the  figure,  and  z  and  W  positive  up.  (This  z  is  not  the  z  used  in  altitude  but 
one  referenced  to  sensor  head.)  For  a  magnetic  field  of  B  (positive  in  direction  out  of  the 
sensor,  i.e.,  out  of  the  figure)  and  velocity  components  of  U  and  W,  the  potentials 
(e  1 , . . . ,  e  4)  at  the  ends  of  the  electrode  lines  are 

—e  I  =  BU I  -  BWi  _  _  Al  (2.91a) 

/  a  a 

—€2  =  -BU2  +  BW2  +  —  +  ^  (2.91b) 

/  a  a 

— eg  =Bt/3 +BW3  +  — (2.91c) 
I  a  a 

=  -5[/4  -BW4  -  —  +  —  ,  (2.91d) 

I  o  a 

where  /  is  the  magnitude  of  the  distance  vector  I;  between  the  center  and  ends  of  indivi¬ 
dual  electrode  tubes,  e,-  =  V<t>-1,-,  B  is  the  magnitude  of  the  magnetic  field  (assumed  to  be 
uniform  over  the  sensor),  Ui  and  W,-  are  the  velocity  components  across  each  electrode 
line,  and  Jxi  and  are  the  electric  current  density  components  along  each  electrode  line. 

A  difficulty  is  that  no  independent  measurement  of  J/a  is  obtained.  That  is,  measur¬ 
ing  Vtj)  along  a  given  vector  line  (1)  and  knowing  B  is  not  enough  to  uniquely  determine 
U  (Shercliffe,  1962).  What  is  determined  by  the  potential  measurements  is  the  difference 
between  the  source  term  U  x  B  and  the  response  term  J/a.  Because  the  induction  is  a 
linear  problem,  the  response  is  proportional  to  the  source.  If  the  full  magnetic  field  and 
flow  are  known,  this  response  function  can  be  computed.  Usually,  it  is  determined 
empirically  from  calibrations  in  a  tow  tank  or  flume.  Both  methods  have  been  applied  in 
this  work.  The  interpretation  equation  used  is 

V(j)-l  =  A(l)  =  CUxB-l.  (2.92) 

From  theory  and  experiment,  we  determined  that  the  head  factor,  C  (denoted  as  Hd  ear¬ 
lier),  is  equal  to  0.5  for  flow  perpendicular  to  the  applied  magnetic  field. 
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Now,  if  we  write  Eq.  2.92  along  the  z  axis, 


Similarly,  for  the  x  axis. 


Jz 

BU-  —  =  CBU 

a 


BW-  —  =  CBW. 

G 


Hence,  the  combinations  of  e  j  and  e  3  yield 

V2  . 


(e3-^l)=(W3  +  Wi)  +  (t/3-[/i) 


(2.93) 


(2.94) 


(2.95) 


-(e3+^l)=(f/3  +  t/l)  +  (1^3-Wi), 


(2.96) 


where  t/,-  and  W,-  are  spatially  averaged  velocity  components  along  the  ith  electrode  line. 
The  combinations  of  ^2  and  yield 

/  \  yxTT  TTT  \  /•XT  XT  \ 


■ie2-e4)  =  iW2  +  W4)-(U2-U4) 


(2.97) 


-(^2  +  C4)  =  -  (f/2  +  U4)  +  {W2  -  W4)  . 


(2.98) 


The  four  independent  voltages  can  be  combined  to  determine  four  flow  characteris¬ 
tics.  The  four  that  are  computed  are  two  velocity  components,  one  vorticity  component, 
and  a  strain  component.  The  relations  are 


^/2,  Ui+U2  +  U-i  +  U4 

(^  ..  —  U  , 

4CBI  1  ^  ^  4 


(2.99) 


^  W1+IT2  +  W3  +  W4 

-^^(ei -^2-^3 +  ^4)  =  =W ,  (2.100) 


ei  +  g2~^3~^4  _ 

2CBl^ 

‘/2(Wi  +  W3)  -  y2iW2  +  W4)  -  YijUi  +  U4)  -  V2{U2  +  U3) 

iaIi 

^  dw  _du  _  dw 

dz  dx  dz 


(2.101) 
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ei  +  €2  +  ej  +  €4 

Bf 

VtjUl  +  U2)  -  V2{U2  +  U4)  -  ‘/2(Wl  +  W4)  +  ’/2(W2  +  W3) 

/W2 


(2.102) 


du  _  aw 

"  az  a;c  ■ 

These  relations  assume  that  the  same  velocity  components  estimated  at  the  intersec¬ 
tion  of  orthogonal  lines  are  equal.  The  following  relations  were  used  in  the  Eqs.  2.99  and 
2.100; 


(Wi  +W2)-(W2  +  W4)  =  0 


(2.103) 


(t/l  +f/2)-(f/3  +  f/4)  =  0. 


(2.104) 


In  addition,  the  2-D  form  of  the  continuity  equation  was  invoked  in  Eq.  2.101. 
There  is  a  subtle  point  that  should  be  made  about  the  absence  of  the  head  factor  in  the 
expression  for  vorticity:  There  is  no  head  factor  for  vorticity;  vorticity  does  not  induce 
any  electric  currents,  so  C  =  1  (as  shown  analytically  in  Section  2.4.2).  The  basis  for  this 
can  be  seen  by  summing  the  expressions  for  the  individual  electrode  voltages  in  terms  of 
velocity  components  and  electric  currents.  The  sum  of  the  latter  forms  an  expression  of 
the  form  V-J/a,  which  is  zero.  This  fact  would  have  been  clear  if  the  individual  terms 
such  as  Jxi  and  were  retained.  This  was  not  done  in  the  interest  of  simplifying  the 
expressions. 
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3.  EARLY  INSTRUMENTS 

Early  EMVMs  used  a  synchronous  demodulation  scheme  to  separate  the  signal  from 
the  noise  by  modulating  the  magnetic  field  and  applying  narrow  frequency-band  signal 
processing  to  the  induced  voltages.  The  magnetic  field  was  changed  in  a  sinusoidal  or 
square-wave  modulation  from  positive  to  negative  values  at  a  frequency  that  was  higher 
than  the  dominate  electrode  drift.  The  electrode  voltages  were  synchronously  demodu¬ 
lated  at  the  magnetic  chopping  rate.  This  section  describes  the  first  four  versions  of  the 
EMVM. 

3.1  First  EMVM 

The  initial  EMVM  was  built  in  1982.  The  electronics  were  based  on  a  modified  cir¬ 
cuit  board  from  an  expendable  current  profiler,  or  XCP  (Sanford  et  al.,  1982).  A  block 
diagram  of  the  electronics  and  processing  is  shown  in  Figure  3.1.  The  electrode  voltages 
were  broadband  amplified  with  a  gain,  G  i ,  of  250;  they  were  then  combined,  high-pass 
filtered,  and  further  amplified  by  a  gain,  G2,  of  200.  This  approach  and  the  subsequent 
voltage-to-frequency  (V/F)  conversions  paralleled  those  used  in  the  XCP.  However,  the 
modified  electronics  were  susceptible  to  ambient  60-Hz  noise  related  to  grounding  of  test 
and  auxiliary  equipment.  This  effect  was  minimized  by  isolating  the  electronics  through 
the  use  of  a  fiber-optic  data  link  and  battery  power. 

The  simple  uniaxial  sensor  (Figure  3.2)  consisted  of  silver/silver-chloride  electrodes 
surrounded  by  agar.  The  length  of  the  electrode  arms  on  each  side  was  approximately 
2.5  cm.  The  test  facilities  are  shown  in  Figure  3.3.  External  Helmholtz  coils  were  used 
to  supply  a  known  vertical  magnetic  field  at  a  specified  frequency.  Vorticity  was  gen¬ 
erated  by  an  impeller-driven,  rigid-lid  tank.  The  sensor  assembly  was  sealed  into  the 
tank,  with  the  electrical  connections  extending  up  through  the  rigid  lid.  Though  provid¬ 
ing  only  an  approximation  to  solid-body  rotation,  the  impeller  system  was  judged  ade¬ 
quate  for  the  initial  phase  of  investigation.  Because  the  electronic  circuitry  was  sensitive 
to  variations  in  vorticity,  not  steady-state  vorticity,  the  impellor  was  driven  by  a 
variable-speed  motor.  By  varying  the  speed  of  the  motor  in  a  controlled  way,  known 
variations  in  the  angular  rotation  rate  were  produced.  These  tests  were  performed  in 
mid-September  1982. 

This  initial  phase  verified  the  concept  of  measuring  the  vorticity-induced  electric 
field  in  the  manner  envisioned.  However,  the  detected  signals  were  two  to  three  times 
less  than  the  expected  values.  This  deviation  may  have  been  due  to  inadequacies  of  the 
simple  sensor.  The  tubular  shape  of  the  sensor  disrupted  both  the  flow  field  and  the 
resulting  electric  field.  The  magnetic  field  was  not  uniform  throughout  the  entire  flow 
volume.  The  effects  of  the  insulating  tank  walls  were  not  understood.  Finally,  the  fluid 
was  not  in  truly  solid-body  rotation,  and  the  vorticity  was  therefore  difficult  to  quantify. 
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Figure  3.1.  Block  diagram  of  the  electronics  for  the  first  version  of  the  EMVM. 
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The  sensitivity  of  the  first  instrument  was  limited  by  the  strength  of  the  magnetic 
field  generated  by  the  Helmholtz  coils.  The  maximum  magnetic  field  available  was  0.3 
to  0.4  X  10“^  T.  Moreover,  the  magnetic  field  was  external  and  not  part  of  a  more  com¬ 
pact  sensor.  A  projection  of  the  sensitivity  was  used  as  the  basis  for  the  design 
refinements  incorporated  in  later  versions.  Detection  of  a  vorticity  signal  of  10“^  s“^ 
(corresponding  to  a  rotation  period  of  2  min)  was  thought  possible  using  a  magnetic  field 
of  10'^  T. 

3.2  Second  EMVM 

The  next  stage  of  development  was  to  increase  the  sensitivity  of  the  vorticity  meter 
and  to  incorporate  the  source  of  the  magnetic  field  into  the  sensor  head.  The  steps 
included  improving  the  electronics  and  processing  scheme  and  making  a  three- 
dimensional  measurement.  A  test  facility  was  built  to  generate  low-level  vorticity  in  the 
form  of  solid-body  rotation.  This  second  version  of  the  vorticity  sensor  was  completed 
on  1  December  1982. 

A  major  problem  with  using  a  rapidly  changing  magnetic  field  is  that  it  induces 
large  voltages.  A  time-varying  magnetic  field  acting  on  a  residual  loop  area  produces 
inductive  coupling.  In  the  first  version  of  the  vorticity  meter,  the  magnetic  field  was  dis¬ 
tributed  over  a  relatively  large  volume  by  the  Helmholtz  coils  (see  Figure  3.3).  In  the 
second  version  (Figure  3.4),  the  magnetic  field  was  provided  by  a  coil  in  the  sensor  itself. 
This  confined  the  influence  of  the  magnetic  field  to  a  smaller  volume,  which  facilitated 
isolation  of  the  electronics. 


Figure  3.4.  Second  version  of  the  EMVM.  The  magnetic  field  is  provided  by  a  coil  in  the  sen¬ 
sor  itself. 
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Most  of  the  residual  loop  area  was  in  the  electrical  connections  between  the  elec¬ 
trodes  and  the  electronics.  However,  the  magnitude  of  the  undesired  signal  was  rela¬ 
tively  low  because  of  the  low  magnetic  field  and  the  large  signal  provided  by  the  vorti- 
city.  At  the  large  magnetic  field  level  needed  to  detect  small  vorticity,  such  a  level  of 
continuous  induction  noise  is  intolerable. 

This  sensor  was  therefore  improved  and  expanded  to  make  a  three-axis  measure¬ 
ment,  as  shown  in  Figure  3.5.  The  radial  symmetry  minimized  interference  with  the  elec¬ 
tric  field.  The  disc  shape  was  intended  to  give  the  least  disturbance  to  the  flow  field.  The 
coil  axis  was  fixed  parallel  to  that  of  the  z  coordinate,  so  all  measurements  were  of  the 
z  component  of  vorticity. 

To  minimize  the  problem  of  electronic  noise  caused  by  induction  coupling,  the  drive 
coil  on  the  sensor  was  modulated  with  a  square-wave  signal.  This  confined  the  induction 
coupling  to  short  noise  spikes  (delta  functions)  superimposed  on  an  otherwise  low  noise 
output.  A  blanking  procedure  was  then  used  to  eliminate  these  noise  spikes  from  the  out¬ 
put.  The  scheme  is  illustrated  in  Figure  3.6.  The  track-and-hold  circuit  was  driven  by 
the  blanking  generator,  and  a  lock-in  amplifier  was  used  to  sample  the  output  synchro¬ 
nously  with  the  magnetic  field  modulation. 


Z-AXIS 


3-AXIS 

SENSOR 


Figure  3.5.  Tri-axial  configuration  of  the  second  version  of  the  EMVM. 
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1st  stage:  gain*  100 


Figure  3.6.  Block  diagram  of  the  electronics  for  the  second  version  of  the  EMVM. 

Because  of  our  increased  sensitivity,  we  needed  a  new  facility  to  impart  a  known 
amount  of  vorticity  into  the  water  surrounding  the  sensor.  In  February  1983,  we  con¬ 
structed  a  tank  that  could  be  rotated  slowly  on  a  lazy-Susan  base.  After  some  time,  this 
rotating  tank  would  spin  the  water  up  into  solid-body  rotation.  Blades  were  attached  to 
the  walls  of  the  tank  to  speed  up  this  process.  The  instrument  was  mounted  so  that  the 
sensor  was  in  the  center  of  the  tank,  where  it  would  be  surrounded  by  water  that  was  in 
solid-body  rotation.  The  period  of  rotation  of  the  water  was  the  same  as  that  of  the  tank 
and  could  be  varied  over  a  wide  range.  This  facility,  shown  in  Figure  3.7,  had  a  lower 
vorticity  limit  of  0. 1  s“^ . 

In  July  1983,  we  improved  the  sensitivity  by  changing  the  input  preamplifiers  from 
Precision  Monolithics  OP-lOs  to  PMI  OP-227s.  This  decreased  the  input  noise  from 
about  10  nV/^fc  to  about  3.5  nV/^ST.  We  also  tried  driving  the  coil  at  higher  currents 
to  increase  the  strength  of  the  magnetic  field.  The  drive  levels  for  the  three-axis  disc  sen¬ 
sor  were  then  of  the  order  of  30  G  peak  to  peak.  Tank  tests  on  this  sensor  indicated 
development  had  progressed  to  the  point  where  we  had  exceeded  our  sensitivity  goal  of 
0. 1  s”* .  The  measured  noise  level  in  still  water  was  0.02  s“^  rms  over  the  band 
0.01-1  Hz.  We  achieved  this  level  by  decreasing  the  noise  level  of  the  electronics  to 
12  nV  rms  over  a  1-Hz  bandwidth.  This  represents  approximately  5  nV  rms  per  amplifier 
for  the  six  amplifiers  used  (one  for  each  electrode),  which  equals  the  manufacturer’s 
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Figure  3.7,  Test  facility  used  for  the  second  version  of  the  EMVM. 

specification  for  this  device.  This  noise  level  was  achieved  only  when  the  impedance 
connected  across  the  amplifiers  inputs  was  1000  Q  or  less.  To  obtain  this  low  impedance, 
we  used  silver/silver-chloride  electrodes  and  relatively  large  connecting  ports  out  to  the 
point  of  measurement.  These  electrodes  exhibited  very  low  impedance  for  their  size  and 
had  stable  and  small  self-potentials. 

3.3  Third  EMVM 

In  August  1983,  we  constructed  an  underwater  housing  for  the  EMVM  electronics 
using  thick-walled  PVC  pipe.  The  sensor  was  attached  to  one  of  the  endcaps,  and  the 
wiring  for  signals  and  power  was  brought  out  through  an  underwater  connector.  In  Sep¬ 
tember  1983,  we  constructed  a  new  deck  box  that  could  interface  to  a  1000-ft  underwater 
tow  cable.  The  new  deck  box  provided  a  16-bit  parallel  interface  to  an  HP-9845  com¬ 
puter.  A  computer  program  was  written  to  perform  the  synchronous  demodulation  and 
graphics  display  of  the  data. 

A  block  diagram  of  the  drive  and  signal  processing  electronics  is  shown  in  Figure 
3.8.  The  magnetic  field  was  modulated  and  applied  at  a  15-Hz  rate.  The  modulation  was 
square  wave,  and  the  voltage  spikes  induced  were  ignored.  The  15-Hz  commutation  fre¬ 
quency  was  chosen  because  it  was  a  subharmonic  of  60  Hz,  which  greatly  decreased  our 
sensitivity  to  power-line  pickup.  It  was  also  high  enough  so  that  l/f  noise  from  the 
amplifiers  was  not  a  problem. 
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Figure  3.8.  Block  diagram  of  the  electronics  for  the  third  version  of  the  EMVM. 
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Our  deck  box  performed  two  basic  functions:  driving  the  coil  at  the  15-Hz  commu¬ 
tation  frequency  and  synchronously  demodulating  the  sensor  signals.  The  synchronous 
demodulation  was  performed  by  counters  that  were  allowed  to  count  up  during  the  part 
of  the  coil  drive  cycle  after  the  switching  transients  had  died  out.  The  counters  were  then 
read  and  reset  during  the  blanking  periods  when  the  switching  transient  was  present. 
Several  different  counting  schemes  were  switch  selectable  from  the  front  panel.  The  best 
results  were  obtained  when  one  coil  was  driven  at  15  Hz,  the  counters  counted  for  1/60  s, 
and  the  blanking  interval  was  1/60  s. 

The  first  field  tests  of  the  integrated  system  were  conducted  from  the  APL  acoustics 
barge  in  October  1983.  For  these  tests,  the  sensor  and  electronics  were  attached  to  the 
ram  mounted  in  the  center  of  the  vessel.  However,  the  noise  level  was  very  large 
because  of  60-Hz  interference.  (The  ac  current  distribution  on  the  vessel  appears  to  have 
numerous  leaks  into  the  surrounding  seawater.)  As  a  result,  nanovolt-level  measure¬ 
ments  were  not  easily  performed.  An  attempt  was  made  to  avoid  60-Hz  interference  by 
installing  the  sensor  and  electronics  on  a  fiberglass  boat  towed  by  the  barge.  Power  and 
data  were  conveyed  by  a  multiconductor  cable.  The  60-Hz  noise  was  much  less,  but  we 
observed  high  noise  levels  due  to  the  vigorous  motion  of  the  light  boat  in  the  wake  of  the 
barge.  Because  our  intent  was  to  determine  more  about  the  low-noise  performance,  we 
decided  to  install  the  system  on  a  platform  on  the  bottom  of  Puget  Sound  for  the  noise 
tests.  At  the  same  time,  we  decided  to  develop  a  new  sensor  with  the  goal  of  reducing 
the  disturbance  of  the  flow  past  the  sensor. 

The  new  sensor  was  completed  on  9  January  1984.  This  sensor,  shown  in  Figure 
3.9,  consisted  of  electrodes  mounted  on  long  struts  extending  into  the  flow.  This  orienta¬ 
tion  presents  the  least  disturbance  to  the  flow,  assuming  the  sensor  is  aligned  with  the 
oncoming  velocity  as  on  a  towed  body.  The  magnetic  field  was  provided  by  Helmholtz 
coils  surrounding  the  sensing  volume.  Although  this  configuration  induced  less  interfer¬ 
ence  with  the  ambient  flow,  the  magnetic  field  strength  was  less  since  the  coils  were 
rather  far  away  from  the  electrodes. 

In  February  1984,  we  installed  the  new  system  on  a  tripod  platform  off  Pier  C  in 
Shilshole  Bay  in  Seattle  and  collected  several  hours  of  observations.  The  goals  of  the 
experiment  were  to  determine  the  in  situ  noise  of  the  system  and  to  verify  the  calibration. 
The  equipment  is  shown  in  Figure  3.10.  A  plastic  garbage  can  was  used  to  shroud  the 
sensor  for  the  noise  tests;  it  did  not  interfere  with  the  sensor  but  excluded  any  velocity 
and  vorticity  signals  from  entering  the  sensing  volume.  An  example  of  the  data  obtained 
with  and  without  the  shroud  is  presented  in  Figure  3.11.  The  variables  presented  are, 
from  top  to  bottom,  vorticity,  the  x  and  y  components  of  velocity,  the  x  and  y  components 
of  acceleration  (from  a  separate  sensor),  and  the  carrier  frequencies  of  the  five  FM  data 
channels.  The  velocity  measurements  are  obtained  from  the  same  information  used  to 
compute  vorticity,  except  that  the  electrode  differences  are  added,  rather  than 
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Figure  3.9.  Third  version  of  the  EMVM.  This  version  is  a  three-axis  sensor  with  Helmholtz 
coils  for  field  modulation. 
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Figure  3.10.  Tripod-mounted  EMVM  (three-axis  open  sensor). 
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Figure  3.11.  Data  taken  with  tripod-mounted  three-axis  open  sensor  with  and  without  the 
shroud. 
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differenced,  and  expressed  as  velocity.  The  accelerations  were  monitored  to  ensure  that 
the  platform  was  not  being  moved  or  disturbed  by  the  flow.  This  figure  illustrates  the 
performance  in  the  field.  The  noise  level  is  about  3  x  10~^  s”^ .  Note  that  the  vorticity 
level  seemed  to  be  about  the  same  with  or  without  the  shroud.  Note  also  that  the  x  velo¬ 
city  channel  observes  some  disturbances  that  are  not  seen  in  vorticity.  This  is  one  of  the 
major  benefits  of  this  sensor — ^it  responds  to  vorticity,  not  just  some  isolated  velocity  dis¬ 
turbance.  Figure  3.12  shows  another  example  of  vorticity  measurements  from  the  plat¬ 
form,  in  this  case,  with  a  diver  creating  a  disturbance  upstream  from  the  sensor.  There  is 
a  sizable  vorticity  signature,  whereas  there  is  little  or  no  disturbance  in  the  velocity  chan¬ 
nels.  Compare  this  example  with  the  previous  figure. 

It  is  not  possible  to  generate  vorticity  structure  of  known  size  and  strength.  Having 
divers  swim  upstream  of  the  sensor  hardly  constitutes  a  calibration  of  the  system.  We 
therefore  mounted  the  sensor  and  electronics  in  such  a  way  that  both  could  be  rotated 
through  a  known  angle  (see  handle  in  Figure  3.10).  In  this  case,  the  time  integral  of  the 
vorticity  is  just  twice  the  rotation  angle.  An  example  of  slow  rotation  over  a  duration  of 
about  10  s  is  shown  in  Figure  3.13;  the  vorticity  is  shown  in  the  top  panel  and  the  velo¬ 
city  components  in  the  next  two  panels.  There  is  a  local  current  of  about  3  cm  s  ^  that 
changes  sign  on  the  y  current  sensor  as  the  unit  is  rotated.  The  time  integral  of  one  of 
these  rotations  is  displayed  in  Figure  3.14.  The  area  is  1.9  rad,  or  1 10°,  the  angle  through 
which  the  sensor  was  in  fact  rotated.  This  agrees  with  the  results  of  similar  laboratory 
experiments. 

3.4  Fourth  EMVM 

It  was  realized  that,  to  make  an  order  of  magnitude  improvement  in  the  sensor  sen¬ 
sitivity,  the  size  of  the  sensor  would  have  to  substantially  increase.  Figure  3.15  shows 
the  fourth  configuration  of  the  instrument  while  placed  in  a  support  frame  for  use  on  the 
seafloor.  The  EMVM  itself  was  just  under  2  m  long  and  about  0.5  m  in  diameter.  In  this 
version,  the  electronics  were  housed  in  a  larger  pressure  case  which  allowed  a  new 
high-power  coil  driver  to  operate  underwater.  In  addition,  a  new  data  transmission  tech¬ 
nique  was  used  to  send  the  data  over  the  underwater  cable.  The  underwater  unit  also 
grew  to  accommodate  additional  sensors  such  as  a  compass  and  pitch,  roll,  and  yaw-rate 
sensors.  The  magnetic  field  was  generated  by  Helmholtz  coils,  producing  a  field  flux 
density  of  50  G  (50  x  lO""^  tesla).  A  block  digram  of  the  electronics  is  shown  in  Figure 
3.16.  The  EMVM  was  previously  operated  by  means  of  a  HP-9845  computer.  New 
software  allowed  us  to  use  an  HP-9020  computer,  giving  faster  acquisition,  more  power¬ 
ful  real-time  processing,  and  a  color  graphics  display. 

First  tests  of  the  completed  instrument  were  performed  in  a  test  tank  at  APL  in  the 
spring  of  1985.  During  that  year,  numerous  laboratory  experiments  were  performed  to 
calibrate  the  EMVM  and  determine  the  noise  levels,  and  a  new  "paddle-type"  sensor  was 
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Figure  3.12.  Vorticity  measurements  with  a  diver  creating  a  disturbance  upstream  of  the 
EMVM  (three-axis  open  sensor). 
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Figure  3.13.  Vorticity  data  taken  while  the  third  version  of  the  EMVM  was  slowly  rotated 
(three-axis  open  sensor). 
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Figure  3.14.  Time  integral  of  one  rotation  of  the  tnpod-mounted  sensor 


Figure  3.15.  The  fourth  version  of  the  EMVM  mounted  in  a  support  frame  for  use  on  the 
seafloor. 
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constructed  during  the  summer  of  1985.  This  was  a  three-axis  sensor  with  the  measure¬ 
ment  ports  arranged  on  the  upper  and  lower  surfaces  of  a  thin  disc.  The  electrodes  and 
measurement  ports  for  the  vertical  axis  were  supported  on  "horns"  which  protruded  up 
from  the  rear  perimeter  of  the  disc. 

On  28  February  1986,  this  EMVM  was  field  tested  south  of  Port  Jefferson  in  Puget 
Sound  in  about  15  m  of  water.  The  deployment  served  as  an  engineering  test  of  the 
instrument  as  well  as  an  opportunity  to  obtain  a  time  series  of  vorticity  in  a  tidal  regime. 
The  instrument  was  placed  on  the  bottom  in  a  support  frame  with  the  help  of  divers.  The 
site  was  fairly  level,  with  a  smooth,  sandy  bottom.  While  in  the  support  frame,  the  vorti¬ 
city  and  current  sensors  were  2  m  above  the  bottom.  Apparent  vorticity  was  generated 
by  the  divers  slowly  rotating  the  EMVM,  providing  a  means  to  check  calibration  as  well 
as  to  determine  a  "zero"  for  the  current  measurements. 

After  the  divers  left  the  water,  a  3-hour  time  series  of  vertical  vorticity  and  horizon¬ 
tal  velocity  was  obtained  (Figure  3.17).  The  coordinates  were  rotated  so  that  velocities 
were  displayed  in  a  "flood-ebb"  reference.  The  EMVM  recorded  1.5  hours  on  either  side 
of  lower  low-water,  and  about  2  hours  before  and  1  hour  after  the  slack  water  before 
flood  tide.  The  vorticity  level  was  about  0.08  s“^  peak  to  peak  for  data  sampled  at  0.5-s 
intervals.  Just  before  slack  tide,  the  vorticity  levels  became  very  erratic,  in  part  because 
of  the  tidal  current  turning  and  approaching  the  EMVM  from  an  oblique  angle.  The 
divers  descended  again  and  rotated  the  EMVM  into  the  current.  The  vorticity  signal  lev¬ 
eled  off  at  about  0.40  s“^  peak  to  peak  after  the  adjustment. 

On  27  March  1986,  another  test  of  the  EMVM  was  made  in  Puget  Sound,  this  time 
2  n.mi.  north  of  Shilshole  Bay;  again  the  sensor  was  mounted  on  a  tripod  in  15  m  of 
water.  The  test  consisted  of  two  trials,  each  about  50  min  long.  The  first  was  at  the  time 
of  the  maximum  flood  tidal  current  in  the  middle  of  Puget  Sound  (0.25  m/s);  the  second 
was  during  a  decelerating  tidal  current  averaging  about  0.19  m/s.  These  values  were 
estimated  from  the  1986  NOAA  tidal  current  tables  for  the  Pacific  Coast  of  North  Amer¬ 
ica.  The  actual  velocities  measured  by  the  EMVM  were  about  6  cm/s.  The  two  trials 
compared  two  different  electromagnetic  field  sensors.  The  first,  the  same  type  used  in  the 
February  test,  was  disk  shaped  whereas  the  second  was  finger  shaped.  It  was  hoped  that 
the  finger-shaped  sensor  would  be  less  intrusive  than  the  disk-shaped  sensor  and  thus 
would  create  less  contamination. 

Further  testing  was  undertaken  in  the  laboratory  to  establish  reference  levels  of 
self-noise  under  the  quietest  conditions  possible.  For  the  laboratory  tests,  the  instrument 
was  suspended  so  that  the  magnetic  coils  and  the  sensors  were  immersed  in  a  salt-water- 
filled  plastic  tank.  The  tank  water  was  electrically  isolated  from  the  floor  and  all  sur¬ 
rounding  objects.  Velocity  and  vorticity  were  simulated  by  oscillating  the  instrument, 
using  a  small  motor  and  a  pivot  arm.  First,  the  arm  was  placed  to  the  side  of  the  EMVM 
to  simulate  the  vorticity  seen  by  the  EMVM  when  it  was  moving  in  pure  rotation  in  still 
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Figure  3.17.  Time  series  of  horizontal  velocity  and  vertical  vorticity  from  EMVM  deployment 
south  of  Point  Jefferson  in  Puget  Sound. 
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water.  Then  the  arm  was  moved  to  the  center  of  the  EMVM  to  simulate  the  velocity 
while  the  instrument  was  in  pure  translation  through  the  water.  In  a  quiescent  tank  with 
no  instrument  motion,  the  spectral  levels  of  vorticity  variance  were  10“^  s"^/Hz  (Figure 
3.18),  and  the  energy  levels  were  6  x  10”^  s”^/Hz.  When  the  instrument  was  set  up  for 
oscillation  at  a  specified  period,  the  level  of  the  vorticity  variance  spectra  remained  at 
10“^  s”^/Hz  (Figure  3.19)  except  for  the  spectral  band  around  the  oscillation  frequency. 


Figure  3.18.  Power  spectrum  of  vorticity  for  sensor  operated  in  still  water  in  an  electrically  iso¬ 
lated  test  tank.  A  noise  floor  of  10“^  (rad/s)^/Hz  is  reached  at  frequencies  above  5  x  10“^  Hz. 
Below  5  X  10“^  Hz,  the  noise  power  increases  100  times  for  every  decade  decrease  in  frequency. 


FREQUENCY,  Hz 

Figure  3.19.  Power  spectrum  of  vorticity  for  sensor  oscillated  in  pure  rotation  in  a  still-water 
test  tank  at  a  frequency  of  2  x  10“^  Hz.  The  noise  floor  remains  unchanged,  with  the  dominant 
signal  at  the  oscillation  frequency. 
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3.5  Experiences  with  Tow  Bodies 

Soon  after  the  first  successful  EMVM  tests  on  a  tripod  early  in  1984,  we  wanted  to 
install  the  instrument  on  a  vehicle  for  towing.  We  investigated  existing  tow  bodies  and 
found  two  that  might  meet  oiir  requirements.  One  was  the  system  used  by  the  Navy 
known  as  the  C  Mk  1.  The  other  was  the  commercial  version  of  the  Batfish,  known  as 
Seasoar.  The  former  could  be  had  on  loan;  the  latter  cost  over  $50,000.  Initially,  we 
sought  to  use  the  former;  ultimately,  we  purchased  the  latter. 

3.5.1  C  Mk  1  tow  body 

In  the  spring  of  1984,  we  received  a  C  Mk  1  towed  vehicle  on  loan  from  the  Naval 
Coastal  Systems  Center.  Figure  3.20  shows  the  unit  as  received.  The  depressor  wings 
have  a  total  wingspan  of  about  5  ft.  The  central  body  is  a  fiberglass  pressure  case  with  an 
outside  diameter  of  about  9  in.  The  bottom-following  electronics  with  pitch  and  roll  sta¬ 
bilization  utilize  about  one-half  of  the  available  internal  volume.  The  sidescan  sonar  sys- 


Figure  3.20.  C  Mk  1  tow  body  as  received  from  the  Naval  Coastal  Systems  Center  in  1984. 
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tern  takes  up  another  one-third  of  the  available  volume,  leaving  a  small  space  forward 
near  the  nose  for  an  auxiliary  payload.  The  internal  electronic  circuits  utilized  vacuum 
tubes  which  made  them  obsolete.  We  were  able  to  find  replacement  tubes,  so  we  used 
the  existing  electronics  rather  than  building  solid-state  replacements.  As  received,  the 
electronics  were  designed  to  fly  the  vehicle  at  a  constant  height  off  the  bottom.  The 
bottom-following  circuitry  was  modified  to  allow  the  vehicle  to  fly  at  a  controlled  depth, 
and  the  APL  winch  used  for  towed  targets  was  modified  to  accommodate  the  C  Mk  I’s 
tow  cable.  Several  test  tows  were  made  in  April  1984  in  Lake  Washington  to  evaluate 
the  dynamic  stability  of  the  vehicle.  Because  of  the  large  surface  area  of  the  depressor 
wings,  large  transient  loads  were  applied  to  the  winch,  cable,  and  wooden  deck  plates  of 
the  towing  vessel.  The  C  Mk  1  was  able  to  fly  at  a  constant  depth,  although  pitch  stabil¬ 
ity  was  obtained  only  under  fair  weather  conditions.  The  tow  characteristics  drifted 
slowly  over  time,  and  upon  retrieval  it  was  found  that  the  fiberglass  wings  were  filling 
with  water. 

3.5.2  GW  Fairing 

During  the  spring  of  1987,  we  designed  and  constructed  a  fiberglass  fairing  for  the 
vorticity  meter  which  allowed  the  whole  instrument  to  be  towed  underwater.  This  was 
done  rather  than  use  the  C  Mk  1  vehicle  because  the  C  Mk  1  would  have  been  unable  to 
accommodate  the  additional  volume  and  weight  of  EMVM  version  4.  Figure  3.21  shows 
the  instrument  out  of  the  water  hanging  from  the  towing  bridle.  The  orientation  of  the 


Figure  3.21.  Vorticity  meter  hanging  from  the  tow  bridle,  in  normal  operating  configuration. 
The  fiberglass  fairing  designed  by  G.  Welsh  (GW  fairing)  is  prominent. 
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instrument  shown  in  this  photograph  is  the  normal  operating  configuration.  The  pressure 
case  of  the  EMVM  is  in  essence  the  leading  edge  of  an  airfoil  which  is  towed  in  an 
upright  orientation  by  a  bridle  attached  to  the  the  upper  endcap  of  the  instrument.  The 
height  of  the  airfoil  is  36  in.  Above  and  below  this  wing  section  are  flat  flow-splitting 
plates  which  taper  from  a  width  of  16  in.  at  the  leading  edge  to  24  in.  at  the  trailing  edge. 
At  the  rear  are  adjustable  elevators  which  are  used  to  trim  the  tow  body  for  level  flight  at 
the  nominal  tow  speed.  The  electromagnet  coil  and  a^ay  of  sensor  measurement  tubes 
extend  below  the  lower  flow-splitting  plate  of  the  fairing.  The  first  test  tows  of  this 
instrument  were  accomplished  in  June  1987  in  Puget  Sound.  Minor  tweaking  of  the  bri¬ 
dle  attachment  point  and  elevator  angle  resulted  in  a  stable  towed  instrument.  However, 
the  roll  stability  was  greatly  affected  by  the  tow  speed.  At  tow  speeds  above  3  kn,  the 
instrument/fairing  combination  began  to  roll  back  and  forth  with  amplitudes  of  ±15°  or 
more  at  a  period  of  a  few  tens  of  seconds.  This  behavior  limited  the  tow  speed  to  less 
than  3  kn. 

3.5.3  Seasoar  tow  body 

The  Seasoar  vehicle  was  purchased  from  Chelsea  Instruments  Ltd.,  U.K.,  in  July 
1988  as  a  replacement  for  the  C  Mk  1.  Figure  3.22  shows  the  Seasoar  vehicle  as  shipped 
from  Chelsea  Instmments.  This  particular  vehicle  was  chosen  on  the  basis  of  a  com¬ 
parison  study  commissioned  by  the  Naval  Ocean  Research  and  Development  Activity 
which  showed  that  the  Seasoar  was  superior  in  several  categories  of  performance.  For 
our  application,  the  reported  low  vehicle  vibration  levels  and  low  yaw  rates  were  key 
specifications.  Because  of  the  high  cost  of  this  unit  and  its  foreign  manufacture,  the  time 
required  to  purchase  and  take  delivery  of  this  device  was  almost  1  year.  By  the  time  the 
Seasoar  was  received,  we  were  successfully  using  the  GW  fairing  and  obtaining  good 
results  with  our  field  observations  in  Puget  Sound. 

3.5.4  Comparison  of  Seasoar  and  GW  fairing  performance 

The  Seasoar  tow  body  was  tested  in  Lake  Washington  on  24  April  1991  for  1  hour 
and  again  on  25  April  for  2  hours.  Tow  speeds  ranged  from  3.5  kn  to  about  5  kn  max¬ 
imum.  Because  the  water  depth  was  at  most  200  ft,  only  150-175  ft  of  cable  were  payed 
out.  Enough  experience  was  gained  from  these  tests  to  warrant  testing  in  Puget  Sound  at 
higher  speeds  and  with  more  tow  cable  extended.  On  29  April  1991,  the  Seasoar  was 
deployed  from  APL’s  RfV  Miller  in  Puget  Sound  for  5.5  hours.  The  maximum  tow  speed 
was  7.5  kn,  and  the  performance  of  the  Seasoar  was  evaluated  with  various  amounts  of 
tow  cable  payed  out. 

The  fiberglass  fairing  designed  by  G.  Welsh  and  used  to  tow  the  EMVM  since  1987 
was  most  recently  used  during  16-18  April  1991  for  tows  through  the  wake  of  Blakely 
Rock.  The  tow  speed  found  most  effective  for  this  work  was  about  2.5  kn.  At  higher  tow 
speeds,  the  EMVM  pitched  down  such  that  the  flow-splitting  plate  was  flying  with  a  0°  to 
-5°  angle  of  attack.  This  evidently  caused  some  flow  separation  at  the  leading  edge,  as  a 
large  increase  in  random  noise  was  observed  in  the  vorticity  data. 
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Figure  3.22.  Seasoar  oceanographic  tow  body.  The  wing  angle  can  be  changed  dynamically 
via  an  internal  hydraulic  actuator.  The  actuator  extracts  its  power  from  the  passing  seawater 
using  a  small  propeller  assembly  coupled  to  a  hydraulic  pump.  This  system  allows  the  tow  body 
to  be  controlled  remotely:  it  can  be  ordered  to  dive,  surface,  or  fly  at  a  prescribed  depth. 
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To  evaluate  the  performance  of  the  towed  underwater  vehicles,  a  special  self- 
contained  sensor  and  data-acquisition  package  was  constructed  to  collect  data  on  vehicle 
dynamics.  This  system  consisted  of  a  pressure  vessel  6  in.  in  diameter  and  22  in.  long 
housing  a  suite  of  three  accelerometers  mounted  orthogonally  to  each  other,  one  three- 
axis  magnetometer,  one  two-axis  angular  rate  sensor  (providing  pitch  rate  and  yaw  rate), 
and  a  pressure  transducer  for  depth  measurement.  These  sensors  were  anti-alias  filtered 
at  5  Hz  and  sampled  at  10  Hz  with  a  10-bit  digitizer.  The  digital  data  were  buffered  by  a 
TattleTale  Model  5  microcomputer  system  and  sent  up  the  tow  cable  as  an  asynchronous 
bit-serial  data  stream  at  9600  baud.  The  serial  data  were  optically  coupled  at  each  end  of 
the  cable,  and  the  power  to  the  sensor  package  was  isolated  from  the  ship  power  by  dc-dc 
converters.  This  eliminated  the  dc  ground  loop  path  allowing  electric  currents  to  flow 
down  the  tow  cable,  through  the  sensor  package,  and  back  up  through  the  seawater  to  the 
towing  vessel. 

Given  the  sensitivity  of  the  fourth  EMVM  to  static  pitch  angles  and  the  inherent 
sensitivity  to  yaw  rate,  it  is  these  parameters,  along  with  roll  angle  stability,  that  have 
been  analyzed  for  both  vehicles.  Figure  3.23  shows  time  series  for  yaw  rate,  pitch  angle, 
and  roll  angle  for  the  Seasoar  while  it  was  being  towed  in  Puget  Sound  at  a  depth  of  70  m 
with  170  m  of  cable  winched  out.  The  added  vertical  lines  indicate  the  times  when  the 
tow  speed  was  changed.  The  tow  speed  started  at  3.5  kn  and  was  increased  in  0.5-kn 
increments  up  to  a  final  to  speed  of  6.5  kn.  At  6.5  kn,  the  Seasoar  vehicle  becomes  less 
stable,  exhibiting  larger  roll  excursions  of  ±4°.  Also  at  6.5  kn,  the  yaw  rate  becomes 
active,  indicating  the  vehicle  is  fishtailing  with  a  maximum  rate  of  0.05  s“^  Figure  3.24 
is  a  similar  time  series  for  the  GW  fairing  while  it  was  being  towed  at  a  depth  of  10  m 
with  about  30  m  of  cable  out.  The  GW  fairing  shows  much  less  roll  stability,  with  excur¬ 
sions  of  ±5°.  This  is  probably  due  to  the  tow  cable  attachment  point  being  located  at  the 
top  of  the  body,  about  0.5  m  above  the  center  of  mass  of  the  combined  tow  body  and  sen¬ 
sor  package.  The  pitch  stability  is  similar  to  that  of  Seasoar,  with  maximum  excursions 
of  ±2°. 

Figure  3.25  compares  Seasoar  and  GW  fairing  yaw-rate  spectra.  The  solid  line  is 
the  spectrum  for  data  taken  with  the  GW  fairing  towed  at  2.5  kn;  the  dotted  line  is  for 
data  taken  with  the  Seasoar  towed  at  3.5  kn.  Below  0.1  Hz,  the  Seasoar  vehicle  exhibits 
yaw  rates  an  order  of  magnitude  lower  than  those  for  the  GW  fairing.  The  GW  fairing 
seems  to  have  an  energy  peak  at  a  20-s  period.  For  frequencies  above  0.2  Hz,  both  tow 
bodies  have  similar  yaw  rates.  The  Seasoar  spectrum  is  quite  a  bit  noisier  at  high  fre¬ 
quencies  because  it  is  realized  from  a  single  2(X)-s  time  series.  The  GW  fairing  spectrum, 
on  the  other  hand,  is  the  average  of  17  spectra,  each  based  on  a  200-s  time  series.  It  is 
therefore  difficult  to  compare  the  two  spectra  at  the  high  frequency  end.  Figures  3.26 
and  3.27  are  similar  plots  for  the  pitch  and  roll  data.  Again,  the  solid  line  is  the  spectrum 
for  the  GW  fairing  and  the  dotted  line  is  for  the  Seasoar. 
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Figure  3,23.  Time  series  of  yaw  rate,  pitch  angle,  and  roll  angle  obtained  while  the  Seasoar 
was  being  towed  in  Puget  Sound  at  various  speeds.  Time  scale  is  in  two-digit  units  indicating, 
from  left  to  right,  year,  month,  day  -  hour,  minutes,  seconds. 
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Figure  3.24.  Time  series  of  yaw  rate,  pitch  angle,  and  roll  angle  of  the  EMVM  with  GW  fairing 
being  towed  at  2.5  kn.  Time  scale  as  in  Figure  2.23. 
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Figure  3.25.  Comparison  of  yaw-rate  spectra  for  the  Seasoar  tow  body  and  GW  fairing. 


Figure  3.26.  Comparison  of  pitch  spectra  for  the  Seasoar  tow  body  and  GW  fairing. 
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Figure  3.27.  Comparison  of  roll  spectra  for  the  Seasoar  tow  body  and  GW  fairing. 


Figure  3.28  shows  the  yaw -rate  spectrum  for  Seasoar  when  the  vehicle  was  being 
towed  at  6.5  kn.  Note  the  large-amplitude  yaw  oscillation  at  about  0.19  Hz.  This  yaw- 
rate  amplitude  would  make  the  Seasoar  unacceptable  for  use  with  the  EMVM  at  this  tow 
speed. 

In  conclusion,  the  Seasoar  outperforms  the  GW  fairing  at  low  tow  speeds  in  all 
three  categories  of  comparison  (yaw  rate,  pitch  angle,  and  roll  angle).  At  tow  speeds 
below  3.5  kn,  however,  the  Seasoar  vehicle  becomes  unstable  in  depth  keeping  owing  to 
insufficient  hydraulic  pressure  to  activate  the  wing.  At  tow  speeds  above  6kn,  the 
Seasoar  again  becomes  unstable  and  exhibits  excessive  yaw  rates  and  roll  rates.  Since 
the  instability  is  narrowband  at  0.19  Hz  and  this  frequency  is  fairly  low  compared  with 
the  resonant  frequencies  of  the  control  surfaces  on  the  Seasoar,  it  may  be  due  to  a 
dynamic  interaction  of  the  tow  body  with  the  tow  cable.  It  may  be  possible  to  damp  out 
this  yaw  oscillation  with  an  appropriate  cable  fairing  or  a  longer  length  of  cable  payed 
out.  In  its  present  configuration,  the  GW  fairing  cannot  be  trimmed  to  maintain  positive 
pitch  angle  for  tow  speeds  above  2.7  kn,  limiting  its  use  to  tow  speeds  of  2.5  kn  or  lower. 

3.5.5  Seasoar  problems  and  manufacturer  repairs  and  upgrades 

After  only  8  hours  of  towing,  the  Seasoar  became  sluggish  in  its  response  to  control 
signals  from  the  deck  box.  In  addition,  control  settings  that  had  been  previously  deter¬ 
mined  were  no  longer  appropriate.  We  decided  to  bring  the  Seasoar  back  up  on  deck  for 
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Figure  3.28.  Yaw-rate  spectrum  for  the  Seasoar  vehicle  towed  at  6.5  kn. 

a  visual  inspection.  At  first  glance  nothing  seemed  to  be  wrong;  however,  when  the  pro¬ 
peller  was  turned  by  hand,  we  noticed  an  excessive  amount  of  radial  play  in  the  rear 
bearing  for  the  propeller  shaft.  We  returned  the  unit  to  the  laboratory  for  a  closer  inspec¬ 
tion. 

We  removed  the  propeller  from  the  tail  shaft.  There  was  no  significant  play 
between  the  propeller  and  the  shaft.  The  shaft  did  have  lateral  bearing  play,  so  we 
decided  to  remove  the  hydraulic  unit  from  the  vehicle.  With  the  unit  removed  from  the 
vehicle,  we  drained  the  rear  tail-shaft  housing.  Approximately  600  ml  of  fluid  was 
drained  from  the  housing.  It  was  a  frothy  mixture  of  oil  and  salt  water  with  a  slightly 
greenish  tint.  Some  brownish  to  black  particles  settled  out  of  the  solution.  The  most 
obvious  problem  was  that  the  Teflon  tail  shaft  bearing  was  free  to  spin  inside  the  PVC 
housing.  The  bearing  measured  1.170  in.  long  with  an  outside  diameter  of 
0.740-0.745  in.  The  inside  diameter  of  the  bearing  measured  0.500-0.501  in.  Both  of 
the  lip  seals  showed  wear  and  tearing  on  their  inside  diameters.  After  the  lip  seals  were 
removed,  it  was  possible  to  measure  the  diameter  of  the  hole  that  receives  the  Teflon 
tail-shaft  bearing.  The  inside  diameter  of  the  hole  measured  0.750  in.  We  have  con¬ 
cluded  that  the  bearing  was  always  free  to  wobble  in  the  PVC  housing  while  the  tail  shaft 
was  spinning.  This  could  cause  the  lip  seals  to  fail  owing  to  the  large  amount  of  lateral 
shaft  movement. 
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In  further  examining  the  PVC  rear  housing,  we  noticed  that  there  was  cracking 
adjacent  to  two  of  the  helicoil  inserts.  One  of  these  cracks  extended  radially  inward 
across  the  0-ring  face  seal  surface.  This  was  another  potential  source  of  seawater  intru¬ 
sion  into  the  oil-filled  cavity.  The  0.132-in.  diameter  (cross  section)  0-ring  that  provided 
the  seal  between  the  PVC  housing  and  the  flanged  stainless  steel  endcap  was  examined 
and  found  to  be  in  perfect  condition,  as  was  the  0.120  in.  deep  by  0.156  in.  wide  O-ring 
groove  in  the  stainless  endcap.  The  design,  however,  allows  for  only  a  10%  compression 
of  the  O-ring.  O-ring  manufacturers  recommend  16%  to  23%  compression  for  a  static 
face-seal  gland.  The  O-ring  was  therefore  another  potential  source  of  seawater  intrusion. 

In  conclusion,  the  Seasoar  developed  serious  mechanical  problems  after  only 
8.5  hours  of  operation.  These  problems  were  repaired,  but  a  redesign  of  the  rear  bearing 
assembly  was  required  for  increased  reliability.  In  addition,  a  redesign  of  the  rear  face 
seal  between  the  PVC  housing  and  the  stainless  steel  endcap  was  recommended.  Other 
users  of  the  Seasoar  had  reported  similar  problems,  and  Chelsea  was  making 
modifications  to  increase  the  reliability  of  the  vehicle.  Moreover,  the  present  EMVM  was 
too  large  to  install  on  the  Seasoar. 

We  therefore  decided  to  postpone  further  tow-body  work  in  favor  of  focusing  on 
improving  sensor  performance.  This  decision  lead  to  construction  of  a  fifth  version  of 
the  EMVM. 
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4.  DESIGN,  CONSTRUCTION,  USE,  AND  EVALUATION  OF  THE  FIFTH 
EMVM 

We  sought  to  implement  a  number  of  needed  improvements  in  a  new  system  while 
retaining  as  much  of  the  main  electronics  as  possible.  The  design  concentrated  on  a  new 
sensor  unit,  but  one  physically  and  electrically  isolated  from  the  main  electronics.  Most 
of  what  was  learned  about  earlier  versions  was  incorporated  into  the  new  sensor,  elec¬ 
trode  array,  and  magnet. 

4.1.  Design  Goals 

The  principal  design  goals  were  to 

•  eliminate  the  electric  measurements  along  the  axis  of  the  magnetic  field 

•  produce  less  disturbance  and  be  more  symmetric  to  the  flow 

•  protect  the  electrodes  for  better  noise  performance 

•  construct  the  sensor  out  of  plastic 

•  separate  the  sensor  from  the  main  electronics 

•  isolate  the  sensor  head  with  optical- fiber  data  links  and  batteries 

•  keep  flow  away  from  the  magnet  except  in  the  region  of  measurement 

•  incorporate  a  permanent  magnet  into  the  sensor. 

One  important  change  was  the  decision  to  forgo  3-D  measurements.  The  electrodes 
along  the  axis  of  the  magnet  were  not  providing  an  important  signal  in  the  situations 
encountered,  such  as  in  the  laboratory  tests.  It  was  decided  to  concentrate  on  an  array  of 
five  electrodes  in  the  x-y  plane,  normal  to  the  principal  axis  of  the  magnet.  The  indivi¬ 
dual  voltages  between  each  outer  electrode  and  the  central  reference  electrode  were 
recorded.  This  arrangement  permitted  us  to  compute  velocity  components  normal  to  the 
magnet’s  axis  and  the  vorticity  aligned  with  the  axis.  One  reason  for  this  design  was  to 
eliminate  asymmetry  in  the  sensor.  Earlier  designs  had  structures  that  disturbed  the  flow 
and  might  have  produced  different  results,  depending  on  the  orientation  of  the  incoming 
flow.  In  seeking  to  achieve  symmetry  in  the  sensor,  we  were  responding  to  problems  that 
were  not  predominantly  from  the  lack  of  symmetry  but  from  a  tendency  for  the  V/F  con¬ 
verters  to  influence  each  other.  Nonetheless,  we  continued  to  improve  symmetry  even 
after  the  V/F  problem  was  identified  and  corrected  because  the  goal  seemed  appropriate 
on  general  principles. 

The  most  noticeable  design  change  was  making  the  measurements  above  a  smooth 
circular  plate  (see  Figure  4.1).  The  circular  plate  has  a  tapered  leading  edge  to  reduce 
turbulence  and  is  intended  to  slide  through  the  relative  flow  with  little  angle  of  attack. 
Ceramic  straws  are  used  to  connect  the  electrodes  to  the  point  where  the  potential  is  to  be 
sampled.  The  electrodes  are  enclosed  in  chambers  beneath  the  surface  of  the  plate. 
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Figure  4.1.  The  splitter  plate  and  ceramic  straws  of  the  fifth  version  of  the  EMVM. 


called  a  splitter  plate.  A  seawater- filled  ceramic  tube  in  a  screw-in  fitting  electrically 
connects  the  electrode  to  the  point  of  measurement.  The  electrode  chamber  and  ceramic 
straw  are  filled  with  a  solution  of  agar  and  seawater.  The  resulting  jello-like  substance  is 
electrically  conducting  and  avoids  problems  from  air  bubbles  entrapped  in  tubes  filled 
with  just  seawater.  The  encapsulation  of  the  electrodes  is  intended  to  damp  the  tempera¬ 
ture  fluctuations  and  eliminate  the  salinity  changes  experienced  by  the  electrodes.  With 
a  sensitivity  of  about  300  |iV/°C  between  similar  electrodes  at  different  temperatures,  it  is 
important  to  reduce  this  effect  in  the  bandwidth  of  interest. 

The  electrode  cups  and  straws  are  4  in.  apart  on  the  plate  and  arranged  in  the  form 
of  a  plus  sign.  The  straws  are  4.5  in.  long  and  of  the  smallest  diameter  possible  for  the 
flow  conditions.  The  splitter  plate,  electrode  plenum,  and  magnet  chamber  are  con¬ 
structed  of  plastic  to  reduce  electrical  noise  from  corrosion  currents. 

The  electromagnet  is  housed  behind  the  splitter  plate  and  electrode  plenum.  It  is 
placed  within  a  plastic  enclosure  to  contain  any  stray  currents  which  might  leak  from  the 
potted  magnet  winding  and  to  reduce  flow  directly  on  the  coil,  which  might  produce 
strong,  erratic  signals.  The  size,  power,  inductance,  and  amount  of  time  needed  for  the 
inductive  spike  to  decay  required  that  we  operate  at  lower  frequency.  For  the  most  part, 
this  version  was  operated  at  frequencies  of  0.5  to  1  Hz.  This  chamber  is  also  designed  to 
hold  the  permanent  magnet  constructed  of  barium  ferrite  ceramic.  It  is  open  to  the  sea¬ 
water.  The  magnets  are  kept  in  place  by  packing  with  glass  microbeads.  The  intention  is 
to  keep  the  magnet  from  moving  with  respect  to  the  electrodes  and  cables  and  thereby 
inducing  spurious  voltages. 
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Experiences  with  earlier  instruments  had  demonstrated  problems  with  ground  loops 
and  other  inadvertent  paths,  even  high-impedance  ones,  between  circuits.  It  was  decided 
that  the  way  to  guarantee  adequate  isolation  was  to  operate  the  electronics  on  batteries 
and  use  optical  data  links.  It  Was  impossible  to  power  the  electromagnet  on  batteries,  but 
it  was  possible  to  isolate  the  magnetic  coils  electrically  from  the  preamplifiers  and  other 
sensor  electronics.  Moreover,  this  scheme  allowed  the  sensor  to  be  located  separately 
from  the  main  electronics,  easing  difficulties  in  placing  the  sensor  on  carriages  and 
operating  it  in  a  test  tank. 

The  low-noise  preamplifiers,  V/F  converters,  optical-fiber  drivers,  and  batteries  are 
installed  in  a  short  electronics  rack.  This  rack  is  housed  in  an  aluminum  pressure  vessel 
that  is  of  the  same  material  and  diameter  as  the  main  electronics  tube  but  is  only  21  in. 
long.  This  chamber  contains  the  essential  electronics,  batteries,  and  optical-link  com¬ 
ponents  and  supports  through-hull  connectors  for  recharging  the  batteries  and  attaching 
the  optical  fibers.  The  data  links  and  magnet  power  are  connected  to  the  main  electron¬ 
ics  package  which  remains  outside  the  test  tank  in  the  laboratory  or  on  the  boat  when 
making  field  observations.  (See  Figure  4.2  for  a  block  diagram  of  the  electronics.) 

The  design  achieved  many  of  the  goals  but  presents  some  difficulties.  The  straws 
are  fragile  and  contribute  a  rather  high  impedance  to  the  circuits.  This  high  impedance 
results  in  significant  Johnson  noise  and  additional  voltages  resulting  from  the 
preamplifier  bias  current  and/or  current  noise.  However,  the  AD  743  amplifiers  have 
voltage  and  current  noise  characteristics  that  are  favorable  for  this  application. 

4.2  Construction 

4.2.1  Electrode  array 

In  previous  sensors,  the  electrodes  were  placed  close  to  the  measurement  points. 
This  kept  the  inter-electrode  impedances  low,  but  at  the  same  time  the  thermal  environ¬ 
ment  of  the  the  electrodes  was  not  well  controlled.  In  the  new  design,  each  electrode  is 
housed  in  a  thick-walled  plastic  housing  with  a  0.5-in.  inside  diameter  and  a  1-in.  outside 
diameter.  These  housings  are  threaded  on  one  end  and  have  an  0-ring  seal  so  that  they 
can  be  installed  and  removed  easily.  The  underwater  cable  connecting  the  electrode  to 
the  preamplifier  is  sealed  to  the  electrode  housing  with  polyurethane  potting  compound. 
As  a  result  of  testing  various  electrodes,  it  was  determined  that  silver/silver-chloride 
electrodes  from  In  Vivo  Metrics  Inc.  gave  the  best  performance.  We  use  ceramic  tubes 
filled  with  a  sodium-chloride  solution  to  connect  the  electrode  chambers  to  the  measure¬ 
ment  points  above  the  splitter  plate.  The  ends  of  these  tubes  must  be  colinear  and  equally 
spaced  for  a  given  measurement  axis.  The  tubes  are  removable  to  allow  replacement  in 
case  of  damage.  An  aluminum- oxide  ceramic  was  chosen  for  its  dimensional  stability, 
great  stiffness,  and  high  electrical  insulation  resistance.  The  ceramic  tubes  are  glued  to  a 
threaded  insert  using  a  jig  which  assures  that  all  the  tube  assemblies  will  be  interchange¬ 
able.  Dimensional  tolerances  are  such  that,  when  the  threaded  inserts  with  tube  are 
screwed  into  the  splitter  plate  from  the  side  opposite  the  electrodes,  the  tips  of  the 
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Figure  4.2.  Block  diagram  of  the  electronics  for  the  fifth  version  of  the  EMVM. 
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ceramic  tubes  are  aligned  and  equally  spaced  within  ±0.005  in.  To  reduce  the  Johnson 
noise  associated  with  the  resistance  of  the  salt-water  bridge  formed  within  the  ceramic 
tubes,  we  filled  them  with  a  water  solution  of  saturated  sodium  chloride  gelled  with 
agarose.  The  resistance  of  this  solution  is  about  three  times  lower  than  that  of  standard 
seawater. 

Agarose  is  a  polysaccharide  isolated  from  agar,  which  is  derived  from  marine  red 
algae.  We  purchase  this  product  as  a  dry  white  powder  from  FMC  BioProducts,  Rock¬ 
land,  Maine.  The  agarose  is  added  to  the  sodium-chloride  solution  in  the  amount  of  2% 
by  weight.  The  mixture  is  stirred,  starting  at  room  temperature,  and  slowly  raised  to 
100°C.  After  boiling  for  about  10  min,  the  solution  becomes  clear  and  will  gel  when 
cooled  to  room  temperature.  The  warm,  ungelled  solution  is  workable  for  many  minutes 
when  kept  on  a  warm  hot  plate.  The  saturated  solution  forms  a  salt-crystal  skin  on  the 
surface  if  exposed  to  air  owing  to  evaporation  of  water.  We  cover  the  mixing  container 
to  impede  evaporation  and  open  it  only  to  extract  the  ungelled  solution  as  needed. 

The  procedure  for  filling  the  electrode  chambers  and  ceramic  tubes  was  developed 
over  a  period  of  time  and  is  designed  to  eliminate  problems  associated  with  trapped  bub¬ 
bles.  After  mounting  the  electrode  chambers  to  the  backside  of  the  splitter  plate,  we  fill 
the  electrode  chambers  with  a  saturated  solution  of  sodium  chloride  (no  agarose).  Next  a 
partial  vacuum  is  drawn  from  the  front  side  of  the  splitter  plate.  This  draws  the  bubbles 
out  from  the  screw  threads  and  other  crevices  in  the  electrode  chamber.  After  the  bub¬ 
bles  are  removed,  the  ceramic  tubes  with  their  threaded  inserts  are  installed.  A  chamber 
filled  with  a  saturated  sodium-chloride  solution  is  placed  over  each  ceramic  tube  and, 
with  the  tip  of  the  tube  within  the  solution,  a  partial  vacuum  is  again  drawn.  At  this  point 
any  trapped  air,  either  in  the  threads  or  in  the  sodium-choride  solution,  is  drawn  out 
through  the  open  end  of  the  ceramic  tube.  The  heated,  ungelled  agarose  and  saturated 
sodium-chloride  solution  is  then  inserted  into  the  ceramic  tubes  with  a  hypodermic 
syringe.  This  must  be  done  carefully  so  that  additional  bubbles  are  not  created  inside  the 
tube.  The  ungelled  material  is  injected  while  the  syringe  needle  is  withdrawn  slowly  and, 
after  cooling,  forms  a  solid  plug  of  gelled  electrolyte  inside  the  ceramic  tube.  The  object 
of  this  procedure  is  to  produce  a  bubble-free  salt  bridge  within  the  tube.  While  not 
affecting  performance  in  the  test  tank,  small  bubbles  collapse  under  pressure  when  the 
instrument  is  operated  in  the  field.  The  repeated  pressure  cycling  works  the  gelled  plug 
loose  from  inside  the  ceramic  tubes  and  causes  the  sensor  to  fail. 

4.2.2  Splitter  plate 

What  we  refer  to  as  the  splitter  plate  is  a  flow-splitting  plate  20  in.  in  diameter  and 
3/4  in.  thick.  Within  the  outer  3  in.  of  the  plate’s  radius,  the  thickness  tapers  from  3/4  in. 
to  0  in.,  leaving  a  sharp  edge  at  the  outer  circumference.  The  taper  is  approximately  par¬ 
abolic  from  the  front  surface  to  the  back  surface,  which  remains  flat  put  to  the  edge.  The 
plate  is  machined  from  PVC  sheet  stock.  Fiberglass  studs  glued  into  blind  holes  on  the 
back  side  provide  the  means  for  securing  the  sensor  to  a  plastic  flange  attached  to  the  coil 
housing. 
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4.2.3  Magnets 

The  electromagnet  used  with  this  sensor  is  a  solenoid  consisting  of  1500  turns  of 
#16  AWG  copper  wire.  The  wire  was  wound  on  a  fiberglass  mandrel  chucked  in  a  lathe 
set  to  turn  very  slowly.  After  winding,  the  coil  assembly  was  potted  in  electronics  grade 
epoxy  potting  compound  selected  for  its  low  water  absorption.  The  epoxy  used  was 
Dexter  Hysol  RE2039/HD3561  mixed  in  the  ratio  of  100  parts  RE2039  to  30  parts 
HD3561  by  weight.  Because  the  curing  process  is  highly  exothermic,  castings  where  the 
weight  of  the  epoxy  is  greater  than  500  g  are  not  recommended.  The  encapsulation  of 
the  magnet  coil  required  three  separate  pours.  For  each  pour,  the  individual  components 
were  first  debubbled  using  a  vacuum,  and  after  mixing  a  final  vacuum  debubbling  was 
performed.  The  material  was  then  poured  into  the  mold  holding  the  coil  and  mandrel. 
The  outer  diameter  of  the  coil  is  11.5  in.,  the  inner  diameter  is  8.00  in.,  and  the  length  is 
4.00  in.  The  resistance  of  the  coil  is  15.8  Q,  and  the  inductance  is  0.597  H.  This  coil 
produces  a  magnetic  field  with  a  strength  of  9.5  G/A  at  a  distance  of  8.75  in.  from  the 
center  of  the  coil,  which  coincides  with  the  center  of  our  sensor  tube  array.  A  map  of  the 
magnetic  field  is  found  in  Figure  4.3.  The  coil  has  a  time  constant  (L/R)  of  38  ms,  which 
limits  the  upper  frequency  at  which  it  can  be  driven.  The  maximum  frequency  of  opera¬ 
tion  is  about  2  Hz  when  using  square-wave  excitation  and  requiring  that  the  field  remain 
constant  during  at  least  75%  of  the  cycle. 


Figure  4.3.  Map  of  the  magnetic  field  of  the  fifth  version  of  the  EMVM,  under  a  1-A  excitation 
current. 
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The  permanent  magnet  was  assembled  from  four  6  x  4  x  1  in.  ceramic  magnets  of 
commercial  manufacture.  The  magnetic  material  is  a  sintered  barium  ferrite  ceramic 
with  a  magnetic  remanence  of  2200  G.  These  magnets  are  magnetized  along  their  thin 
dimension,  yielding  pole  faces  that  are  6  x  4  in.  The  maximum  field  strength  measured  at 
the  pole  face  is  about  1 100  G.  This  is  about  half  the  textbook  value  for  the  remanence  of 
this  material.  The  edges  exhibited  the  strongest  fields,  and  the  center  of  the  pole  face 
dipped  to  about  900  G.  The  magnets  were  glued  together  with  epoxy  (Emmerson  & 
Cummings  Ecco  Bond  286)  which  has  a  very  high  (2500  psi)  shear  strength.  The  indivi¬ 
dual  magnets  were  held  together  with  wooden  clamps  while  the  epoxy  was  curing.  The 
final  result  was  a  magnet  6  x  8  x  2  in.  with  pole  faces  6  x  8  in.  This  magnet  assembly 
was  placed  in  a  mold  whose  outside  diameter  was  the  same  as  that  of  the  electromagnet. 
The  magnet  was  then  potted  in  an  epoxy  potting  compound  (Dexter  Hysol 
RE2039/HD3561,  100:30  by  weight)  to  protect  it  from  corrosion.  This  permanent  mag¬ 
net  assembly  produces  about  140  G  at  the  center  of  the  sensor  tube  array, 

4.2 A  Electronics,  optical  links,  batteries,  and  cabling 

The  electronics  are  now  separated  into  three  main  parts.  The  computer  interface 
that  was  used  with  the  previous  versions  of  the  instrument  remains  the  same.  What  had 
been  the  underwater  cable  on  the  previous  version  is  now  much  shorter  and  is  a  deck 
cable.  The  large  pressure  case  that  housed  the  underwater  electronics  in  most  of  previ¬ 
ous  versions  now  sits  on  the  deck.  The  low-noise  preamplifiers  and  their  associated  V/F 
converters  have  been  removed  to  their  own  small  pressure  case.  Also  inside  this  small 
case  are  a  +16-V  battery  pack  and  fiber-optic  cable  drivers  and  receivers.  The  fiber-optic 
cable  drivers  place  four  channels  of  FM  data  on  the  underwater  fiber-optic  cable:  One 
channel  for  vorticity,  two  for  orthogonal  velocity,  and  one  for  a  low-gain  dc-coupled 
average  of  the  outer  four  electrode  potentials  with  respect  to  the  center  electrode.  (The 
latter  is  a  diagnostic  channel  used  to  check  the  electrodes.)  A  single  fiber-optic  receiver 
provides  a  synchronization  signal  for  blanking  the  output  of  the  preamplifier  during  the 
switching  transient  when  the  electromagnet  is  used. 

The  low-noise  preamplifiers  are  the  single  most  important  part  of  the  instrument 
since  they  determine  the  noise  floor  below  which  the  data  are  no  longer  visible.  A 
schematic  diagram  of  the  preamplifier  for  a  single  axis  is  shown  in  Figure  4.4.  The 
configuration  is  a  standard  three  operational-amplifier  implementation.  The  input 
amplifiers  (AD  743s)  were  chosen  for  their  combination  of  low  voltage  noise  and  low 
current  noise.  Since  the  impedance  of  the  sensor  source  is  about  5(KX)  the  amplifier 
needs  this  combination  of  properties  to  be  effective.  A  comparison  of  the  noise  voltage 
densities  of  this  and  some  of  the  alternative  amplifiers  considered  is  shown  in  Figure  4.5. 
The  comparison  is  somewhat  unfair  to  the  Keithley  181  since  the  noise  shown  for  the 
other  operational  amplifiers  is  for  a  single  amplifier,  not  an  instrumentation 
configuration.  The  AD  743  has  the  lowest  total  noise  at  1  Hz  for  source  impedances 
greater  than  3(X)0  Q.  We  chose  Intersil  ID-KK)  dual  diodes  for  the  input  overload  protec¬ 
tion  circuitry  because  of  their  extremely  low  leakage  current  (1(X)  fA).  This  low  current 
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Figure  4.4.  Schematic  diagram  of  preamplifier  electronics  for  one  of  three  single-axis  chan¬ 
nels.  Analog  Devices  AD743  amplifiers  were  chosen  because  of  their  combined  low  voltage 
noise  and  low  current  noise  specifications.  The  input  protection  diodes  were  chosen  for  their 
extremely  low  leakage  currents.  The  potentiometers  are  adjusted  for  maximum  common-mode 
rejection  at  60  Hz. 
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Figure  4.5,  Total  noise,  referred  to  the  input,  for  three  candidate  EMVM  preamplifiers  as  a 
function  of  source  impedance.  The  source  impedance  of  the  EMVM  sensors  typically  ranges 
from  5x10^  to  lx  10^  Cl. 


flowing  through  the  source  impedance  generates  less  than  300  pV,  well  below  the  voltage 
noise  of  the  input  amplifiers.  Provision  is  made  for  adjusting  the  dc  common-mode  rejec¬ 
tion  with  a  balance  potentiometer  in  the  second  stage.  The  velocity-induced  signals  are 
equal  and  opposite  in  sign  at  the  input  to  amplifiers  Y1  and  Y2  (see  Figure  4.4),  and  we 
require  these  to  cancel  out  for  the  vorticity  calculation.  This  means  that  the  gains  of 
preamplifiers  Y1  and  Y2  must  be  the  same.  We  want  to  observe  vorticity-induced  signals 
that  are  as  small  as  1  x  10“"^  times  the  velocity-induced  signals.  The  gain  balance 
required  between  channels  is  therefore  1  part  in  10"^.  To  provide  the  required  gain 
matching  and  stability,  we  have  used  precision  Vishay  resistors  with  a  tolerance  of 
0.01%  and  a  temperature  coefiScient  of  2  ppm/°C.  We  have  chosen  a  relatively  low  value 
of  510  for  the  overall  preamplifier  gain  so  that  we  can  accommodate  the  large  voltages 
(10  mV)  produced  by  the  combination  of  electrode  offset  potentials  and  velocity-induced 
signals. 

The  amplified  signals  are  next  presented  to  track-and-hold  amplifiers.  These 
amplifiers  are  synchronized  with  the  coil  drive  of  the  electromagnet  and  hold  the 
preamplifier  output  starting  16.66  ms  before  the  switching  transient  and  lasting  for 
150  ms  thereafter.  The  amplifier  then  switches  back  and  tracks  the  preamplifier  output 
for  the  rest  of  the  time.  This  has  the  effect  of  blanking  the  switching  transient  from  the 
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next  stage,  which  for  vorticity  is  a  summing  amplifier  and  a  final,  ac-coupled  amplifier 
with  a  gain  of  400.  With  the  switching  transient  removed,  the  ac  amplifier  can  be  cou¬ 
pled  at  a  much  lower  frequency.  We  use  a  simple  one-pole  high-pass  filter  with  a  comer 
frequency  of  0.05  Hz.  For  velocity,  the  output  of  the  track-and-hold  amplifiers  is  ac  cou¬ 
pled  to  differencing  amplifiers  with  a  much  lower  gain.  The  gain  of  the  differencing 
amplifiers  is  only  10  so  that  the  large  signals  (1  mV  per  meter  per  second)  present  at  the 
electrodes  can  be  accommodated. 

The  amplified  vorticity  signals  and  the  two  orthogonal  velocity  channels  are  next 
applied  to  V/F  converters.  The  choice  of  converters  is  critical,  since  a  great  deal  of  our 
earlier  problems  can  be  traced  to  the  abnormal  but  subtle  behavior  of  the  previous  V/F 
converters.  Our  previous  problems  appeared  to  be  a  slight  cross-coupling  of  the  velocity 
and  vorticity  measurements.  It  turned  out  that  the  output  frequency  of  the  converters 
previously  used  (Burr  Brown  VFC  320s)  was  pulled  very  slightly  toward  the  operating 
frequency  of  the  adjacent  V/F  converter.  The  effect  was  small  and  subtle  and  proved 
difficult  to  correct  with  the  circuit-board  layout  already  in  place.  Because  space  is  at  a 
premium,  it  was  not  possible  to  separate  the  converters  enough  physically  to  eliminate 
the  problem  completely.  Our  new  circuits  use  VFC  KX)  integrated  circuits  also  from 
Burr  Brown.  The  VFC  KX)  uses  a  completely  different  conversion  technique.  The  con¬ 
verter  requires  an  external  clock  which  sets  the  full-scale  frequency  to  one-half  the  clock 
frequency.  In  addition,  output-frequency  pulse  transitions  are  forced  to  coincide  with 
clock  transitions.  This  introduces  a  considerable  amount  of  phase  jitter  in  the  output;  as  a 
result,  the  instantaneous  frequency  may  not  be  correct,  but  after  many  cycles  the  average 
frequency  is  correct.  We  utilize  a  single  system  clock  to  drive  all  four  V/F  converters  so 
that  all  four  converters  have  outputs  with  transitions  that  are  synchronous  with  the  same 
clock.  Under  these  conditions,  the  converters  are  not  able  to  influence  an  adjacent  chan¬ 
nel,  since  their  output  transitions  are  occurring  either  at  the  same  time  or  at  times  when 
the  adjacent  converter  is  not  sensitive.  Operating  at  a  clock  frequency  of  150  kHz,  the 
converters  have  a  center  frequency  of  37.5  kHz  and  a  deviation  of  7.5  kHz/V.  The 
linearity  is  excellent.  The  upper  panel  of  Figure  4.6  shows  the  straight-line  transfer  func¬ 
tion  of  the  converter  as  measured.  The  lower  panel  shows  the  deviation  from  a  straight- 
line  fit.  The  deviation  is  less  than  2  Hz  from  nominal  for  the  entire  70-kHz  range. 

The  output  of  the  V/F  converters  is  a  frequency-modulated  signal  in  the  5-75  kHz 
band.  Each  of  these  signals  is  applied  to  a  fiber-optic  cable  driver  which  modulates  a 
660-nm  LED  which  in  turn  illuminates  the  cable  fiber.  There  are  five  fiber-optic  cables, 
one  for  each  of  the  three  data  channels  and  one  for  the  synchronization/blanking  signal. 
Each  cable  consists  of  a  plastic  fiber  1  mm  in  diameter  and  a  plastic  jacket  with  an  out¬ 
side  diameter  of  2  mm.  The  fiber-optic  cable,  drivers,  and  receivers  are  manufactured  by 
the  Hewlett  Packard  Co.  and  sold  under  the  name  Versatile  Link.  These  components 
were  chosen  because  of  their  low  cost  and  their  suitability  for  our  low  data  rates  and  rela¬ 
tively  short  distances.  Our  previous  experience  with  the  instrument  determined  that  opti¬ 
cal  isolation  of  our  underwater  sensor  unit  and  the  rest  of  the  electronics  from  the  hull  of 
the  towing  ship  was  essential. 
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Figure  4.6.  Straight-line  transfer  function  of  V/F  converters  (top)  and  deviation  from  a 
straight-line  fit  (bottom).  The  upper  panel  shows  calibration  data  for  one  of  the  VFC  100 
voltage-to-frequency  converters  used  to  FM  encode  the  data  for  transmission  over  the  fiber-optic 
cable  link.  The  output  frequency  is  linearly  related  to  the  input  frequency.  The  fit  to  the  data  is 
represented  by  y,  which  is  given  at  the  upper  left.  The  bottom  panel  graphs  the  difference 
between  the  best-fit  straight  line  and  the  actual  data  as  a  function  of  input  voltage.  The  worst- 
case  deviation  is  2  Hz  out  of  75  kHz,  or  0.003%  of  full  scale. 
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The  four  fiber-optic  cables  are  bundled  together  with  the  two-conductor  neoprene- 
jacketed  cable  that  supplies  the  excitation  current  for  the  electromagnet  coil.  An  outer 
sheath  of  braided  plastic  fiber  provides  protection  against  abrasion  while  the  cable  is 
being  handled  on  deck.  The  overall  diameter  of  the  cable  bundle  and  protective  braid  is 
about  0.7  in.  In  addition,  a  Spectra  strength  member  with  an  outside  diameter  of 
0.375  in.  is  attached  to  the  cable  bundle.  This  synthetic  line  has  a  working  strength  of 
1000  lb,  more  than  sufficient  for  either  lifting  and  deploying  the  tripod  or  towing  the 
instrument  with  the  fiberglass  fairing. 

The  power  source  for  the  underwater  unit  is  a  bipolar  battery  pack  of  plus  and 
minus  15  V.  The  capacity  of  the  battery  pack  is  4.5  A-h.  The  batteries  are  of  the  gelled 
lead-acid  type.  These  were  chosen  for  their  reliability  and  relatively  high  energy  density. 
They  are  easily  recharged  overnight.  An  underwater  connector  provides  access  to  the 
battery  leads  for  charging  purposes.  A  special  jumper  plug  on  the  connector  is  used  to 
switch  on  the  power  to  the  internal  electronics. 

4.2.5  Pressure  vessels,  housings,  and  structures 

The  main  pressure  case  for  the  underwater  electronics  is  constructed  of  6061-T6 
aluminum.  It  has  an  outside  diameter  of  12.00  in.,  a  wall  thickness  of  0.625  in.,  and  an 
overall  length  of  15  in.  The  endcaps,  also  constructed  of  6061-T6,  are  1.90  in.  thick. 
They  seal  with  both  a  single  radial  0-ring  and  a  facial  0-ring.  One  endcap  has  a  special 
penetrator  which  provides  a  continuous  silver- wire  feedthrough  for  the  low-level  elec¬ 
trode  signals.  The  other  endcap  has  an  underwater  connector  which  supports  the  fiber¬ 
optic  cable  interface.  The  fiber-optic  cable  drivers  and  receivers  are  potted  in  an  alumi¬ 
num  tube  which  is  attached  to  the  back  of  a  24-pin  underwater  connector.  The  plastic 
optical  fibers  are  permanently  potted  into  the  driver/receiver  assembly,  which  operates  at 
hydrostatic  pressure.  The  aluminum  case  and  stainless  steel  fasteners  produce  a 
significant  amount  of  corrosion  with  associated  corrosion  currents.  These  currents  are 
not  significant  when  we  use  the  oscillating  magnetic  field  and  synchronous  demodula¬ 
tion.  They  are  troublesome,  however,  when  we  use  a  permanent  magnet  for  the  exciting 
field  and  low-frequency  ac-coupled  amplifiers  for  making  the  measurements.  This  prob¬ 
lem  was  alleviated  somewhat  by  coating  the  housing  with  vinyl  paint.  In  the  future,  this 
underwater  housing  will  be  changed  to  an  all  plastic  design,  perhaps  one  utilizing  glass- 
filled  Delrin. 

Several  tripods  have  been  constructed  to  allow  the  underwater  unit  to  operate  at  a 
fixed  height  off  the  bottom.  The  first  tripod  was  constructed  of  solid  steel  rod  0.75  in.  in 
diameter.  The  completed  tripod  with  instrument  attached  is  shown  in  Figure  4.7.  The 
structure  was  extremely  solid  and  relied  on  its  mass  to  keep  it  from  tipping  over.  During 
some  deployments  with  this  tripod,  unexplained  noise  was  present  which  we  felt  might 
be  due  to  corrosion  currents  associated  with  the  tripod  itself.  Consequently,  the  structure 
was  coated  with  vinyl  paint.  Later,  a  second  tripod  was  constructed,  this  time  entirely  of 
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Figure  4.7.  The  first,  all-steel  tripod  used  to  place  the  EMVM  near  the  bottom.  The  tripod  is 
shown  with  the  fifth  version  of  the  EMVM  mounted  so  that  measurements  of  the  vertical  com¬ 
ponent  of  vorticity  can  be  obtained. 


composite  glass  and  plastic  materials.  Joints  were  constructed  with  layed-up  fiberglass 
reinforcements  and  fiberglass  bolts.  Granite  curb  blocks,  each  weighing  100  lb  in  air, 
were  used  at  each  of  the  three  comers  to  provide  restoring  force  against  tip-over.  Fig¬ 
ure  4.8  shows  this  tripod  with  the  instrament  attached,  sitting  on  the  aft  deck  of  the  RA^ 
Miller. 

This  sensor  was  used  exclusively  with  the  GW  fairing  described  in  Section  3.5.2. 
This  fairing,  previously  shown  in  Figure  3.21,  allows  towing  the  sensor  in  a  vertical 
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upright  position  only.  The  size  and  shape  of  the  fairing  were  dictated  by  the  original 
underwater  sensor  and  electronics,  which  were  much  larger.  No  provision  was  made  to 
tow  the  sensor  horizontally,  since  it  was  envisioned  that  the  Seasoar  tow  body  would  sup¬ 
plant  this  fairing.  However,  this  fairing  remains  in  use  to  this  date  because  it  provides 
lower  vehicle  vibration  levels  at  low  (<2.5  kn)  tow  speeds. 

4.2.6  Shipboard  electronics 

The  EMVM  electronics  are  separated  into  three  main  assemblies.  The  first,  contain¬ 
ing  the  battery-isolated  preamplifiers  and  V/F  converters,  is  in  the  small  underwater 
housing  previously  described.  What  we  now  call  the  deck  electronics  was  previously 
underwater  with  the  preamplifiers.  Now  mounted  on  the  deck,  it  performs  the  frequency 


Figure  4.8.  All-plastic  tripod  with  fifth  version  of  the  EMVM  attached.  The  tripod  is  sitting  on 
the  aft  deck  of  the  research  vessel  WW  Miller.  The  structure  is  made  up  entirely  of  fiberglass. 
The  fasteners  are  glass-reinforced  plastic.  The  ballast  weights  (not  shown)  were  granite  blocks 
that  weighed  about  150  lb  in  water.  This  constraction  technique  was  chosen  to  eliminate  the 
generation  of  electric  fields  due  to  corrosion  of  the  tripod. 
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digitizing  and  digital  data  buffering  and  creates  the  asynchronous,  bit-serial, 
Manchester-encoded  data  stream  that  is  sent  over  a  twisted  pair  to  the  computer  inter¬ 
face.  It  also  provides  all  the  necessary  timing  signals  and  power  drivers  for  driving  the 
electromagnet.  Several  auxiliary  sensors  are  housed  with  these  electronics,  including  a 
triaxial  force-balance-accelerometer  set,  an  angular  rate  sensor,  and  a  compass.  A  block 
diagram  of  the  deck  electronics  is  shown  in  Figure  4.2. 

The  data  rate  is  adjustable  via  software  programming  of  the  timing  generator;  how¬ 
ever,  we  have  settled  on  using  a  basic  data  rate  of  two  samples  every  1.067  s.  One  of 
these  samples  is  taken  while  the  north  of  the  applied  magnetic  field  is  pointed  up;  the 
other  is  taken  when  the  field  is  reversed.  Since  an  up  and  a  down  measurement  are 
required  to  generate  an  estimate  of  vorticity  or  velocity,  this  yields  a  new  data  point 
every  1.067  s.  The  timing  generator  also  provides  a  variable  blanking  period  during 
which  the  FM  data  from  the  underwater  unit  are  ignored.  This  is  used  to  blank  out  the 
transients  that  occur  when  the  magnetic  field  switches  direction.  The  transients  are  pro¬ 
duced  in  whatever  residual  loop  area  is  left  in  the  sensor  wiring.  After  the  blanking  time, 
the  counters  used  to  digitize  the  FM  data  are  gated  on  for  an  exact  multiple  of  the  60-Hz 
line  frequency  period.  This  has  the  effect  of  integrating  the  signals  over  an  n-cycle  period 
of  the  local  60-Hz  field.  Signals  present  in  the  system  at  the  60-Hz  line  frequency  are 
greatly  rejected.  This  allows  us  to  test  the  instrument  in  the  laboratory  where  large 
60-Hz  electric  fields  are  present.  The  same  data  rate  applies  to  all  data  channels,  of 
which  there  are  17  plus  a  synchronization  word.  The  total  is  then  36  16-bit  words  every 
second.  These  data  are  applied  to  an  asynchronous,  serial  Manchester  adapter.  These 
bit-serial  data  are  clocked  out  over  a  twisted  pair  of  wires  at  100  kbit/s.  The  twisted-pair 
drivers  and  receivers  for  this  serial  link  are  transformer  coupled  to  the  line  at  each  end, 
thereby  reducing  common-mode  noise  and  grounding  problems. 

The  computer  interface  consists  of  two  parts:  the  serial  Manchester  decoder  with 
4096  word  FIFO  buffer,  and  a  16-bit  parallel  I/O  interface  with  proper  handshaking  logic 
for  an  HP-9020  computer.  The  4096- word  FIFO  provides  roughly  100  s  of  data- 
buffering  capacity  in  the  computer  interface.  This  may  seem  like  an  extreme  amount; 
however,  we  originally  anticipated  having  a  data  rate  15  times  faster.  The  HP-9020  runs 
UNIX,  which  is  not  designed  to  be  a  real-time  operating  system.  Consequently,  there  can 
be  periods  of  up  to  several  seconds  when  the  computer  is  unable  to  service  the  data- 
acquisition  task.  The  buffering  allows  the  computer  to  be  occupied  with  other  tasks  while 
the  data  are  saved  in  the  interface  unit. 

4.3  Performance  Evaluation 

4.3. 1  Sensitivity  and  self -noise  in  APL  test  tank 

Many  tests  have  been  conducted  in  which  the  sensor  was  operated  underwater  in  a 
salt-water  test  tank  measuring  6  ft  x  8  ft  x  4  ft  deep.  This  tank  is  wooden  framed  with  a 
fiberglass  liner.  It  is  electrically  nonconducting  and  is  isolated  from  ground  by  virtue  of 
resting  on  a  polyethylene  plastic  sheet  1/8  in.  thick  which  in  turn  sits  on  a  smooth,  flat 
concrete  floor.  The  tank  is  fitted  with  a  movable  carriage  which  can  translate  the  full 
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length  of  the  tank,  one  on  each  side  of  the  6-ft  dimension.  The  plenums  extend  to  the  full 
depth  of  the  tank.  Each  is  drilled  with  80  small  holes  in  a  regular  pattern  on  its  inner 
face.  The  interior  volume  of  each  plenum  houses  a  very  large  silver/silver-chloride  elec¬ 
trode  which  can  be  driven  with  up  to  20  mA  of  current.  This  arrangement  provides  a  uni¬ 
form  transverse  electric  field  of  known  magnitude. 

Figure  4.9  shows  the  noise  level  achieved  in  the  test  tank  with  the  instrument  sub¬ 
merged  to  the  back  of  the  splitter  plate.  No  signals  were  being  injected  via  the  plenums. 
The  water  was  extremely  still,  having  been  left  undisturbed  for  several  hours  after  the 
instrument  was  placed  into  the  tank  before  data  were  collected.  The  upper  panel  shows 
the  vorticity  time  series  for  120  s  of  data.  The  rms  variability  of  the  signal  over  the  120-s 
window  is  1.9  x  10"^  s“^ .  The  largest  part  of  the  signal  is  a  slow  drift  upward  after  30  s 
into  the  test.  The  lower  panel  shows  the  power  spectrum  of  vorticity  for  the  same  inter¬ 
val. 


Figure  4.9.  Fifth  version  EMVM  noise  levels  in  still  water.  The  upper  panel  shows  a  time 
series  of  vorticity  from  the  EMVM  while  it  was  operated  in  still  water  in  an  electrically  isolated 
all-plastic  test  tank.  The  lower  panel  shows  a  power  spectrum  of  the  data  from  the  upper  panel. 
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4.3.2  Summary  of  low  signals  in  Puget  Sound  experiments  ( upper  bound  on  noise) 

Several  tests  were  conducted  in  Blakely  Harbor,  a  small,  shallow  inlet  in  Bain- 
bridge  Island  on  Puget  Sound.  The  inlet  measures  0.4  n.mi.  across  at  its  mouth  and  tapers 
back  to  0.1  n.mi.  at  a  point  0.8  n.mi.  inside  the  harbor.  The  depth  is  about  70  ft  at  the 
mouth  and  shallows  up  to  6  ft  at  a  point  0.8  n.mi.  inside  the  harbor.  The  bottom  on  the 
inner  third  of  the  harbor  is  smooth  and  sandy.  Some  of  our  quietest  in  situ  data  have  been 
acquired  at  this  location.  Figure  4.10  shows  a  power  spectrum  of  vorticity  for  data  taken 
on  14  December  1990.  The  quietest  300  s  of  data  from  the  ebb  flow  portion  of  the 
deployment  were  used.  The  mean  velocity  of  water  flowing  past  the  tripod-mounted  sen¬ 
sor  during  this  time  was  0.08  m/s.  The  vorticity  signal  level  is  two  orders  of  magnitude 
greater  than  that  observed  in  the  still  water  of  the  test  tank. 


FREQUENCY.  Hz 

Figure  4.10.  Power  spectrum  of  vorticity  data  taken  in  a  quiet  location  in  Puget  Sound  on 
14  December  1990.  The  raw  data  were  filtered  to  remove  the  effects  of  intermittent  spikes  (simi¬ 
lar  to  those  shown  later  in  Figure  5.2).  The  raw  data  were  filtered  with  a  median  filter  whose 
length  was  10  points  of  the  original  raw-data  time  series.  The  filter  window  was  then  slid  by  five 
data  points,  and  the  10-point  median  again  calculated.  The  output  median-filtered  time  series  has 
a  data  rate  of  l/5th  the  original  sample  rate.  Since  the  original  sample  rate  was  0.9375  Hz,  the 
median-filtered  data  rate  is  0.1875  Hz.  This  makes  the  Nyquist  frequency  for  the  spectrum  of  Fig¬ 
ure  4.10  equal  9.375  x  10~^  Hz.  The  quietest  300-s  portion  of  the  ebb  flow  firom  a  tripod- 
mounted  deployment  in  Blakely  Harbor  was  used.  During  this  time,  the  mean  flow  past  the  sen¬ 
sor  was  0.08  m/s.  The  spectral  levels  are  two  orders  of  magnitude  greater  than  those  observed 
when  the  instrument  was  operated  in  the  still  water  of  a  test  tank. 
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5.  OBSERVATIONS  OF  VORTICITY  IN  PUGET  SOUND 

The  present  sensor  system  was  used  to  observe  background  levels  of  natural  vertical 
vorticity  at  various  sites  in  Puget  Sound  during  1990  and  1991.  (Earlier  EMVM  deploy¬ 
ments  were  used  to  estimate  instrument  performance.)  These  observations  were  designed 
to  study  the  strength  and  distribution  of  vorticity  in  strong  tidal  flow,  in  a  protected  har¬ 
bor,  and  in  the  wake  of  an  island  or  a  promontory. 

The  principal  measurement  sites  are  shown  in  Figure  5.1.  Point  Jefferson  is  at  the 
northeast  entrance  to  the  bight  of  Port  Madison.  Blakely  Rock  is  on  the  east  side  of 
Bainbridge  Island  near  the  entrance  to  Blakely  Harbor.  Restoration  Point  and  Decatur 
Reef  are  off  the  southeast  extremity  of  Bainbridge  Island.  Toliva  Shoal  is  in  the  middle 
of  Puget  Sound  off  Fox  Island  southwest  of  the  Tacoma  Narrows. 

5.1  Point  Jefferson 

On  13  December  1990,  the  EMVM  was  deployed  off  Point  Jefferson  on  a  steel  tri¬ 
pod  which  rested  on  a  sandy  bottom  about  0.3  n.mi.  offshore  at  a  depth  of  43  ft.  The 
instrument  was  put  on  the  bottom  just  after  slack  water  at  high  tide,  and  about  2  hours  of 
data  were  taken  during  the  ebb  flow.  The  flow  was  from  the  west,  and  the  rms  vorticity 
variability  was  quite  high  during  the  entire  deployment.  Figure  5.2  shows  a  representa¬ 
tive  time  series  for  the  13  December  1990  deployment  at  Point  Jefferson. 

The  data  from  this  deployment  (and  the  one  the  next  day  at  Blakely  Rock)  were 
contaminated  with  random  intermittent  spikes.  These  spikes  occurred  relatively  infre¬ 
quently  (one  to  two  per  minute)  and  mostly  when  the  instrument  was  cold.  It  was  later 
discovered  that  the  cause  of  these  spikes  was  an  intermittent  phase-lock  condition 
between  two  or  more  of  the  V/F  converters.  A  circuit  utilizing  newer  V/F  integrated  cir¬ 
cuits  was  designed,  constructed,  and  used  for  the  subsequent  deployments  off  Point 
Jefferson  during  mid-February  1991.  No  spiking  was  observed  at  any  time  during  these 
deployments. 

For  the  February  deployments,  the  EMVM  was  mounted  on  a  steel  tripod  which 
rested  on  the  bottom  of  Puget  Sound  at  a  depth  of  about  50  ft.  Data  were  collected  on 
three  days  for  about  6.5  hours  each  day. 

The  first  of  the  three  deployments  was  on  7  February  about  0.33  n.mi.  offshore  in 
water  45  ft  deep.  The  instrument  was  put  on  the  bottom  just  after  slack  water  at  high 
tide.  About  6.5  hours  of  data  were  taken  during  the  ebb  flow.  The  flow  was  from  the 
west,  and  the  rms  vorticity  variability  was  quite  high  during  the  entire  deployment. 

The  instrument  was  operated  with  the  oscillating  electromagnet  field,  and  the  sensor 
assembly  consisted  of  an  array  of  thin  ceramic  tubes  filled  with  a  gelled  agarose  and  con¬ 
centrated  sodium-chloride  solution.  These  were  the  smallest  diameter  tubes  used  to  date 
and  were  only  0.10  in.  in  outside  diameter. 
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Figure  5.1.  Map  of  Puget  Sound  with  EMVM  field  operations  marked  as  large  dots.  PJ  =  Point 
Jefferson;  BR  =  Blakely  Rock;  RP  =  Restoration  Point^ecatur  Reef  area;  TS  =  Toliva  Shoal. 
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Figure  5.2.  Time  series  of  data  from  a  deployment  off  Point  Jefferson  showing  random  inter¬ 
mittent  spikes.  The  spiking  was  caused  by  a  phase-lock  condition  between  two  or  more  of  the 
V/F  converters.  The  upper  panel  is  vertical  vorticity,  the  center  panel  is  current  direction  (CW 
from  north),  and  the  bottom  panel  is  current  speed.  Time  scale  is  in  two-digit  units  indicating 
from  left  to  right,  year,  month,  day  -  hour,  minutes,  seconds  (UTC). 


On  8  February,  the  instrument  was  deployed  at  nearly  the  same  spot  as  the  previous 
day  and  a  depth  of  38  ft.  The  instrument  was  deployed  during  slack  water  at  high  tide, 
and  6.5  hour  of  data  were  recorded  mainly  during  the  ebb  flow  from  the  west. 

For  this  deployment,  the  sensor  assembly  was  changed  to  one  having  no  ceramic 
tubes.  Instead,  the  measurement  ports  were  flush  with  the  large-diameter  plastic  plate 
separating  the  sensor  head  from  the  rest  of  the  instrument. 

The  third  deployment,  on  1 1  February,  was  at  the  same  location  off  Point  Jefferson 
in  44  ft  of  water.  During  this  deployment,  the  instrument  was  again  operated  with  the 
ceramic  tubes,  but  a  permanent  magnet  was  installed  in  place  of  the  electromagnetic  coil. 
The  strength  of  the  magnetic  field  generated  at  the  sensor  array  by  the  permanent  magnet 
was  about  three  times  the  peak  amplitude  produced  by  the  electromagnetic  coil.  About 
4  hours  of  data  were  collected,  starting  at  slack  water  at  low  tide  and  including  about 
3  hours  of  observations  during  the  flood  flow.  The  processing  software  was  changed  to 
acconunodate  the  lack  of  signal  chopping.  The  signals  were  ac  coupled  after  the  low- 
noise  preamplifiers,  providing  a  rolloff  below  3  x  10“^  Hz. 
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Figure  5.3  compares  the  power  spectra  for  the  different  magnet  configurations.  The 
solid  line  represents  the  spectrum  produced  from  2  hours  of  data  taken  when  the  elec¬ 
tromagnet  coil  was  energized.  The  dashed  line  represents  2  hours  of  data  taken  when  the 
instrument  was  operated  with  the  permanent  magnet.  Because  the  signal  is  not  synchro¬ 
nously  demodulated  with  the  instrument  in  the  permanent-magnet  configuration,  we 
expect  a  rolloff  of  the  signal  at  a  frequency  of  0.03  Hz.  For  this  comparison,  an  attempt 
was  made  to  use  data  taken  under  similar  environmental  conditions.  It  is  encouraging 
that  the  spectral  levels  are  comparable  over  a  wide  range  of  frequencies,  down  to  and 
including  the  roll-off  point.  Since  there  is  no  excess  noise  at  low  frequencies  in  the 
permanent-magnet  data,  we  may  be  able  to  lower  the  roll-off  frequency  in  the  electronics 
without  having  problems  with  excess  noise  due  to  electrode  drift.  Lowering  the  low- 
frequency  cut-off  point  in  the  electronics  should  be  relatively  easy  and  would  increase 
the  measurement  bandwidth. 


FREQUENCY,  Hz 

Figure  5.3.  Comparison  of  spectra  obtained  during  two  different  deployments  of  the  EMVM  on 
a  tripod  in  Puget  Sound.  A  different  electronic  configuration  was  used  for  each  deployment.  The 
solid  line  is  the  spectrum  of  vertical  vorticity  when  an  electromagnet  was  used.  TTie  dashed  line 
is  the  spectrum  when  a  permanent  magnet  and  an  electronic  rolloff  at  0.03  Hz  were  used. 
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Figure  5.4  compares  the  spectra  obtained  with  the  ceramic  straws  in  place  and  with 
the  straws  removed  and  replaced  with  small  orifices  mounted  flush  with  the  splitter  plate. 
Again,  an  attempt  was  made  to  use  data  taken  under  similar  environmental  conditions. 
The  flush-mounted  configuration  seems  to  have  excess  energy  at  frequencies  greater  than 
0.10  Hz.  This  may  be  due  to  fluctuations  in  the  boundary  layer  flow.  A  flush-mounted 
sensor  is  desirable  from  the  standpoint  of  mechanical  integrity,  but  if  the  excess  noise 
above  0.1  Hz  is  due  to  boundary  layer  fluctuations,  then  this  configuration  is  unaccept¬ 
able.  More  deployments  will  be  needed  to  determine  whether  the  spectra  continue  to 
exhibit  high  levels  at  0.01  Hz  and  higher. 

Figure  5.5  shows  two  spectra  taken  during  the  same  deployment;  the  lower  one  was 
with  the  electromagnet  turned  off  and  the  upper  one  with  it  turned  on.  There  is  a 
difference  of  nearly  two  orders  of  magnitude  (in  power)  between  the  two.  The  spectral 
level  of  the  lower  trace  is  somewhat  larger  than  had  been  observed  in  situ  during  previ¬ 
ous  deployments.  We  believe  this  is  because  the  data  used  for  this  spectrum  were  taken 
just  after  the  instrument  reached  the  bottom,  and  the  electrodes  were  still  coming  to 
equilibrium. 


FREQUENCY,  Hz 

Figure  5.4.  Comparison  of  spectra  obtained  from  the  EMVM  when  two  different  electrode  port 
configurations  were  used:  one  flush  mounted  (dashed  line)  and  one  at  the  ends  of  ceramic  straws 
(solid  line). 
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Figure  5.5.  Comparison  of  spectra  obtained  from  the  EMVM  with  the  electromagnet  turned  on 
and  with  it  turned  off. 


5.2  Blakely  Rock  and  Blakely  Harbor 

The  EMVM  was  deployed  6  ft  off  the  bottom  on  a  tripod  at  a  depth  of  46  ft  in 
Blakely  Harbor,  Bainbridge  Island,  on  14  December  1990.  The  purpose  of  this  deploy¬ 
ment  was  to  obtain  observations  in  waters  less  influenced  by  direct  tidal  flow  than  the 
measurements  at  Point  Jefferson  the  previous  day.  The  instrument  was  deployed  during 
the  flood  tide,  and  the  4-hour  record  spans  1.5  hours  before  slack  water  at  the  high  tide 
and  2.5  hours  after.  Both  flood  and  ebb  flows  were  observed,  with  velocities  up  to 
0.18  m/s  measured  during  the  beginning  of  our  observation  of  the  flood  flow.  Figure  5.6 
shows  the  entire  4-hour  record.  The  upper  panel  shows  the  vorticity  time  series.  The  first 
10  min  of  the  record  are  quiet  because  a  plastic  shroud  covered  the  sensor  and  protected 
it  from  motion  of  the  ambient  water.  The  mean  vorticity  during  the  shrouded  operation  is 
about  0.1  s~*.  This  is  due  to  an  offset  in  the  instrument  and  has  not  been  removed. 
Divers  removed  the  shroud  in  situ  and  calibrated  the  sensor  by  rotating  it  through  a 
known  angle  several  times.  The  rms  vorticity  was  1.6  x  10“^  s“'  during  the  flood  flow 
and  8  X  10~^  s“'  during  the  later  part  of  the  ebb  flow.  Thus  the  rms  vorticity  observed 
during  the  flood  flow  is  at  least  twice  that  of  the  ebb  flow.  During  the  flood  flow,  water  is 
advected  from  the  mouth  of  Blakely  Harbor  and  Puget  Sound  where  more  activity  is 
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Figure  5.6.  Data  taken  with  the  EMVM  deployed  on  a  tripod  at  a  depth  of  15  m  in  Blakely 
Harbor.  The  upper  panel  shows  the  vorticity  time  series  for  the  first  4  hours  of  deployment. 
There  was  a  flood  flow  during  the  first  2  hours  which  tapers  off  and  becomes  slack  water  at  the 
high  tide.  The  middle  panel  is  the  flow  direction  measured  CW  from  north.  The  bottom  panel 
shows  the  flow  speed  past  the  sensor.  The  flow  shows  a  dramatic  decrease  at  2120  UTC  near  the 
middle  of  the  record.  The  ebb  flow  slowly  builds  toward  the  later  part  of  the  record.  The  quietest 
part  of  the  record  occurs  during  the  ebb  flow,  when  water  that  was  sitting  in  the  inner  reach  of  the 
harbor  was  flowing  past  the  sensor. 


90  TR9503 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY. 


expected.  The  middle  panel  shows  the  direction  that  the  water  was  flowing  with  respect 
to  the  sensor  head.  The  direction  becomes  quite  variable  during  the  20  min  or  so  of  slack 
water  at  the  middle  of  the  record.  The  bottom  panel  shows  the  water  velocity  relative  to 
the  sensor.  The  highest  velocity  observed  is  0.2  m/s  at  the  beginning  of  the  deployment 
during  the  flood  flow.  Although  the  tides  cause  flow  into  and  out  of  the  harbor,  the  flow 
was  less  energetic  than  that  off  Point  Jefferson.  The  tripod  provides  the  most  stable  plat¬ 
form  for  the  instrument,  and  observations  from  it  are  thought  to  exhibit  the  lowest  noise 
level.  These  observations  are  compared  with  later  ones  obtained  by  towing  the  EMVM 
in  the  harbor. 

\ 

The  EMVM  was  towed  around  Blakely  Rock  and  Blakely  Harbor  for  several  days 
in  mid- April  1991.  The  instrument  was  housed  inside  a  fiberglass  fairing  with  the  sensor 
end  protruding  about  1  ft  below  the  tow  body.  The  instrument  was  towed  behind  the  R/V 
Miller  at  low  speeds  (2.5  kn)  and  a  depth  of  25  to  30  ft.  The  tows  were  primarily  in  a 
straight  line  from  within  Blakely  Harbor  to  a  point  within  several  hundred  yards  of 
Blakely  Rock  and  out  into  the  central  portion  of  Puget  Sound  adjacent  to  Elliot  Bay.  The 
first  tows,  on  16  April  1991,  were  to  the  south  of  Blakely  Rock,  as  we  were  observing 
during  a  flood  flow  (i.e.,  south-flowing  current).  The  tows  on  17  April  were  to  the  north 
of  the  rock,  as  we  were  observing  during  an  ebb  flow.  During  this  series  of  deployments, 
a  vehicle  dynamics  package  intended  to  evaluate  the  Seasoar  tow  body  was  attached  to 
the  GW  fairing  for  the  EMVM.  In  addition,  an  Eagle  8  GPS  receiver  was  installed  on 
board  R/V  Miller,  and  position  data  were  logged  at  a  1-Hz  rate. 

For  the  tows  on  16  and  17  April,  the  EMVM  was  operated  with  the  electromagnet 
and  with  ceramic  sensor  tubes  0.10  in.  in  diameter  extending  2.5  in.  beyond  the  splitter 
plate.  These  tubes  were  1  in.  shorter  than  in  previous  deployments.  This  change  was 
made  because  the  longer  tubes  were  extremely  fragile,  often  breaking  after  only  15  min 
of  towing.  Because  the  measurements  made  in  February  with  the  flush-mounted  sensor 
ports  produced  higher  noise  levels,  we  decided  to  try  an  intermediate-length  tube  to  gain 
durability  without  sacrificing  noise  performance.  An  in  situ  calibration  was  performed  at 
the  end  of  the  second  day  of  operations,  in  which  we  rotated  the  instrument  while  it  was 
suspended  from  the  overhead  hoist  off  the  aft  end  of  R/V  Miller  while  the  vessel  drifted 
slowly  in  Blakely  Harbor.  Even  though  we  were  not  under  way,  there  was  an  increase  in 
the  vorticity  noise  level,  presumably  because  of  interference  from  electric  fields  associ¬ 
ated  with  the  zinc  anodes  and  the  bronze  propeller. 

Figure  5.7  shows  a  time  series  of  the  vorticity  data  for  the  first  4  hours  of  towing 
near  Blakely  Rock  on  16  April.  The  lower  panel  shows  the  vorticity  corrected  for  yaw 
rate  and  bandpass  filtered  as  a  function  of  time.  The  upper  panel  shows  the  rms  vorticity 
variability  computed  over  20-s  windows.  During  the  least  active  portions  of  the  tow,  the 
rms  variability  was  slightly  less  than  5  x  10“^  s“^ .  This  is  slightly  lower  than  what  was 
achieved  with  the  sensor  mounted  on  a  tripod  moored  at  the  bottom  of  Blakely  Harbor. 
The  quietest  rms  vorticity  levels  observed  on  the  tripod  during  the  14  December  1990 
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Figure  5.7.  Data  from  the  EMVM  towed  with  the  GW  fairing  on  16  April  1991  in  the  wake  of 
Blakely  Rock.  The  lower  panel  shows  vorticity  corrected  for  yaw  rate  and  bandpass  filtered  as  a 
function  of  time.  The  upper  panel  shows  the  rms  vorticity  computed  over  20-s  windows.  There 
are  several  instances  where  the  rms  vorticity  was  >  0.04  s~*  while  traversing  the  wake. 

deployment  were  8  x  10"^  s"* .  Figure  5.7  shows  six  separate  events  when  the  rms  varia¬ 
bility  is  greater  than  4  x  10“^  s"^ .  The  most  active  event  has  an  rms  variability  of 
4.5  X  10“^  s“' .  These  events  occurred  when  the  sensor  passed  just  south  of  Blakely 
Rock.  The  wake  of  the  rocks  was  to  the  south,  as  we  were  towing  during  flood  flow  of 
the  tidal  cycle.  Figure  5.8  shows  the  tide  heights  for  Restoration  Point  on  16-18  April 
1991.  The  heights  at  the  times  that  we  were  towing  have  been  highlighted. 

Figures  5.9  through  5.16  show  vorticity  variance  computed  for  20-s  intervals  from 
data  taken  on  16  April.  The  vorticity  variance  is  plotted  as  a  variable-sized  dot  whose 
center  corresponds  to  the  latitude  and  longitude  of  the  average  position  for  that  particular 
interval.  The  largest  dots  correspond  to  a  vorticity  variance  of  3  x  10“^  s“^,  the  smallest 
dots  to  a  variance  <2  x  10~^  s~^.  It  can  be  seen  from  these  plots  that,  in  addition  to  the 
large  spatial  variability  associated  with  Blakely  Rock,  there  is  a  temporal  variability 
which  is  not  entirely  explained  by  the  variation  in  mean  flow  due  to  tidal  forcing.  In 
addition  to  the  area  just  south  of  the  rock,  there  is  another  area  with  elevated  vorticity 
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Figure  5.8.  Tide  heights  at  Restoration  Point  for  16-18  April  1991.  Data  for  the  times  of  the 
tows  are  highlighted. 


levels  between  Blakely  Rock  and  the  entrance  to  Blakely  Harbor  (see  Figures  5.9  and 
5.16).  The  source  may  be  Tyee  Shoal  2  km  to  the  north  or  the  point  just  to  the  north  of 
the  entrance  to  Blakely  Harbor. 

Figure  5.17  shows  the  ensemble  average  of  five  vertical  vorticity  spectra  taken 
when  the  sensor  was  passing  through  the  wake  of  Blakely  Rock.  Also  shown  in  dashed 
lines  is  the  ensemble  average  of  three  spectra  taken  when  the  sensor  was  in  deep  water  in 
the  middle  of  Puget  Sound.  The  spectra  have  similar  shapes,  but  the  spectrum  for  the  lee 
of  the  rock  is  elevated  by  a  factor  of  3.  The  rolloff  at  small  wavenumber  (high  fre¬ 
quency)  is  due  to  low-pass  filtering  during  the  processing,  which  runs  a  four-point 
weighted  window  filter  through  the  time-series  data. 

Figure  5.18  shows  similar  time  series  data  from  the  second  day  of  towing.  Again, 
there  are  several  periods  with  elevated  vorticity.  These  occurred  during  passage  of  the 
instrument  through  the  lee  of  Blakely  Rock.  During  this  time  series,  the  wake  of  the  rock 
was  to  the  north,  as  we  were  towing  during  the  ebb  flow  of  the  tidal  cycle.  During  the 
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Figure  5.9.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  2120-2143  UTC  on  16  April  1991.  The  dots  represent  variance  between  3  x  10“^  s“^  and 
2x  10"^  s"^ 


Figure  5.10.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  2146-2220  UTC  on  16  April  1991. 
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Figure  5,11.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  2223-2256  UTC  on  16  April  1991. 


Figure  5.12.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  2258-2328  UTC  on  16  April  1991. 
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Figure  5.13.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  2331-2352  UTC  on  16  April  1991. 


Figure  5.14.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  2355-0017  UTC  on  16-17  April  1991. 
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Figure  5.15.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  0020-0042  UTC  on  17  April  1991. 


Figure  5.16.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  0044-01 1 1  UTC  on  17  April  1991. 
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Figure  5.17.  Comparison  of  vertical  vorticity  spectra  from  active  and  quiet  regions  of  Puget 
Sound  near  Blakely  Rock.  Solid  line:  average  of  five  spectra  taken  while  towing  the  sensor 
through  the  wake  of  Blakely  Rock.  Dashed  line:  average  of  three  spectra  taken  while  towing  the 
sensor  in  deep  water  in  the  middle  of  Puget  Sound. 


TIME,  s 

Figure  5.18.  Time  series  of  vertical  vorticity  on  a  tow  with  the  GW  fairing  on  17  April  1991. 
The  lower  panel  shows  vorticity  corrected  for  yaw  rate  and  bandpass  filtered  as  a  function  of 
time.  The  upper  panel  shows  the  variability  in  rms  vorticity  computed  over  20-s  windows. 
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later  part  of  the  tow,  the  flow  became  weak,  and  slack  water  occurred  at  about 
2040  UTC.  This  corresponds  with  the  quietest  part  of  the  record,  which  occurs  at 
2030  UTC.  Figures  5.19  through  5.24  show  vorticity  variance  for  20-s  intervals  of  data 
from  17  April.  Again,  the  size  of  the  dots  corresponds  to  the  vorticity  variance,  and  the 
average  latitude  and  longitude  of  the  20-s  interval  determines  the  center  of  the  dot.  The 
highest  vorticity  levels  were  observed  north  and  considerably  east  of  Blakely  Rock.  It  is 
difficult  to  attribute  this  active  area  to  Blakely  Rock  because  the  area  is  so  far  to  the  east. 
The  currents  for  this  area  are  not  well  determined,  however,  so  we  cannot  rule  out  the 
possibility  of  Blakely  Rock  as  the  source.  There  are  also  periodic  elevated  vorticity  lev¬ 
els  in  the  area  between  Blakely  Rock  and  the  mouth  of  Blakely  Harbor.  Some  of  this 
activity  is  apparently  due  to  the  promontory  just  to  the  north  of  the  harbor  entrance. 

Figure  5.25  shows  two  spectra  of  vorticity  data  taken  in  Blakely  Harbor.  The  solid 
line  is  the  spectrum  computed  from  data  taken  on  14  December  1990  when  the  sensor 
was  mounted  on  a  tripod  resting  on  the  bottom.  The  dashed  line  is  the  spectrum  com¬ 
puted  from  data  taken  on  16  April  1991  when  the  sensor  was  being  towed  at  a  depth  of 
30  ft  and  a  speed  of  approximately  2.5  kn.  The  spectra  have  been  computed  with  respect 
to  wavenumber,  so  effects  of  the  speed  of  the  sensor  are  accounted  for. 

In  conclusion,  the  tests  of  the  EMVM  found  clear  evidence  of  environmental  signals 
with  rms  amplitudes  as  great  as  4.5  x  10“^  s~^  Blakely  Rock  seems  to  be  the  source  for 
most  but  not  all  of  the  high-vorticity  areas.  The  sporadic  high  levels  observed  away  from 
the  rock  are  probably  due  to  distant  promontories  and  shoals,  vessel  traffic,  and  biologi¬ 
cal  activity.  The  data  for  a  given  area  show  a  great  deal  of  temporal  variability,  only 
some  of  which  is  correlated  with  the  tidal  cycle. 

The  EMVM  was  deployed  again  in  Puget  Sound  on  two  consecutive  days  in  mid- 
June  1991.  The  instrument  was  housed  inside  the  GW  fairing  with  the  sensor  end  pro¬ 
truding  about  1  ft  below  the  tow  body  and  was  towed  behind  R/W  Miller  at  low  speeds 
(2.5  kn)  and  a  depth  of  about  30  ft.  The  tows  on  both  days  were  to  the  south  of  Blakely 
Rock,  as  we  were  observing  during  a  flood  flow.  During  this  series  of  deployments,  the 
vehicle  dynamics  package  for  Seasoar  was  attached  to  the  GW  fairing.  In  addition,  an 
Eagle  8  GPS  receiver  was  installed  on  board  RfV  Miller,  and  position  data  were  logged 
at  a  1-Hz  rate. 

The  EMVM  was  operated  with  a  permanent  magnet  and  ceramic  sensor  tubes 
0.10  in.  in  diameter  extending  2.50  in.  beyond  the  flow-splitting  plate.  An  in  situ  calibra¬ 
tion  was  performed  at  the  beginning  and  end  of  each  day  of  operations.  The  calibration 
consisted  of  rotating  the  instrument  while  it  was  suspended  from  the  overhead  hoist  off 
the  aft  end  of  RA^  Miller  as  the  vessel  was  drifting  slowly  in  Blakely  Harbor. 

Figure  5.26  shows  a  time  series  of  vorticity  data  from  the  4  hours  of  towing  near 
Blakely  Rock  on  12  June.  The  lower  panel  shows  vorticity  corrected  for  yaw  rate  and 
bandpass  filtered  as  a  function  of  time.  The  upper  panel  shows  the  rms  vorticity  variabil¬ 
ity  computed  over  20-s  windows.  Although  we  use  the  term  rms  variability,  we  are  actu¬ 
ally  calculating  the  standard  deviation  of  the  data  about  the  mean  over  an  interval  of  k 
points.  The  calculation  we  perform  to  yield  the  "rms"  vorticity  is 
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Figure  5.19.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  1728-1748  UTC  on  17  April  1991.  The  size  of  the  dots  corresponds  to  vorticity  variance, 
and  the  average  latitude  and  longitude  of  the  20-s  interval  determines  the  center  of  the  dot. 
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Figure  5.20.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  1750-1806  UTC  on  17  April  1991. 
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Figure  5.21.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  1808-1846  UTC  on  17  April  1991. 


Figure  5.22.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  1848-1932  UTC  on  17  April  1991. 
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Figure  5.23.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  1935-2014  UTC  on  17  April  1991. 
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Figure  5.24.  Vorticity  variance  along  ship’s  track  computed  at  20-s  intervals  for  data  collected 
from  2017-2054  UTC  on  17  April  1991. 
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Figure  5.25.  Two  spectra  of  vertical  vorticity  data  collected  in  Blakely  Harbor.  The  solid  line 
is  the  spectrum  computed  from  data  taken  when  the  sensor  was  mounted  on  a  tripod  resting  on 
the  bottom.  The  dashed  line  is  the  spectrum  computed  from  data  collected  while  the  sensor  was 
towed  at  a  depth  of  30  ft  at  approximately  2.5  kn. 
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Figure  5.26.  Time  series  data  from  a  tow  on  12  June  1991 .  The  lower  panel  shows  vertical  vor¬ 
ticity  corrected  for  yaw  rate  and  bandpass  filtered  as  a  function  of  time.  The  upper  panel  shows 
the  variability  in  rms  vorticity  computed  over  20-s  windows. 
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rms  vorticity  = 


^  2  ^ 

X  xl-kx^ 

n=\ 


k-\ 


(5.1) 


where  x  and  n  are  the  individual  vorticity  measurements  and  x  is  the  mean  of  the  vorticity 
measurements  over  the  interval  Qc  =  20  for  a  20-s  interval). 

During  the  least  active  portions  of  the  tow,  the  rms  variability  was  slightly  less  than 
5  X  10“^  s“^ .  This  is  slightly  lower  than  what  was  achieved  with  the  sensor  mounted  on  a 
tripod  moored  at  the  bottom  of  Blakely  Harbor  on  14  December  1990.  The  quietest  vor¬ 
ticity  levels  observed  on  the  tripod  during  the  14  December  deployment  were 
8  X  10“^  s“*  rms.  Figure  5.26  shows  nine  separate  events  where  there  is  a  signal  greater 
than  6  x  10“^  s"^  rms.  The  most  active  area  has  an  rms  variability  of  8.3  x  10“^  s“^ . 
These  events  occurred  when  the  sensor  passed  just  south  of  Blakely  Rock.  The  wake  of 
the  rock  was  to  the  south,  as  we  were  towing  during  flood  flow  of  the  tidal  cycle.  Figure 
5.27  shows  the  tide  heights  and  average  current  for  Restoration  Point  on  12-14  June 
1991.  The  tow  times  have  been  highlighted.  Figure  5.28  shows  the  vorticity  variance 
computed  for  20-s  intervals  for  data  taken  on  12  June.  The  vorticity  variance  is  plotted 
as  a  variable-sized  dot  whose  center  corresponds  to  the  latitude  and  longitude  of  the 
average  position  for  the  particular  interval.  The  legend  in  the  lower  left  shows  the 
correspondence  of  dot  size  to  vorticity  variance.  In  this  case,  vorticity  variance  is  just 
the  square  of  the  rms  vorticity  signal  levels  computed  over  the  20-s  window.  We  prefer 
plotting  vorticity  variance  because  it  is  more  indicative  of  the  power  contained  in  the 
signal. 
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Figure  5.27.  Tide  heights  and  average  current  for  Restoration  Point  on  12-14  June  1991.  Tow 
times  have  been  highlighted. 
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Figure  5.28.  Vertical  vorticity  variance  computed  at  20-s  intervals  for  data  collected  on 
12  June  1991.  The  variance  is  plotted  as  a  variable-sized  dot  whose  center  corresponds  to  the 
latitude  and  longitude  of  the  average  position  for  a  particular  interval.  The  legend  in  the  lower 
left  shows  the  correspondence  of  dot  size  to  vorticity  variance. 


Figure  5.29  shows  a  similar  time  series  of  data  from  the  second  day  of  towing. 
Again,  there  are  several  periods  with  elevated  vorticity  levels.  The  most  active  has  a  sig¬ 
nal  level  of  1.1  X  10“*  s“^  rms.  Figures  5.30  and  5.31  show  the  vorticity  variance  (square 
of  the  rms  vorticity  over  a  20-s  interval)  for  data  from  13  June.  Figure  5.30  is  for  the  first 
part  of  the  tow  when  the  flow  was  weak.  (The  surface  flow  seemed  to  increase  as  we 
towed  farther  south  and  around  the  mid-channel  marker  buoy.  We  could  only  estimate 
the  flow  based  on  the  set  taken  by  our  vessel  as  we  tried  to  steer  an  east/or/west  course.) 
Figure  5.31  shows  data  taken  later  in  the  day  when  the  flood  flow  was  the  greatest.  Both 
Figures  5.30  and  5.31  show  a  shadow  zone  directly  behind  Blakely  Rock  where  the  vorti¬ 
city  level  is  low.  The  highest  vorticity  levels  occur  in  two  bands  just  south  of  Blakely 
Rock  during  the  peak  of  the  flood  flow.  One  band  is  slightly  to  the  west  of  the  rock  and 
the  other  is  slightly  to  the  east,  areas  where  the  water  has  flowed  past  the  edge  of  the 
rock. 

Figure  5.32  shows  a  blowup  of  one  particular  event  in  relation  to  signals  observed 
during  the  entire  day  of  towing.  The  upper  panel  shows  rms  vorticity  levels  for  the  entire 
deployment  on  13  June.  An  arrow  points  to  the  signal  that  is  expanded  in  time  in  the 
lower  panel.  This  time  series  was  taken  as  the  EMVM  passed  through  an  active  area  just 
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Figure  5.29.  Time  series  data  from  a  tow  on  13  June  1991.  The  lower  panel  shows  vertical  vor- 
ticity  corrected  for  yaw  rate  and  bandpass  filtered  as  a  function  of  time.  The  upper  panel  shows 
the  variability  in  rms  vorticity  computed  over  20-s  windows. 
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Figure  5.30.  Vertical  vorticity  variance  for  data  from  13  June  1991  taken  during  a  time  when 
the  tidal  flow  was  weak. 
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Figure  5.31.  Vertical  vorticity  variance  for  data  from  13  June  1991  taken  during  the  time  of 
maximum  flood  tide. 


to  the  southeast  of  Restoration  Point.  This  area  corresponds  to  the  large  dot  at  the  begin¬ 
ning  of  the  track  near  Decatur  Reef  on  Figure  5.31.  The  tow  speed  at  this  point  was 

1.2  m/s  and  the  event  spans  a  time  period  of  60  s.  The  physical  extent  of  this  feature  is 
therefore  about  70  m  along  the  direction  of  the  tow,  which  was  perpendicular  to  the  flow 
from  the  north. 

In  conclusion,  the  EMVM  tests  in  June  1991  found  clear  evidence  of  environmental 
signals  with  rms  amplitudes  as  great  as  1.1  x  10“*  s~* .  Blakely  Rock  and  the  point  at  the 
north  of  the  entrance  to  Blakely  Harbor  seem  to  be  the  source  for  most  of  the  vorticity 
observed  in  the  area  north  of  Restoration  Point.  Restoration  Point  also  seems  to  be  a 
source  for  vorticity  observed  in  the  area  to  its  sOuth.  The  data  for  a  given  area  show  a 
great  deal  of  temporal  variability,  only  some  of  which  is  correlated  with  the  tidal  cycle. 

5.3  Restoration  Point  and  Decatur  Reef 

In  July  1991,  the  EMVM  was  deployed  again  for  2  days  in  Puget  Sound.  The 
instrument,  equipped  with  the  GW  fairing,  was  towed  behind  R/V  Miller  at  low  speeds 
(2.5  kn)  and  a  depth  of  about  30  ft.  On  both  days,  the  tows  were  to  the  south  of  Restora¬ 
tion  Point,  just  south  of  Blakely  Rock,  during  a  flood  flow.  During  these  tows,  a  vehicle 
dynamics  package  was  attached  to  the  instrument.  In  addition,  an  Eagle  4  GPS  receiver 
was  installed  on  board  R/V  Miller,  and  position  data  were  logged  at  a  1-Hz  rate. 
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Figure  5.32.  Blowup  of  one  particular  event  in  relation  to  signals  observed  during  the  entire 
day  of  towing.  Upper  panel  shows  the  rms  vertical  vorticity  signal  levels  for  the  entire  deploy¬ 
ment  on  13  June  1991.  Arrow  points  to  the  data  segment  that  is  expanded  in  time  in  the  lower 
panel. 

For  these  deployments,  the  EMVM  was  operated  with  a  permanent  magnet  which 
provided  an  average  field  strength  of  1.4  x  10"^  T  (140  G)  within  the  sensing  volume. 
This  is  approximately  three  times  the  field  strength  of  the  previous  electromagnet 
configuration.  As  in  the  June  tests  at  Blakely  Rock,  the  instrument  was  fitted  with 
ceramic  sensor  tubes  0.10  in.  in  diameter  extending  2.50  in.  beyond  the  flow-splitting 
plate.  The  small-diameter  tubes  seemed  to  produce  a  lower  level  of  self-induced  noise, 
and  the  short  length  provided  greater  durability.  The  towing  bridle  was  adjusted  from  the 
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June  deployments  to  provide  a  pitched-up  condition  of  about  7°,  which  seemed  to  give  a 
lower  background  noise  over  a  wider  range  of  operating  conditions,  including  ship  heave 
and  turning  comers.  An  in  situ  calibration  was  performed  at  the  beginning  of  each  day  of 
operations  by  rotating  the  instrument  while  it  was  suspended  from  the  overhead  hoist  off 
the  aft  end  of  the  RA^  Miller  as  the  vessel  was  slowly  drifting. 

Figure  5.33  shows  a  time  series  of  the  vorticity  data  for  10  July.  The  lower  panel 
shows  the  vertical  vorticity  corrected  for  yaw  rate  and  bandpass  filtered  as  a  function  of 
time.  The  upper  panel  shows  the  rms  vorticity  variability  computed  over  20-s  windows. 
During  the  least  active  portions  of  the  tow,  the  rms  vorticity  variability  was  slightly  less 
than  4.5  x  10“^  s“* .  In  the  most  active,  it  was  8  x  10“^  s~^ .  Figure  5.34  shows  the  tide 
heights  and  average  current  for  Restoration  Point  on  10-12  July  1991.  The  tow  times 
have  been  highlighted. 

For  this  series  of  measurements,  two  6-in.-diameter  by  12-in.-long  buoys  with  4  ft 
by  8  ft  "window  shade"  drogues  suspended  by  a  thin  line  to  a  depth  of  25  ft  were 
launched  in  areas  of  high  activity  near  Restoration  Point.  To  aid  visibility,  a  small  black 
flag  was  attached  to  the  top  of  a  3-ft-high  pole  running  through  the  center  of  each  buoy. 


910710-1800  910710-1900  910710-2000  910710-2100  910710-2200  910710-2300  910711-0000 

time;  s 


Figure  5.33.  Time  series  of  vertical  vorticity  data  for  10  July  1991.  Lower  panel:  Vorticity 
corrected  for  yaw  rate  and  bandpass  filtered  as  a  function  of  time.  Upper  panel:  variability  in 
rms  vorticity  computed  over  20-s  windows. 
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HOURS  SINCE  10  JULY  1991.  PST 

Figure  5.34.  Tide  heights  and  average  currents  for  Restoration  Point  during  10-12  July  1991. 
Tow  times  have  been  highlighted. 

The  ship’s  track  and  the  tracks  of  the  two  drogue  buoys  are  shown  in  Figure  5.35.  The 
positions  of  the  boat  are  plotted  every  minute  with  an  additional  tick  mark  and  time  label 
every  tenth  minute.  Four  positions  of  the  drogue  buoys  are  also  plotted.  These  four  posi¬ 
tions  were  calculated  from  a  series  of  range  and  bearing  measurements  taken  with  the 
Miller’s  radar.  The  motion  of  the  buoys  should  correspond  closely  to  that  of  the  water  at 
a  depth  of  25  ft. 

The  average  velocity  for  the  two  buoys  was  0.84  kn.  The  average  velocity  between 
the  first  and  second  positions  measured  for  buoy  #2  was  1.17  kn,  the  highest  observed. 
We  placed  buoy  #1  in  an  area  of  moderate  vorticity  activity  during  our  first  westbound 
towing  leg.  It  turned  out  that,  as  we  progressed  more  to  the  west  and  closer  to  Restora¬ 
tion  Point,  the  vertical  vorticity  signal  levels  increased.  After  towing  through  the  area 
with  the  peak  vorticity  signal,  we  turned  back  to  the  east  and  deployed  buoy  #2  in  a  par¬ 
cel  of  water  that  had  much  higher  signal  levels. 

Figure  5.36  shows  data  from  the  first  part  of  the  day,  when  we  tagged  a  piece  of 
water  with  the  drogue  buoys  and  followed  them  as  best  we  could  by  towing  a  serpentine 
course  which  brought  the  instrument  alongside  the  buoys  periodically.  The  highest  verti¬ 
cal  vorticity  levels  were  observed  in  an  area  southeast  of  Restoration  Point.  This  is  the 
expected  location  for  large  vorticity  signals  given  the  velocities  observed  with  the  drogue 
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Figure  5.35.  Ship’s  track  (dots)  and  drogue  buoy  tracks  (lines)  in  areas  of  high  vertical  vorti- 
city  activity  near  Restoration  Point.  The  positions  of  RA^  Miller  are  plotted  every  minute,  with 
an  additional  tick  mark  and  time  label  every  tenth  minute. 


buoys.  If  we  use  the  initial  velocity  vector  of  buoy  #2  to  extrapolate  back  to  the  source, 
we  find  that  the  parcel  of  water  passed  along  the  edge  of  Restoration  Point  and  over 
Decatur  Reef.  It  is  interesting  to  note  that  we  were  able  to  observe  the  signature  of  the 
Decatur  Reef  source  on  our  southernmost  transect.  The  peak  vorticity  signal  levels  on 
this  transect  correlate  well  with  our  point  of  closest  approach  to  buoy  #2.  At  this  time 
(2050  UTC),  the  parcel  of  water  had  been  traveling  for  about  2  hours  since  its  passage 
over  Decatur  Reef.  The  maximum  vertical  vorticity  signal  level  decreased  only  slightly 
during  this  interval.  We  observed  a  maximum  vorticity  signal  level  of  6.4  x  10“^  s~^  rms 
during  our  first  transect  near  Restoration  Point  and  a  maximum  signal  level  of 
6.1  X  10“^  s“^  rms  during  our  southernmost  transect  80  min  later.  What  may  be  more 
significant  is  that  the  size  of  the  area  with  high  vertical  vorticity  signal  levels  seems  to 
have  decreased.  During  the  passes  near  Restoration  Point,  the  area  with  high  vorticity 
signal  levels  was  about  300  m  across  in  the  direction  of  the  tow  (east-west).  During  the 
southernmost  transect  80  min  later,  the  area  with  high  vorticity  signals  was  about  150  m 
across. 
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Figure  5.36.  Vertical  vorticity  data  from  the  first  half  of  10  July  1991.  The  variance  is  plotted 
as  a  variable-sized  dot  whose  center  corresponds  to  the  latitude  and  longitude  of  the  average 
position  for  a  particular  interval.  The  legend  in  the  lower  left  shows  the  correspondence  of  dot 
size  to  vorticity  variance. 


During  the  second  half  of  the  day,  we  towed  the  EMVM  through  a  pattern  designed 
to  map  the  vertical  vorticity  levels  in  the  area  just  to  the  south  of  Restoration  Point.  For 
these  transects,  we  decided  to  proceed  farther  to  the  west  before  turning  back  to  an  east¬ 
erly  course.  These  data  were  collected  later  in  the  tidal  cycle  near  high  tide  when  the 
current  was  decreasing.  Figure  5.37  shows  that  the  area  with  high  signal  levels  to  the 
north  near  Decatur  Reef  splits  into  two  distinct  bands.  One  curls  back  to  the  west  around 
Restoration  Point  heading  southwest,  and  the  other  heads  out  into  the  middle  of  Puget 
Sound  to  the  southeast. 

Results  from  the  second  day  of  towing  on  1 1  June  are  shown  in  Figure  5.38.  This 
day  consisted  of  repeated  transects  both  east  and  west  along  47°34'40"'  N  latitude.  The 
most  active  area  during  these  transects  has  an  rms  level  of  1.1  x  10“^  s“^  (over  a  20-s 
interval).  Figure  5.39  shows  vorticity  variance  for  20-s  intervals  for  data  from  1  July. 
The  highest  vorticity  levels  are  southeast  of  Restoration  Point  during  the  peak  of  the 
flood  flow.  Although  closer  to  Decatur  Reef,  the  presumed  source,  the  area  just  south  of 
the  point  shows  lower  vorticity  signals.  We  believe  this  is  because  the  mean  flow  is 
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Figure  5.37.  Vertical  vorticity  data  from  the  second  half  of  10  July  1991.  Note  the  two  bands 
of  high  vorticity  variance:  one  curling  back  to  the  west  around  Restoration  Point,  and  one  head¬ 
ing  out  into  the  middle  of  Puget  Sound. 
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Figure  5.38.  Time  series  of  vertical  vorticity  data  from  11  July  1991.  The  lower  panel  shows 
vorticity  corrected  for  yaw  rate  and  bandpass  filtered  as  a  function  of  time.  The  upper  panel 
shows  the  variability  in  rms  vorticity  computed  over  20-s  windows. 
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Figure  5.39.  Vertical  vorticity  variance  at  20-s  intervals  for  data  from  11  July  1991.  The  size 
of  the  dots  corresponds  to  vorticity  variance  and  the  average  latitude  and  longitude  of  the  20-s 
interval  determines  the  center  of  the  dot. 

parallel  with  the  shoreline  to  the  north  of  the  point,  causing  the  vorticity  signals  to  be 
advected  the  southeast.  The  water  velocities  measured  with  the  drogue  buoys  the  previ¬ 
ous  day  near  Restoration  Point  indicated  just  this  type  of  water  motion. 

Figure  5.40  shows  nine  along-track  plots  of  vertical  vorticity  taken  as  we  passed 
back  and  forth  along  47°34"40""  N  just  to  the  south  of  Restoration  Point.  The  zero  for  the 
along-track  axis  is  a  point  due  south  of  Restoration  Point.  The  earliest  transect  is  shown 
at  the  top  of  the  figure  and  the  latest  at  the  bottom.  The  start  and  stop  times  for  each  pass 
are  printed  on  the  plot.  Figure  5.41  shows  nine  along-track  plots  of  rms  vorticity  calcu¬ 
lated  for  20-s  intervals.  Again,  the  earliest  transect  is  on  the  top.  It  can  be  seen  from 
these  plots  that  the  area  containing  the  highest  vorticity  signals  shifts  back  and  forth  dur¬ 
ing  the  observing  period. 

Figure  5.42  shows  the  data  from  pass  1  (the  top  panel  in  Figure  5.41)  formed  into 
100-m  averages;  that  is,  all  the  data  within  a  100-m  band  have  been  averaged  together 
and  displayed  as  a  single  point.  The  pass  has  been  broken  into  26  segments,  and  the  aver¬ 
age  is  plotted  at  the  geometric  center  of  each  segment.  (The  first  segment,  for  example. 
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Figure  5.40.  Nine  along-track  plots  of  vertical  vorticity  along  a  line  just  south  of  Restoration 
Point.  The  earliest  transect  is  shown  at  the  top,  and  the  latest  at  the  bottom;  start  and  stop  times 
are  noted  on  each  plot. 
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Figure  5.42.  Data  from  pass  1  (top  panel.  Figure  5.39)  formed  into  100-m  along-track  averages 
and  displayed  as  a  single  point. 


runs  from  -550  m  to  -450  m,  and  the  average  is  plotted  at  -500  m.)  When  this  is  done 
for  all  nine  passes,  the  data  for  a  particular  segment  can  be  plotted  versus  time  (pass 
number).  Figure  5.43  shows  this  type  of  time  series  for  the  1000-m  and  1 100-m  seg¬ 
ments.  The  peak  signal  levels  are  observed  during  pass  6,  whose  median  time  is  about 
2300  UTC.  This  corresponds  to  1500  PST,  near  the  calculated  time  of  maximum  flow 
due  to  tidal  forcing  (see  Figure  5.34).  If  we  combine  the  data  from  all  the  100-m  seg¬ 
ments  and  all  nine  passes,  we  can  plot  a  three-dimensional  map  of  vorticity  signal  levels 
relative  to  distance  along-track  and  time.  Figure  5.44  shows  contours  of  equal  rms  vorti¬ 
city  signal  levels  versus  distance  (horizontal  axis)  and  time  (vertical  axis).  During  passes 
1, 2,  and  3,  the  rms  vorticity  levels  were  high  over  a  large  portion  of  the  tow  track.  Later 
during  passes  8  and  9,  the  high  vorticity  signals  were  localized  to  a  small  area.  The  peak 
levels  shifted  to  the  west  as  time  progressed. 

In  conclusion,  the  EMVM  tows  at  Restoration  Point  and  Decator  Reef  found  clear 
evidence  of  environmental  signals  with  rms  amplitudes  as  great  as  1.1  x  10"^  s"’  over  a 
20-s  interval.  The  locations  of  these  large  signals  were  consistent  with  Decatur  Reef  as 
the  source  and  advection  at  the  water  velocities  measured  with  the  drogue  buoys. 
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5.4  Toliva  Shoal 

On  27-29  August  1991,  the  EMVM  was  deployed  around  Toliva  Shoal  in  South 
Puget  Sound.  The  instrument  utilized  the  GW  fiberglass  fairing  and  was  towed  behind 
RA^  Miller  at  low  speeds  (2.5  kn)  and  at  a  depth  of  about  30  ft.  Toliva  Shoal  in  the  mid¬ 
dle  of  Puget  Sound  just  southwest  of  the  Tacoma  Narrows.  Figure  5.1  shows  the  operat¬ 
ing  area  in  relation  to  the  rest  of  Puget  Sound.  This  underwater  feature  rises  from  a  mean 
depth  of  400  ft  to  within  20  ft  of  the  surface.  The  center  of  the  shoal  is  marked  at  the  sur¬ 
face  with  a  buoy.  The  tows  were  conducted  during  an  ebb  flow  while  circling  the  center 
of  the  shoal  at  a  distance  of  about  0.25  n.mi.  Figure  5.45  shows  the  track  of  the  towing 
vessel.  During  these  tows,  a  vehicle  dynamics  package  was  attached  to  the  instrument. 
In  addition,  an  Eagle  4  GPS  receiver  and  two  Megapulse  4(XX)  LORAN  receivers  were 
installed  on  board  RA^  Miller,  and  position  data  were  logged  at  a  1-Hz  rate. 

The  permanent  magnet  (1.4x10“^  T)  was  again  used.  The  towing  configuration 
was  the  same  as  the  previous  deployments  in  Puget  Sound  in  July. 

Figure  5.46  shows  a  time  series  of  the  vertical  vorticity  data  for  28  August.  The 
lower  panel  depicts  the  rms  vorticity  corrected  for  yaw  rate  and  bandpass  filtered  as  a 
function  of  time.  The  upper  panel  shows  the  variance  in  vorticity  computed  over  20-s 
windows.  During  the  least  active  portions  of  the  tow,  the  vorticity  variance  was  slightly 
less  than  1.5  x  lO"^  s"^  (rms  =  4  x  10“^  s“^).  In  the  most  active  portions,  the  vorticity 
variance  was  1.2  x  10~^  s~^.  The  time  series  is  better  understood  by  breaking  it  up  into 
four  sections,  each  corresponding  to  a  "circle"  of  data  acquired  as  we  towed  the  instru¬ 
ment  around  Toliva  Shoal. 

Figures  5.47  through  5.50  show  the  vertical  vorticity  variance  plotted  as  variable¬ 
sized  dots  whose  centers  are  plotted  at  the  latitude  and  longitude  where  the  data  were 
acquired.  On  the  first  two  passes  (Figures  5.47  and  5.48),  we  steered  a  constantly  chang¬ 
ing  heading  which  produced  a  circular  path  with  a  radius  of  0.25  n.mi.  On  the  last  two 
passes  (Figures  5.49  and  5.50),  we  were  trying  to  produce  a  box  pattern  with  sides 
approximately  0.5  n.mi.  long  centered  on  the  submerged  shoal.  The  fourth  circle  was 
interrupted  when  the  submerged  sensor  hit  an  unknown  object  and  all  five  ceramic  sensor 
tubes  were  sheared  off  at  their  bases.  The  data  show  both  temporal  and  spatial  variabil¬ 
ity;  in  general,  however,  the  geographical  location  with  the  highest  vorticity  signal  levels 
is  to  the  northeast  of  the  buoy  marking  the  center  of  the  shoal.  This  is  consistent  with  the 
observing  period  being  during  ebb  flow,  when  the  bulk  of  the  water  that  was  flowing  past 
(and  over)  Toliva  Shoal  was  progressing  up  the  main  channel  of  Puget  Sound  just  east  of 
Fox  Island.  This  set  of  conditions  puts  the  area  with  highest  vorticity  signals  directly  in 
the  wake  of  the  shoal. 

To  investigate  the  azimuthal  distribution  of  vorticity  with  respect  to  the  center  of 
Toliva  Shoal,  the  data  for  the  first  three  circles  around  the  shoal  were  sorted  into  7.5° 
bins.  Using  coordinate  transformation  algorithms,  we  transformed  the  latitude/longitude 
data  into  radial  distance-vs-angle  bins.  The  data  were  then  sorted  into  7.5°  angle  bins 
and  averaged.  Because  the  vessel’s  speed  over  ground  was  not  constant  owing  to  wind 
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Viewpoint 


Figure  5.45.  Ship’s  track  around  Toliva  Shoal  in  Puget  Sound  in  August  1991.  The  ship  cir¬ 
cumnavigated  the  shoal  four  times.  Two  of  the  circumnavigations  were  circles  at  a  radius  of 
0.25  n.mi.,  and  two  were  boxes  with  sides  0.5  n.mi.  long  and  centered  at  0.25  n.mi.  from  the 
shoal. 
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and  currents,  the  number  of  data  points  in  a  given  angle  bin  varied  from  a  high  of  20  to  a 
low  of  5.  These  binned  and  averaged  vorticity  variance  data  are  shown  in  Figure  5.51. 
This  plot  confirms  our  hypothesis  that  the  area  with  the  highest  vorticity  signals  is  in  the 
wake  of  the  submerged  feature.  The  area  with  the  lowest  vorticity  was  in  the  relatively 
undisturbed  water  coming  from  the  south,  where  the  nearest  flow-disturbing  feature  (an 
island)  was  over  4  n.mi.  away. 

5.5  Horizontal 

Two  field  tests  were  conducted  in  Puget  Sound  in  which  the  vorticity  sensor  was 
operated  horizontally,  one  on  24  February  and  one  on  27  February  1992.  For  both  tests, 
the  vorticity  sensor  was  attached  to  a  fiberglass  tripod  sitting  on  the  bottom  at  a  depth  of 
15-20  ft.  The  EMVM  obtained  high-quality  data  during  both  deployments. 


Circle  1  Circle  2  Circle  3  Circle  4 


Figure  5.46.  Time  series  of  vertical  vorticity  data  for  28  August  1991.  The  lower  panel  shows 
vorticity  corrected  for  yaw  rate  and  bandpass  filtered  as  a  function  of  time.  The  upper  panel 
shows  the  variability  in  rms  vorticity  computed  over  20-s  windows.  The  circle  the  data  come 
from  is  noted  at  the  top. 
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Figure  5.47.  First  circle  of  Toliva  Shoal,  1800-1850  UTC.  The  size  of  the  dots  indicates  the 
magnitude  of  the  vertical  vorticity  variance  (s“^);  the  scale  of  the  dots  is  shown  in  the  lower  left- 
hand  comer. 
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Figure  5.48.  SeconiJ  circumnavigation  of  Toliva  Shoal,  1850-1940  UTC.  The  size  of  the  dots 
indicates  the  magnitude  of  the  vertical  vorticity  variance  (s"^);  the  scale  of  the  dots  is  shown  in 
the  lower  left-hand  comer  of  the  figure. 
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Figure  5.49.  Third  circumnavigation  of  Toliva  Shoal,  1940-2030  UTC.  The  size  of  the  dots 
indicates  the  magnitude  of  the  vertical  vorticity  variance  (s“^);  the  scale  of  the  dots  is  shown  in 
the  lower  left-hand  comer  of  the  figure. 
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Figure  5.50.  Fourth  circumnavigation  of  Toliva  Shoal,  2030-2106  UTC.  The  size  of  the  dots 
indicates  the  magnitude  of  the  vertical  vorticity  v^ance  (s“^);  the  scale  of  the  dots  is  shown  in 
the  lower  left-hand  comer  of  the  figure. 
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Figure  5.51.  Azimuthal  distribution  of  vertical  vorticity  variance  (s"^).  Vorticity  variance  dur¬ 
ing  the  first  three  circumnavigations  of  Toliva  Shoal  was  sorted  into  7.5°  bins  and  averaged.  The 
radial  distance  indicates  the  magnitude  of  the  vorticity  variance  according  to  the  scale  shown. 
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Small  signals  perturbing  the  measurements  were  noted  during  past  experiments 
when  the  EMVM  was  operated  from  a  steel  tripod.  A  new  tripod  was  constructed  espe¬ 
cially  for  these  tests.  It  was  made  entirely  of  fiberglass  and  plastic  so  that  electric  fields 
due  to  the  corrosion  of  dissimilar  metals  were  not  present.  The  new  tripod  allowed  the 
vorticity  sensor  to  pivot  freely  around  the  vertical  axis. 

During  the  experiments,  the  EMVM  was  operated  horizontally.  A  large  plastic 
vane  and  drag  parachute  attached  to  the  sensor  caused  it  to  pivot  about  the  vertical  axis 
and  remain  aligned  with  the  local  flow,  which  had  mean  values  of  up  to  0.4  m/s.  Both 
deployments  were  conducted  during  a  period  that  spanned  the  maximum  ebb  flow.  The 
sensor,  which  was  mounted  at  the  top  of  the  tripod,  was  submerged  only  3-5  m  below  the 
surface,  depending  on  the  height  of  the  tide. 

Several  interesting  features  were  observed.  During  both  experiments,  large  vorti¬ 
city  signals  occurred  when  there  were  only  small  velocity  signatures.  Figure  5.52  shows 
a  portion  of  the  data  recorded  during  the  second  deployment.  The  horizontal  vorticity 
time  series  contains  substantial  signals  with  rms  amplitudes  of  3.2  x  10“^  s~^.  There  are 
no  horizontal  or  vertical  velocity  signals  that  correlate  well  with  the  observed  vorticity. 
There  is  also  a  patch  of  water  near  the  middle  of  the  record  that  has  much  a  much  lower 
vorticity  variability,  less  than  5  x  10”^  s“’  rms.  There  is  no  apparent  change  in  either  the 
mean  velocity  or  the  rms  variability  of  velocity.  Likewise,  there  are  large  velocity  sig- 


Figure  5.52.  Data  acquired  with  the  EMVM  sensor  aligned  so  that  the  sensitive  axis  was  hor¬ 
izontal.  The  top  panel  shows  the  horizontal  velocity  component  measured  along  the  local  stream¬ 
lines.  The  middle  panel  shows  the  vertical  velocity  component.  The  bottom  panel  shows  the  hor¬ 
izontal  component  of  vorticity  aligned  perpendicular  to  the  local  streamlines  and  the  vertical. 
There  are  substantial  signals  in  vorticity  that  do  not  have  corresponding  signatures  in  velocity. 
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nals  due  to  variations  in  the  mean  flow  and  surface  waves  without  corresponding  vorti- 
city  signals.  During  the  second  experiment,  the  EMVM  measured  horizontal  velocities 
of  up  to  20  cm/s  peak  to  peak  due  to  surface  waves  from  a  large  vessel  passing  nearby. 
Vertical  velocities  of  up  to  10  cm/s  were  simultaneously  observed.  This  event  is  evident 
in  Figure  5.53.  The  middle  and  upper  panels  show,  respectively,  the  vertical  and  one 
component  of  the  horizontal  velocity  field.  The  bottom  panel  shows  the  horizontal  vorti- 
city.  The  surface  waves  peak  near  2315  UTC  and  appear  as  a  high  frequency  on  the 
record  because  of  the  time  scale  used  for  plotting.  Visually,  little  of  the  signal  appears 
present  in  the  vorticity  data.  This  is  the  expected  result.  A  more  sensitive  analysis  could 
be  done  using  spectral  techniques. 


920227-2310  920227-2315  920227-2320  920227-2325  920227-2330  920227-2335  920227-2340 

TIME,  s 

Figure  5.53.  Data  acquired  with  the  EMVM  sensor  aligned  so  that  the  sensitive  axis  was  hor¬ 
izontal.  The  uppermost  panel  shows  the  horizontal  velocity  component  measured  along  local 
flow  streamlines.  The  middle  panel  shows  the  vertical  velocity  component.  The  bottom  panel 
shows  the  vorticity  about  an  axis  perpendicular  to  the  local  streamlines  and  the  vertical.  The 
high-frequency  signals  present  in  the  first  half  of  the  record  are  due  to  the  waves  of  a  passing 
vessel.  These  surface  gravity  waves  do  not  contain  vorticity  and  therefore  are  not  present  in  the 
vorticity  data  channel. 
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